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1 Introduction 

Since many decades nanoporous materials, for example zeolites, play an eminently im- 
portant role in the catalysis of oil refining and petrochemistry. On the other hand, 
molecular sieve materials began also to attract some attention as optical material in 
last years. It was realized that their nanometer size pores allow to host guest molecules 
giving so substance to a new class of optical material with properties which neither the 
host, nor the guest alone could ever possess. In this way new pigments and luminophores 
were realized as well as novel optically nonlinear and switching materials. The various 
actually realized materials are reviewed in this book in chapter Nanoporous compound 
materials for optical applications - Material design and properties. 

A closer look at the approach of arranging molecules in an ordered framework of 
pores reveals a series of aspects of fundamental interest: For example, the stereometric 
restrictions which the pore framework imposes on the motional degrees of freedom of 
the enclosed molecules reduces their diffusion to one dimensional diffusion in channel 
pores T]. Or, as with the concentration also the average distance between two guest 
molecules is controlled, their dipolar near field interaction can so be adjusted. In opti- 
mal circumstances optical excitation energy can then be transferred nonradiatively over 
distances of several micrometers [2j. That are just two examples of new phenomena in 
molecular sieves, which at this moment are still studied to achieve full understanding, 
and which will soon find their way into applications in science and technology. 

In this chapter we will be concerned with another new optical application of nano- 
porous compound materials, namely microlasers in which light is generated by organic 
dyes embedded in wavelength size resonators of molecular sieve material. In these laser 
devices the dye molecules are enclosed in nanometer size channel pores of the molecular 
sieve, whereas the crystalline sieve material itself which forms the resonator has exterior 
dimensions on the order of a few micrometers. Figure Q illustrates the arrangement 
of the molecular dye dipoles which are aligned in the pores of a hexagonally shaped 
AIPO4-5 molecular sieve host. In this example of molecular sieve material the channel 
pores point to the direction of the crystal c-axis. If the enclosed species of dye molecules 
have an elongated shape and exhibit a transition dipole moment along their elongation 
axis, then the dipoles end up oriented parallel to the host c-axis as well. 

A dipole can not emit in direction of its axis. Therefore the emission of the dipoles 
shown in Fig. ^ occurs in a plane normal to the host crystal c-axis. Once the light 
emitted by the dipoles arrives at the boundary of the molecular sieve crystal, part of it 
is reflected back into the material. In fact, given the hexagonal geometry there is even 
a bundle of directions for which total internal reflection occurs. Figure |21 illustrates how 
optical rays of this bundle loop around the crystal, and, reflected by the hexagonal sides, 
form a whispering gallery mode. In that way regularly shaped crystals of molecular sieve 
hosts form microresonator environments for the light emission of inclosed dye molecule 
guests, and if the molecules provide sufficient optical gain, laser oscillations will build 
up0. 

Light emission in a microcavity environment, however, is in many respect different 
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Figure 1: Schematic picture of the arrangement of the molecular dye dipoles in the chan- 
nel pores of the hexagonal molecular sieve host crystal AIPO4-5. The spatial 
cosine-square emission efficiency of a dipole is represented as doughnut shaped 
surface around the dipole axes. 
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Figure 2: Light which is confined by total internal refiection at the hexagonal boundaries 
can circulate inside the molecular sieve crystal as whispering gallery mode. 



from the familiar emission into free space, such as fiuorescence. In fact, this difference 
can be striking. For example, to achieve laser action in a conventional millimeter size 
(or larger) laser resonator the pump must overcome a certain threshold. In a microlaser 
in which the resonator is of the size of a few wavelengths, however, lasing can occur 
without threshold. That means that every absorbed pump quantum is transformed into 
a laser photon. Therefore, in order to appreciate the significance of light generation in 
molecular sieve based microlasers, we need to understand the differences between the 
emission processes of a molecule in a free space environment as opposed to emission in 
a micro cavity environment. 

In the following we will review spontaneous as well as stimulated light emission of 



molecules, and we will show that the respective emission rate is not an inherent molecular 
property, but is a function of its environment, or more precisely, of the mode density 
of the electromagnetic field. This becomes apparent when in the distance of a half to a 
few wavelengths of the molecule mirrors exist. But this is indeed exactly the situation 
of here discussed dye molecules which are enclosed in a molecular sieve microresonator. 

On the other hand, we know that emission and absorption of radiation is accompanied 
by a transition of the molecule from a state with energy Ei to a state with energy E2. 
The frequency of the emitted or absorbed light is then u = \Ei — E2\/h, where h is 
Planck's constant. The presence of h clearly shows that emission and absorption of light 
is intrinsically a quantum mechanical process. Therefore our discussion will have to 
involve the quantum aspects of the interaction of the dye molecules with the light field, 
as well as the quantum nature of the light field itself. 

The situation is even more peculiar because the resonators of the molecular sieve 
microlasers have a hexagonal outline. For a hexagonal resonator one can not find an 
orthogonal coordinate system in which the wave equation can be solved by the usual 
method of separation of variables. Thus we will have to discuss the properties of mi- 
croresonators in view of this impediment. 

In consideration of these facts we have organized the discussion in the following way: 
In the first section we introduce the concept of modes of the electromagnetic field as its 
countable degrees of freedom, and based on this, we introduce the quantized optical field. 
In the next section fluorescence, i.e. spontaneous emission in a free space environment 
is discussed, and the frame is set for the treatment of cavity effects in the next section, 
in which spontaneous emission in a resonator environment is examined. Then the effect 
of a resonator on stimulated emission and laser action are surveyed. After this we 
characterize the peculiarities of microresonators, and finally we present the most recent 
achievements and realizations of molecular sieve microlasers. 

2 The concept of "modes" of the electromagnetic field 
and its quantization 

In optics in general, and particularly when lasers are involved, the notion of modes is 
ubiquitous. As this term is used in many different and disparate circumstances it is 
necessary to define the term for further usage. In this section we give a short tutorial 
introduction of the concept of m,odes of the electromagnetic field, and we show the im- 
portant role the mode concept plays in the procedure of the canonical quantization of the 
field. An account of the mode structure of resonators and particularly microresonators 
is then presented in section IHl 

2.1 The dynamics of the classical field 

A convenient way to introduce the mode concept is to consider a simple realization of 
the electromagnetic field, for example a source-free field. We can think of it as the field 



that subsists after a source located far away has stopped to emit. After the emission 
stopped, the classical field evolution is governed by Maxwell's equation in the following 
form (cgs-units): 



V-E = 


V X E = 


IdB 
c dt 


V-B = 


V xB = 


IdE 
c dt 



(1) 



2.2 Discretizing fields - Random fields 



Most practical light sources, such as incandescent lamps or light emitting diodes LEDs 
(though not lasers), emit a nondeterministic, chaotic field, i.e. a field whose spatio- 
temporal evolution can only be described in statistical terms. In communication engi- 
neering terms such a field is referred to as a noise field. We like to point out that we can 
understand important features of those fields in classical terms. It is thus not necessary 
to enter the quantum world in order to encounter nondeterministic fields. In order to 
keep the mathematics simple, it is convenient to deal with a field occupying a discrete, 
instead of a continuous number of degrees of freedom. 

The trick usually used to achieve discretization consists in confining the field into a 
finite volume V. At the end, the limit V^ — >■ cx) can be carried out, if necessary. Con- 
finement to a finite volume V allows us to represent the spatial component of the field 
as a series (superposition) of a discrete number of functions. As (Hj) is a system of hnear 
equations of the fields E and B we can always represent a solution as a linear superposi- 
tion of field functions. Of course it is convenient to choose a complete and orthonormal 
set of functions {<^i(r)}, that also fulfill eventually given boundary conditions. These 
functions are called (spatial) modes. Thus a mode of the classical electromagnetic field 
is characterized by the following properties: 



{ipi{r)} orthonormal: 



ip';^ir)ip4r)dV = 5^n (2) 

V 



{ipi{r)} complete: 



E.(f,t) = ^C,v^,(f)e-'-'* + 5^C;^:(r)e^-»* (3) 

i i 

E.{f,t) = E(+) + E(-) (4) 

{{pi{r)} satisfies the spatial boundary conditions given by the shape of volume V. 



Obviously Ex = (^Ex )* (with * we denote complex conjugation). The sets {Ci}x 
are discrete, and they now represent the complete information about the field. In a 



y,z 



deterministic field the individual {Ci}x,y,z are fixed complex numbers. For a nondeter- 
ministic field, however, the {Ci}x,y,z represent random variables. Thus for a stationary 
field they are defined in terms of probability functions: 

pi{a}) = piCuC2,Cs, ...) (5) 

If we consider a function F that depends on the random field, say E, or E^^\ then 
we can only express F in statistical terms, that means we can only assign expectation 
values: 

{F{E(^^))= f p{{C,})F[EM({C,})]l[d'C, (6) 

^^ k 

In contrast to a thermodynamic situation in which the expectation values are sharply 
peaked, an experimental realization of F in optics can significantly differ from the ex- 
pectation value {F). 

2.3 The classical Hamiltonian of the source-free field 

For a source-free field the classical Hamiltonian H can be interpreted as the total energy 
of the field. As we have constrained the field to the Volume V, the total energy is given 
byi: 

H = l [ [soE^f, t) + -B\r, t)] dV . (7) 

^ Jv /^o 

In the following we will show that the Hamiltonian ((Tj) can be represented as a sum of 
terms that are analogous to a harmonic oscillator (The thoughtful reader anticipates the 
reason...). For that purpose we express the fields E and B in terms of their potential. 

2.3.1 The potential of the free field 

We recall that Maxwell's equations of a source free field are gauge invariant. In the 
case we discuss here we choose the Coulomb gauge, and as a result we obtain a purely 
transverse field potential A: 

V ■ 1 = . (8) 

The electric and magnetic fields are related to the potential according to: 

E{f,t) = --A{r,t) , Bir,t) = V X A{f,t) . (9) 



Inserting this into (P) we obtain the wave equation 



V^l-i|^i=0 . (10) 



2.3.2 Discretization procedure for the potential A 

We follow the thread outlined in section 12.21 and to keep things simple, we consider a 
cube shaped volume V with an edge length of L. With this conditions the discretization 
of A acquires the form of a Fourier series (plane wave expansion) : 

i(r,t) = -=^5^^g(t)exp(i^-f) , (11) 

k 

where the wave vector k has the components 

ki = 27mi/L , rii = 0, ±1, ±2, ... 

k2 = 21:712/1 , ^2 = 0, ±1, ±2, ... (12) 

k^ = 27171?,/ L , ng = 0, ±1, ±2, . . . , 

and ^^ extends over the modes indexed by ni,n2,n3. The chosen normalizing factor 
will soon prove to be useful. Evaluation of the gauge relation (jH)) results in 

Y^ k ■ Af^it) exp(iA^ ■ r) = (13) 



k 

for all r. This is only possible when 

k ■ AS) = , (14) 

thus k±Aj: (i.e. transversal field). As the potential A{r,t) assumes real values, the 
coefficients Aj^ observe 

A_S) = 4(t) . (15) 

In addition, wave equation (fTUj) must be satisfied, resulting in 

^ + -1) Mt) = o, (16) 

with ujj: = ck. The general solution of this ordinary differential equation is represented 
as 

1^)= &,e'-^-n^ + c:^^-^'. (17) 

As ()14|) fixes the transversal character of each plane wave mode, only two components 
(the polarization components) of the vector Aj:{t) are at free disposition. In order to 
simplify the notation let us agree to let index ^ point to ni,n2,n3 (see dH))), as well as 
to the polarization components ^1,^2- To remember this, we will from now on refer to 
the corresponding index as k^- With this notation we can write for the potential 

Ar, t) = --^ Y.l^^.,e--^' + c^^e'^-^^Y'-' , (18) 



and with 

«k,(t) = Ck,e-'"'^«* (19) 

we get 

Ar,t) = -|^EKWe''''><,We"^'1- (20) 



By inserting ()2()|1 in 0, we can express the electric and magnetic field as a sum of 
mode functions 

mt) = -=^$^K(t)(^xeV'^-^-c.c.] (22) 



where with ^ we represent the polarization unit-vector. We note that the total informa- 
tion on the (classical) field is now contained in the functions Ukt(^) (jl^p . 

Note: To avoid confusion, we have to point to a minor inconsistency in the notation: 
As is typical for expressions of the magnetic field, the vector product in (|22|) reshuffles 
spatial and polarization indices so that the correct expression in fact consists of two sums 
containing the polarization vector. The resulting expression looks bulky, and requires 
some consideration. For this tutorial we prefer to emphasize the basic mathematical 
structure. So we choose this visually intuitive representation although the indices are 
not correctly rendered. To obtain the correct result one may work out the procedure in 
the component notation of (p]j) and Q. 

2.3.3 The mode density 

Let us return for a moment to the volume L^ we considered for deriving the mode expan- 
sion (fTT|). which was a cube with edges oriented along the coordinate axes (xi, 0:2,3:3). 
According to (jll|) we can expand the field in this cube in a 3-dimensional set of run- 
ning modes {ki}, {k2}, {^3} [cf. (fT^ ]. Along the Xi-axis we have the modes labeled by 
ki = 27mi/L] {rii = ±1, ±2, . . . ); because we consider only running modes (travelling 
waves), we omit rii = 0. In the interval between ki and ki + dki we find drii modes, 
where drii = dki L/2'k. The analogous applies for the other directions. According to 
Fig.iniwe can specify the number of modes in the interval between \k\ and |/i; + d/c| simply 
by counting the dots in the volume of the corresponding spherical shell, which amounts 
to Ank^ dk ■ 2 (factor 2 because there are two independent polarization states associated 
with each k) . Considering the scaling of the axes in Fig. El we can use this to express 
the number of modes in Volume V^ = L^ as 

Number of modes in \^ = 2 ■ Ank'^dk ( — | . (23) 
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Figure 3: Mode distribution in /c-space scaled by multiplication with ^. Each dot 

corresponds to a wavevector k which represents two independent, orthogonally 
polarized waves. 

The mode density p{k)dk, that is the number of modes in interval |dA;| per Volume 
V = L^, is then given by 



p{k)dk 



# of modes in V k"^ dk 



V 



TT^ 



With A; 



we obtain 



p{uj)duj 



00 



^j^2 



duo . 



(24) 



(25) 



This is known as the free space mode density. Note that p{uo) represents the volume 
normalized density of modes, i.e. the units of p{uo) are number of modes / (frequency x 
volume). 

Let us now consider a linear function C, of the field, which thus can be represented as 
a sum over the field modes k as 



c 






(26) 



For a large cube, L -^ oo, the sum over the discrete set of modes has to be transformed 
into an integral, which has to sum over a measure that is a density 



k E «*) 



L3 



ki,k2,k3 



(2vr) 



d^'A; C{k) . 



(27) 



Switching to spherical coordinates d^k = k'^dk sin^d^ d0, where denotes the az- 
imuthal and 6 the polar angle of k, and transforming this equation to frequency space 
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uj we obtain observing k = uj/c 



ki,k2,k3 

In this representation we take into account that ( can depend on the polarization, 
which can be expressed as C = C(0; ^; k). Also, C allows us to introduce fc-space structure 
functions into the mode counting sum, so that the mode density can be calculated for 
arbitrary boundary conditions. Note that the procedure relies on certain properties of the 
field, one of which is (|7j). This equation expresses energy conservation, thus the above 
calculated mode density refers to undamped modes. The situation with microlasers, 
where sources of the field are present and the modes are damped, is more intricate and 
is discussed in section |H] 

2.3.4 The classical Hamiltonian of the source free field 

Let us insert the electric field (jTIj) and the magnetic field (|^ into the energy ((Tj) (note 
that J^s e^^'^~'^'''^d^r = L^ Sj:g), then we obtain the Hamiltonian 

H = 2J2^'\^^M'- (29) 

The energy appears here as the sum of the energy of each mode. This is intuitive, and 
at this point we could stop satisfied with the result. However, there is one thing that 
experience revealed: As elegant it might be, expression (j^IH) does not lend itself to an 
idea that opens a viable way for the quantization of the electromagnetic field. 

Although at this time it looks utterly artificial, but the following substitution proved 
to be a wonderfully prolific device: 

gk,(t) = K(t) + <,(t)] (30) 

Pke(i) = -iujk^[uk^{t) - ul^{t)] (31) 

At this time this expressions can be regarded as the definition for the functions q and 
p that, as can be shown, satisfy the relations gk^ = pk^ and p^^ = —u'^q^^. Thus the 
functions q and p evolve in analogy to the canonical position and momentum variables of 
the harmonic oscillator: In fact, inserting (J3U|) and (j31|) in (j29|) . we obtain the expression 

which has the same form as the Hamiltonian of a series of independent (not coupled) 
harmonic oscillators. 

What do we have achieved so far? 

1. By restricting the fields to a finite Volume V, the potential of the electromagnetic 
field was decomposed in a discrete, although infinite, series of spatial modes. 
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2. Each spatial mode was associated with a harmonic oscillator. 

3. That this is a sound result can be seen by verifying the canonical equations of 
motion, which indeed show that dH/dpk^ = q^^ and dH/dq^^ = — Pk^. 

To be complete we write the fields expressed in terms of q and p : 

A{f,t) = _^5^{[g,^(t)+^p,^(t)]e'^^-+c.c.} (33) 

mt) = -^5^{Kgk,(t)+ipk,(t)]e^^^--c.c.| (34) 

B{r,t) = -^E{h(t)+::^Pk,(t)](^xfe^^>)-c.c.} (35) 



Note: As in (j2!2j) the same inconsistency arises here with (j35|) . See the note on p.llUI 

Thinking about quantization of the electromagnetic field? Well, quantizing a harmonic 
oscillator - we know how to do this, don't we? 

2.4 Canonical quantization of the electromagnetic field 

We assume that the reader is familiar with the postulates of quantum mechanics in 
general, and with the quantum mechanics of the harmonic oscillator in particular (see for 
example j3]). In anticipation of the following procedures, we have already put together 
the main constituents we need to apply the correspondence principle 5, p. 337ff] for 
quantizing the field. As first step we associate with each of the classical dynamic variables 
a Hilbert-space operator. At this time we choose operators in the Heisenberg picture. 
In this way the correspondence to the classical formulas is more obvious. As always in 
quantum mechanics we also have to carefully consider commutation of the operators. 

According to the correspondence principle intent jHl IZj, we associate the following 
operators with the classical canonical variables (|3U|) and (|31|): 

gk,(t) -^ qk,(t) and (36) 

Pke(t) -^ Pk,(t) • (37) 

The operators q and p of the same mode represent noncompatible operators. We use 
that the commutator of a canonically conjugated pair of operators amounts to ih. On 
the other hand, equation (jH^ indicates that the modes (oscillators) are uncoupled. 
Therefore the associated Hilbert-space operators of different modes will commute. Thus 
we obtain the following commutation relations: 

(38) 
(39) 
(40) 
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[qkeW,Pk^W] 


= i^^kek^ 


qkeW.cik^W^ 


= 


[Pk,(i),Pk^(i)] 


= , 



and accordingly, the following uncertainty relation: 

Agk,Ap^ > ||([qk„Pk^])| (41) 

h 

2 "k^k^ 



> fCki • (42) 



The state of a quantum mechanical system - here the electromagnetic field - is char- 
acterized by a state vector, say the ket \ip). The measurement of the observable O will 
produce a value coinciding with an eigenvalue of the Hilbert space operator associated 
with the observable O. However, only when the state of the field \ip) happens to be an 
eigenstate of the operator 0, we can precisely predict the outcome of the measurement 
of O. If the state of the system does not coincide with an eigenstate of 0, then it is only 
possible to predict the probability to obtain a certain eigenvalue. Thus the outcome 
of the measurement of O is predicted by the probability distribution given through the 
scalar product {ip\0\ip). 

At this point the profit we gained by reinterpreting the classical Hamiltonian ()29j) in 
the form of ()H2|1 becomes apparent. By identifying the classical observables in ()H2|1 with 
the operators ()37|) and (jHTj) . we obtain a Hamiltonian operator that exhibits required 
quantum properties. Thus 

Note that with the transposition of the classical problem into Hilbert space, the variable 
space was transposed as well: The new dynamic variables now are the operators qk£(^)) 

Pk,(t), as well as the field operators A(r, t), E(r, t), B(r, t), etc.. The classical variables 
r and t are now relegated to play the role of mere parameters. 

In the following we try to find out why the analogous transposition of the classical 
Hamiltonian (j^Uj) fails to produce a Hamiltonian operator that is consistent with the 
quantum mechanical postulates. We will start to seek the quantum mechanical equiv- 
alent of mode functions u\^^{t) appearing in (^01) and (jHH). Obviously these must be 
operators. 

2.5 Creation and annihilation operators 

The operators that play the analogous role to the functions 'Ukf(t) in (jSHI) and (|^ 
are called annihilation operator aj^ (t) and creation operator ak, (t)^. They can be ex- 
pressed in terms of the Hermitian operators q k, {t) and p^ (t) as the following pair (the 
normalization will be justified later): 

<At) = ^==Kqk,(t)-iPk,(i)] • (45) 



ak,(t) = yi^Kqk,(t)+iPk,(t)] (44) 



^yrau^ 



^Although they play an analogous role, they do not correspond to the classical mode functions in the 
sense of the correspondence principle, as we will soon discover. 
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The comparison with ()H()|1 and (jHlll also shows that these operators exhibit the the time 
dependence 

ak,(t) = a^(0)e-K* (46) 

at,W = ai^(0)eK* . (47) 

Similarly as the classically equivalent functions in (jHUj) and (j^ are complex conju- 
gates, the annihilation and creation operators are Hermitian adjoints, which is expressed 
by the symbol ^. Since aj^ (t) 7^ ak^(t), they are not themselves Hermitian but Hermitian 
adjoints, and thus can not qualify as physical observables. Like Qk^l^) and Pk^it) they 
are not compatible (do not commute). From ()38|1 and ()4()j) we obtain the commutation 
relations 

[ak,(t),ai,(t)] = Sl^^, , (48) 

[ake(t),a^(t)] = , (49) 

[aI,W,aj,,(t)] = . (50) 

However the dynamic variables Qkfl*^) and Pk, (^) are Hermitian, and can consequently 
be associated with observables 



qk,(t) = ^/l:JJ[a^(t)+al,^(t)] , (51) 



Pk,(t) = iV^k,(t)-ak,(t)] . (52) 

After comparison of these two equations with (IHUj) and (jHH), we are confident that 
akj(t) and aj. (t) play an analogous role as the classical mode functions Mk^(t) and u^^ (t). 
Therefore, it is the set of operators {akj(t)} that in the quantum analogy carries the 
total information about the field. As the information of each quantized field mode is 
now expressed in terms of an operator, instead of simple c-numbers, we must expect a 
correspondingly richer structure of the quantized field - richness in the sense that in the 
quantum description, non classical field configurations arise that do not have a classical 
equivalent. 

We can now express the Hamiltonian operator in terms of annihilation and creation 
operators. We insert (|HT|) and (|H^ into PH|1 and obtain 

H = i5^;z^kjak,(t)ai^(t) + ai^(t)ak,(t)] (53) 

In this equation the annihilation and the creation operators appear in a symmetric 
way. From its visual appearance (j3Hjl can be considered to closely resemble the classical 
Hamiltonian (|^^. An the other side we have to consider that annihilation and creation 
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operators do not commute, and the apparent symmetry may only be a typographical 
one. There is no doubt that the Hamiltonian operator per se is Hermitian and thus 
represents an observable, the energy of the field. The energy of an optical field is usually 
measured by absorbing the light in a photo detector, say for example a photo diode. As 
will become obvious below [c.f. (J58|) - (J6(J|) ] for the description of absorption processes it is 
practical to have the Hamiltonian operator in a form in which the annihilation operators 
stand to the right of the creation operators. This normal ordering form is achieved by 
successive application of the commutation relations ()48 p - (j3m) . As a result we obtain 

H = J];ia;kJaL^(t)a^(t) + l] . (54) 



Obviously the operator product a|^ (t)akj(t) is Hermitian, and represents the number 
operator Nk^ 

Nke = ai^(t)ak,(t) . (55) 

As its name indicates, Nk^ counts the number of photons in each mode k^. Note that 
because of (j^6|) and ()47|1 the number operator Nk^ is constant in time - a mandatory 
property for a conserved quantity of the field, like H ()54|1 . The eigenvalues of the 
Hamiltonian operator (jK^ are hwi^^ x (ra^ + i), with nk^ = 0, 1, 2, 3, . . . , oo. It is evident 
that rikj represents the photon occupation number of mode kg. In (j3^ the contribution 
of ^/icukj to each mode represents the vacuum fluctuations of the field. Clearly, this term 
is not present in the classical Hamiltonian ((221) • 

At this point we should have discussed the eigenfunctions of the annihilation and 
creation operators. We postpone the discussion to the following section, in which 
we introduce the coherent states, and were we are able to better illustrate the special 
meaning of those eigenstates. 

Summarizing we can now reconstruct the reasons why the direct application of the 
correspondence principle to the classical Hamiltonian expression ()29|) fails, but works on 
(El: 

• The correspondence principle refers to physical observables. On the other side, 
from the perspective of classical physics, the concept of observables is not well 
defined, and therefore it is not possible to anticipate that the operators a and a^ 
associated with the classical mode functions Uk^^{t) and m^ (t) are not Hermitian 
and thus do not represent observables. 

• In the representations of the classical Hamiltonian the mode functions u^ (t) and 
u^ (t) enter in a symmetric way. The normal ordering procedure of the operators a 
and a^ in the Hamiltonian operator shows, however, that the equivalent operators 
do not contribute equally to the Hamiltonian. In ()54|) the factor ^/icukf representing 
the vacuum fiuctuations is a result of this asymmetry. Since vacuum fiuctuations 
do not exist in the classical picture, the term could not be anticipated. 
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2.6 States of the optical field 

In classical physics the information of the electrodynamic field is contained in the electric 
and magnetic fields, or their associated potential functions. In the quantum description, 
however, the information is contained in a state vector (wave function). In the following 
we discuss two important states of the field, the Fock states (or number states, states of 
fixed photon number), and the coherent states (harmonically oscillating fields). 

A practical way to represent a state vector is in terms of a series of eigenstates of a 
suitable operator, and this is what we are going to do as next. 

2.6.1 The Fock states (number states) 

The Fock states are defined as states with a fixed photon number. Thus they can be 
represented as eigenstates of the number operator ()55|) . Considering the mode k^, we 
have the following eigenvalue equation 

NkJ^^kj) = aj,^(t)ak^(t)|nk^) = rikj^k^) , (56) 

in which nk^ represents the number of photons in mode k^. The Fock states are stationary 
because of PUJ) and (jlTj) . In fact, we can see that a ground state |0) (vacuum state) is 
associated to each mode by 

ake|0) = . (57) 

The following relations motivate the names of the annihilation operator a and creation 
operator a^ : 



akjj^kj = y^ \nk^-l) (58) 

atj'^ke) = y^^k^n l^kj + 1) (59) 

aLakJ^k.) = rik, |^k.) (60) 



We can construct an arbitrary Fock state l^k^) by an rik^-fold application of the creation 
operator to the ground state. 



it \"k 



^^) = ^i=^(aLr^|0) . (61) 



/"k, 



The Fock states are orthogonal, and with the factor introduced in (jUj) and P^j) they 
are normalized, 

(nkjmkj) = (5„^ , (62) 

and complete 

oo 

Xll^kJ(^kel = 1 • (63) 

They therefore form a complete system of base vectors. The Fock state base is often 
used to describe fields with few photons of high energy, for example 7-radiation. Optical 
fields of visible radiation, like fields emitted by a laser, are more suitably described by 
coherent states. 
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2.6.2 Coherent states 

The coherent state represents the quantum mechanical analogue of a classical, harmon- 
ically oscillating field. To reduce clutter in the following discussion we will pick out 
one spatial mode of the field: |a)k£ -^ \(^)- In the following we adhere to Glauber's 
"classical" presentation Sj. The coherent state |a) of a mode is defined as eigenstate of 
the annihilation operator a: 

ala) = a\a) (64) 

with the complex eigenvalue a (remember, a is not Hermitian) 

a = \a\ e^'f' (65) 

a being not Hermitian, it is not surprising that the associated eigenstates, the coherent 
states I a), are not orthogonal. Thus they can not provide a universally suitable base 
system. Although coherent states can therefore not be used to expand an arbitrary 
electromagnetic field in the conventional way (see below), they represent an important 
class of fields, namely harmonically oscillating fields, like radio frequency fields, or fields 
emitted by well stabilized lasers, and for this reason are interesting to characterize. 

Fock states as eigenstates of the number operator have an intuitive meaning, and 
therefore we will seek a representation of the coherent states as a superposition of Fock 
states. As the coherent state is the quantum analogue of a harmonically oscillating field, 
we shall determine the corresponding electric field operator E as well. 

Formally we find the expansion of the coherent state in the Fock base with the usual 
trick of inserting the unit operator ()63|) : 

\a) = y^\n){n\a) (66) 

\n) 

The scalar product {n\a) therefore represents the expansion coefficients. Their evalua- 
tion is performed with (|Mjl. and after we multiply from left with {n\ we obtain 

(ra|a|a) = a{n\a) (67) 

Inserting the Hermite adjoint of (|59|) gives 

Vn + l(n + 1|q;) = a{n\a) . (68) 

The expansion coefficients are then obtained by the analogous application of the recur- 
sion (jnU as 

a" 
{n\a) = —= {0\a) . (69) 



"n\ 
After inserting this into ()66j) . we obtain the representation 



a) = {0\a)J2^\n) , (70) 
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which now is to be norinahzed. Normahzation dehvers the value of (0|a): 

^) = l(0|«)P$^^ = |(0|«)peH^ ,thus (71) 



(a\a} 



l|^|2 



\a) = e-^l"l . (72) 

As a result we obtain the coherent states in the Fock base expansion as 

\a) = e-^H^j:^|n) (73) 

— y/nl 



a = e 2' ' 



At this point we are able to calculate the expectation value for finding the considered 
mode in a state |a) with n photons: 



p(a,n) = Kn|a)p = ^e-H^ , (75) 

where |ap denotes the average photon number n of the mode, 

n = |ap = (a|a1'a|a) . (76) 

The expression (f73|) represents a Poisson distribution with mean value n = |ap and 
variance (An)^ = lap = n. 

Equations (fTH jI -ffTS j) show that the the coherent state |a) associated with the eigen- 
value a = coincides with Fock state |0), i.e. a state with vanishing photon number 
expectation, thus called vacuum state. 

We derived the coherent states as eigenvectors of a nonhermitian operator. According 
to the postulates of quantum mechanics we can therefore not expect that they span a 
base of the Hilbert space. In fact, their scalar product does not define an orthogonal 
relation 

(^a\p) = y S^J-!L^-kH\-m' = e["*/3-i(l"l^+l/3l^)l . (77) 

n 

(Note: Although the coherent states are not orthogonal this does not automatically 
imply that it is impossible to expand an arbitrary state in a series of coherent states 
[H].) The absolute value of the scalar product is 

|(a|/?)P = e-l"-^l' . (78) 

Loosly speaking, the coherent states \a) and |/3) become increasingly more orthogonal 
the more they are apart. This is due to the overlap of the state vectors, which is given 
by the zero point fluctuations. In fact, the coherent states can also be represented as a 



19 



displaced vacuum state. Intuitively this corresponds to a classical field amplitude with 
added zero point fluctuations. In operator notation: 

|«) = D(a)|0) , (79) 

where the displacement operator D{a) is represented as jH]: 

D(a) = exp (aa"^ - a*a) . (80) 

2.6.3 The operator of the electric field E 

According to the correspondence principle we construct the electric field operator E 
corresponding to (jMj) by replacing the classical canonical variables q and p with the 
operators (|^ and (j^ . For the mode we consider this leads to 

E{t) = ^^ [p{0) cos{k ■ f - ujt) - ujq{0) sm{k ■ f - out)] , (81) 

and when we insert (I44|) and ()45|) 



EW = h/lSb^ Kty'-' - ^\t)e-''l , (82) 



where in the literature the factor 



is often called electric field per photon. Expression ()81|) particularly illustrates the fol- 
lowing points: 

• The field consists of two components oscillating in quadrature. 

• The operators q(t = 0) and p(t = 0) associated with the quadratures are Hermitian 
and constant in time. 

• However, q(t = 0) and p(t = 0) are not compatible. That means the quadratures 
are not simultaneously measurable quantities, a result that from a classical point 
of view can not be anticipated. 

• The quantum properties of the field, including its fluctuations, are determined at 
one time, for example t = 0, and do not change at later times (provided the field 
is not further manipulated). 

• As the two field quadrature operators are not compatible, we should not be sur- 
prised that it is not possible to construct a Hermitian operator for the phase of 
the field. 
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2.6.4 The uncertainties of the coherent state 

Remember that the coherent states are defined as eigenstates of the annihilation operator 
and not of the quadrature operators q and p. Therefore, the expectation value of latter 
observables can not be exactly defined, but will exhibit an amount of uncertainty. In 
addition, we have seen that the quadrature operators are incompatible, thus, they obey 
the uncertainty relation (021) • Let us now determine the uncertainty product AgAp for 
a coherent state. (Note that {AAf = ((A - (A)!)^) = (A^) - (A)^.) 

We calculate the expectation value of the observables q and p for a coherent state |a) 
by inserting (|^ and (|^. Referring to (jM|) . we obtain 

(q)coh.st. = («|q(0)|a) = ^/h/2^ [a + a*] (84) 

(p)coh.st. = («|p(0)|«) = i^/hIJJ2 [-a + a*] , (85) 



operator q^ 



After application of the commutation relations ()38|) - Pm) and ()5H) we find for the 

q2 = A [a2 + at' + aa^ + a+a] (86) 

a' + at^ + 2ata + 1] , and (87) 



A f^2 j_ ^t2 

2uj 



(q')coh.st. = (a|q» = ^ [«^ + a*^ + 2a*a + 1] 
The width of Ag for a coherent state |a) is then given by 



(A?)coh.st. = (q') - (q)' = I; • (89) 



Analogously we obtain for the conjugated observable p 



and for their product 



(Aj>)L.st. = (P^) - (P)^ = ^ , (90) 

(Ag)coh.st.(A|')coh.st. = - • (91) 

This results shows that this uncertainty product assumes the minimum value allowed by 
Heisenberg's famous relation; cf. (jl2|). Thus, the coherent state is a minimum uncertainty 
state, and can be considered to be as close to the equivalent classical state as quantum 
mechanics permits. 

To conclude this section we use fl84|) and (jH3|) to visualize the nature of the coherent 
state. For that purpose we write the expectation values of the quadratures in their time 
dependent form. In (J8T| we can see that already the operators exhibit a purely scalar 
time dependence (c-number). Thus 

(q(i))coh.st. = (tt|q(i)|a) = ^/2h/uj\a\cosuJt (92) 

(p(t))coh.st. = («|p(t)|«) = -V2}tw |a|sina;t . (93) 
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Figure 4: Quadrature representation of a coherent state |a) according to equations (jU^ - 
(jHEl)- Note that the quadrature operators are not compatible, and therefore, 
can not be simultaneously determined. Thus, the uncertainty circle in this 
picture should not be interpreted as the contour of a classical joint probability 
distribution of observables X and F, but as Wigner distribution function. 

The expectation values for the quadrature observables q and p of the coherent state thus 
evolve in the way of a classical, harmonic field. With operators normalized to the photon 
energy we get 



and for the expectation values respectively 

('^/coh.st. 
( ' /coh.st. 

(AX)coh.St. 

(AF)coh.st. 



1 



2huj 



\q\ COSiVt 

\a\ sinurt 
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(94) 



(95) 
(96) 
(97) 
(98) 



We picture these equations in Fig. HI (the quadratures are represented in units of the 
photon energy hu). The uncertainty is equal to the vacuum fluctuation energy ^hw of 
the considered field mode, and the amplitude |a| is characterized by the average photon 
number. In summary. Fig. |3] illustrates that the coherent state can be interpreted as a 
classical, harmonic field (amplitude \a\) with added vacuum fluctuations (uncertainty 
circle). 

We said that the coherent state represents the quantum mechanical analogue of a 
classical, harmonically oscillating field. Most people connect this with images of same 
sort of antenna and radio frequency fields. In fact, we are used to modern electronic 
devices which process radio frequency cycle times of up to some tens of ps (10~^°s). 
For optical frequencies, however, electrons would have to respond on the 10~^^ s-scale. 
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Circuits with such fast responses are not available yet. The standard way to receive 
optical signals is therefore not with an antenna, but with a photo detector, such as a 
photodiode or a photomultiplier. In a photo detector an incident photon liberates an 
electron. The generated photoelectrons may then be processed with available electronic 
devices. For the ideal detector (photon to electron conversion efficiency, or quantum 
efficiency, r^ = 1) we can thus regard the electron current at the detector output as an 
exact image of the photon current at the input. 

In practice we are thus faced with the following question: how does the field state 
translates into practically observable features of the photo current? This question can 
only be answered in statistical terms. If we assume that the optical field consists of one 
mode in a coherent state, then ()75|1 contains the answer. The probability to find the 
detector output releasing n electrons is therefore given by the same Poisson distribution 

p{n) = —e-^' , (99) 

n! 

The Poisson distribution characterizes a point process (here the process of outputting 
photons, or electrons respectively), that means that the quanta are emitted indepen- 
dently and therefore only exhibit delta function like correlations. Loosely speaking, the 
Poisson distribution is the distribution with the "most random" properties - most ran- 
dom also in the sense that the underlying process can be realized in the greatest number 
of different ways. It is interesting to realize that the state that reveals the most random 
distribution if we detect its quanta, reveals a perfectly determined oscillatory evolution 
in case we choose to detect its electric field. Anyway, from the properties of the photo 
electron distribution (jUUI) and the underlying Poisson (or point) process we can deduce 
the statistical measures we are interested in, like photo electron current average over an 
interval r (which is nr), its standard deviation (which is ^/^), or its spectrum (which 
for a point process is the Fourier transform of delta function). 

Let us now consider a field, which was emitted by a single mode laser, say 1 mW 
at 500 nm wavelength. Because the laser is well stabilized this field corresponds to a 
coherent state. According to above discussion the number of emitted photons per second 
n/r corresponding to 1 mW amounts to 2.5 ■ 10^^ photons per second, and the standard 
deviation (RMS shot noise) of the current is 5 ■ 10'^ photons per second. (After photo 
detection the photons appear as electrons.) In this quantum detection picture, what 
makes us saying that the coherent state represents the quantum mechanical analogue 
of a classical, harmonically oscillating field? The quantum mechanical aspect of the 
coherent state is its noise, here 5 ■ 10"^ s^^. The classical aspect of the coherent state 
becomes visible if we imagine adding or removing one photon to the coherent state. In 
principle we disturb the quantum state by such a manipulation. But already for coherent 
states with few photons n, this perturbation is smaller than the RMS shot noise level. 
In above example the perturbation is in the order of 10~^^ of the average current, which 
fiuctuates naturally by a factor of 10^^. Clearly such a perturbation is very hard, if ever, 
to detect. The classical property of the coherent state therefore is its insensitivity for 
disturbances on the level of a few quanta. This quasi classical behavior of the coherent 
state forms the basis on which the semidassical approach for describing the interaction 
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of atoms or molecules (quantum objects) with the optical field (classical object) is built. 
For that purpose we keep only the terms that are significant for a field with a large 
number of photons (limit h ^ 0. That is, we cut off the factor ^ in the Hamiltonian 
(jHH), resulting in the semiclassical Hamiltonian 

H = 5^;^k,-at^(t)ak,(t) . (100) 

Another interesting property is the following: A coherent state can be generated by 
a point process. If we attenuate the field with another point process, the result is still 
a coherent state, although one with a reduced n. Most practical attenuators, like grey 
glass, or semi transparent mirrors, are such attenuators. This is certainly a property 
we are used to expect for a classical, harmonically oscillating field, but which is not 
automatically true for an arbitrary realization of a quantum state. 

Where do we stand now? When we started this discussion we assumed a source 
free field. As we have stated at the beginning, we have to imagine this as the field 
that subsists after a source located far away has stopped to emit. This assumption 
allowed us to learn about some essential features of the basic quantum properties of the 
electromagnetic field. For example, (j8ip illustrates that if the source were harmonic, 
the state of the radiated field consists of a coherent state, and the spectrum an observer 
records shows a single sharp line at frequency u. From a practical point of view such 
fields are not very useful. Who ever wants to wait until some charges in a distant galaxy 
stopped to wiggle so that we can receive their coherent field in our laboratory to start 
the important spectroscopic investigation we plan to? No, we need fields we can turn 
on and off in the way we need them. That means, we need the sources in our hands, 
in our laboratory. The sources we are particularly interested in are molecular sieve 
microcrystals that contain fluorescent dye molecules, and under favorable conditions 
these molecular-sieve-dye compounds can emit laser radiation. 

So let us proceed the discussion with the simplest source of an optical fleld, a single, 
flourescent molecule. This is not an unrealistic situation. It is the very simplicity of 
this arrangement that makes a single molecule an important device for exploring the 
spectroscopic properties of various chemical and biological environments 1101 • 

3 Fluorescence in free space 

With fluorescence we designate the spontaneous emission of photons by atoms or mole- 
cules. The photon is emitted into a mode of the electromagnetic fleld, so augmenting its 
photon number. Simultaneously the atom or molecule transits from an electronic state 
of energy E2 to one of lower energy Ei. The frequency of the emitted photon then is 
u = (£'2 — Ei)/h. The presence of Planck's constant h in the description of this process 
indicates that it is a quantum mechanical process. However, real atoms, even more real 
molecules, are complicated systems, which even in their most simple realization, namely 
hydrogen, exhibit a non-trivial structure of states. We therefore are forced to simplify 
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the reality. For the purpose of this discussion we will restrict ourselves to the simplest 
model of an atom or molecule, the two-level atom or molecule JT]. In the following we 
use the term "system" to spare using bulky "atom or molecule" . 



I up) 
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Figure 5: Energy levels of a quantum system with two energy states llo) and I up) sep- 
arated by the energy Itluq. The energy of the states |lo) and |up) is evaluated 
relative to the reference energy Eq. 



3.1 Two-level system and its variables 

As mentioned above we consider a quantum system with two energy states |lo) and 
I up) separated by the energy Itluq as shown in Fig. |31 According to the postulates of 
quantum mechanics the states llo) and I up) are then eigenstates of the (non-interacting) 
Hamiltonian H^ with eigenvalues Eq ± ^Jtwq: 

HA|lo) = (Eo + i^o)|lo) ^^^^^ 

Ha|up) = (Eo- |^o)|up) . 

The states |lo) and |up) form an orthonormal and complete set: 

{\\\')=6xx', A, y = l(lo), 2(up) (102) 

2 

5]|A)(A| = 1 . (103) 

A=l 

In the last section in which we discussed the quantized field, we introduced the non- 
Hermitian operators a and a^, which lowered, respectively raised the excitation of the field 
mode by one quantum (photon) hu. Similarly we now introduce the atomic operators b 
and b' which lower and raise the energy of the atom (molecule) by hwQ. Unlike the field, 
the energy of the two-level system is restricted, i.e. has a lower bound Ei = Eq — ^Itujo 
and an upper bound E2 = Eq + ^hwo. Therefore, the effect of b on state |lo). as well as 
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of b' on I up) must vanish: 

b|up) = |lo) b"^|up) = 

b|lo) = b^|lo) = |up) . (104) 

Repeated application has the following effects 

bb"''|up) = b"''b|up) = |up) 

bbt|lo) = |lo) btb|lo) = . (105) 

We can see that bb^ and b^b have the effect of number operators with eigenvalues 0, 1 for 
the lower and upper states respectively, whereas the repeated application of the same 
operator always vanishes: 

b' = = b^' . (106) 

Those properties can be summarized in the following anti- commutation rules: 

{b, b} = {b^ bn = , 

^ / ^ ^ (107) 

{b,bt}=l , ^ ^ 

where for the anti-commutator of A and B we used the notation {A, B} = AB + BA. 
These anti-commutator relations are characteristic for fermions, and are analogous to 
the commutator relations ()48|) - (j3m) of a single mode of the electromagnetic field (photon 
field), which is a boson field. 

3.1.1 Analogy to spin-1/2 system 

All quantum systems characterized by only two possible states are mathematically equiv- 
alent. The prototype of such a system is the spin-| particle in a magnetic field. We 
thus can borrow from this formalism fT^ IT^ . Be aware that unfortunately many tech- 
nical terms keep their original names, although here they refer to a completely different 
physical system. 

To describe physical observables of the two-level system at hand, we need Hermitian 
operators. Borrowing from the mentioned theory we define the following set of three 
traceless Pauli spin operators'^ 

Ro = I 1 

Ri = |(bt + b) 

1 + (108) 

R2 = ^(bt-b) 

R3 = i(btb-bb^) . 

In the original case of a spin-| particle in a magnetic field these four operators describe 
the dynamics of the spin-system. In the case of the two-level system discussed here, 

^ Often the following spin operators are defined that are not traceless: (7+ = b , cr =b, cr^ = R3. 
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however, they are not related to any spin. But because they form a hnearly indepen- 
dent, complete set of Hermitian observables, they can fully cover the dynamics in the 
two-dimensional Hilbert space of two-level atoms or molecules. Therefore, any system 
operator can be represented as a series of Pauli spin operators 



= Y,9aRa , (109) 



0=0 



where the coefficients Qa are determined by 0. The operators ()108|) satisfy the following 
commutation and anti-commutation rules^ 



{Ri,Rm} = l^im , (/,m,n = 1,2,3) 

as well as the relations 

R'=i , (« = 0,1,2,3) 



;iio) 



4 
3 

a 
a=0 



E<^'^ 



:iii) 



In short the following relations will be useful 

[b, Ri] = -R3 

[b,R2]=iR3 (112) 

[b,R3] = b . 

We now construct the representation of the operators in terms of the two states | lo) 
and I up) of the two- level system. The procedure consists in multiplying the operator 
with the unit operator ()103p from the left side and from the right side, and observing 
relations ()104|) and ()105|) . At the end we obtain 



;ii3) 



and similarly 

RqIud\ = l|un> 

(114) 

Inspecting ()114|1 we see that the states |lo) and |up) are eigenstates of the Hermitian 
operator R3 that can thus be regarded to measure the amount of inversion in the 2-level 

system. 

'^Qmn is the fully antisymmetric Kronecker symbol whose only nonvanishing values are £123 = £231 = 
£312 = -ei32 = -£321 = -£213 = 1 [m p. 209]. 
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R3 up) 


= iup) 
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= -1 lo) 



3.1.2 System energy and dipole moments 

The energy of the two-level system is represented by the Hamiltonian H^ ()101|) which 
by above procedure can be written as 

Ha = Eo + i/z^o (|up>(up| - |lo)(lo|) , (115) 

and after using pi3|) and pOSj) we obtain 

Ha = Eo + hwoRs . (116) 

If, like in our figure (cf. Fig.E)) the lower state is the system ground state, then Eq = ^hjQ, 
and Ha = fiuJo{R3 + |) = ^ob^b. 

In the following we will elucidate the physical significance of the other Hermitian 
variables Ri and R2. Those operators are closely related to the dipole moment f-i. For 
atoms, ions or molecules the dipole moment can be defined as 

M = $^eri, (117) 

i 

where r^ is the position operator for the i-th charge e in the system. To express the 
dipole operator in terms of b, b -operators we apply the unit operator trick 

/x = (|up)(up| + |lo)(lo|)^i(|up)(up| + |lo)(lo|) 
= /l22b"^b + ^iibb"^ + ;Ui2b + /l2ib"^ . (118) 

The coefficients flij stand for the matrix elements (i|Ai|j) {i,j = 1 (lo), 2 (up)). We 
therefore note that /In and /i22 represent expectation values for the dipole moment 
in the lower and upper system state, respectively. However, we know that the dipole 
moment has an odd parity, and therefore those coefficients must vanish. Because the 
dipole moment is a Hermitian operator we thus obtain 

/x = /iiab + il*2b^ . (119) 

If the system transition from the upper state lup) to the lower state llo) is characterized 
by Am = in the real system, then /ii2 is a real valued vector. On the other hand, for 
a Am = ±l-transition (e.g. induced by polarized light), /2i2 is a complex valued vector. 

Example: Hydrogen atom; let us assume the following correspondence for 
a Am = transition: 



two-level system 


hydrogen atom 


lo) — > 


n = 1, / = 0, m = 0) (s state) 


up) — ^ 


n = 2, Z = 1, m = 0) (p state) 
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If we take the 2;-axis as the quantization axis and evaluate the specified 
hydrogen wave functions, we find |llj : 

1 28\/2 
flu = (lo|/x|up) = ^^g eaoe, , (120) 

where oq is the Bohr radius and ex,ey,ez are unit vectors in direction of 
the axes. We can see here that flu has the properties of a 3-diniensional 
Euchdean vector. With ()119p and p08|) we obtain 

/x = /ii2(b + bt) = 2/ii2Ri. (121) 

However, if the transition involves Am = 1, as in the correspondence 

two-level system hydrogen atom 

|lo) — >• |n = 1, / = 0,m = 0) (s state) 

|up) — > \n = 2,1 = l,m = 1) (p state) , 

then /ii2 is complex: 

1 28 
flu = (lo|/x|up) = — eao(e; + iCy) . (122) 

In the case of a complex flu we can always rewrite ()119p in the following 
form 

M = /^i2(Ri - iR2) + flui^i + iR2) 

= 2Re(/ii2)Ri + 2Im(/Ii2)R2 • (123) 

In this example we have illustrated how the operators Ri, R2 are related to 
the dipole moment operator n. 

As next we will derive an expression for the rate of change of the dipole moment 
operator /x. (In above example of the hydrogen atom the rate of change of n would 
correspond to the electron velocity v multiplied by e.) In general, the rate of change 
of an observable is determined by Heisenberg's equation of motion, which in this case 
reads as 

^4lM.H.l. (124) 

and which with ()llfi|l . ()119|1 . and the commutation relations (jllOj) or ()112j) results in 

^ = T^[(/2i2b + /2*2bt) , [hwoRs + Eo)] 

= -iuoiflub - flu^^) . (125) 
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The operators b and b' are transformed into the interaction picture according to the 
usual rule 



b(t) = exp 



^H^(t-to) 



b(to) exp 



-^^Ait-to) 



t>to. 



(126) 



Inserting ()116p . and observing the operator expansion theorem, as well as commutation 
rules (I112J) . we find 



b(t) = b(to) e-'"o(*-*°) 
bt(t) = b"f(to)e'^°(*-*°) , 

and for the dipole operator in the interaction picture 

/x(t) = 2Re(/Ii2)Ri(t) + 2Im(/ii2)R2(t) 

= /ii2b(to) e-*"°(*-*«) + iI*,2b\to) e''^°(*-*») 



(127) 



(128) 



3.1.3 Bloch-representation of the state 

We have seen [c.f. fll02p and ()103p ] that the states |lo) and |up) form an orthogonal and 
complete set within the two level model. Hence, any pure state of the two-level system 
can be represented as the linear combination 



Ci|lo) + C2|up) 



with 



I 1 2 I 1 2 



1 . 



(129) 



(130) 



When the system state is not pure and has to be specified in statistical terms, it is 
represented by the atomic (or molecular) density operator p^^^: 



P^^^ =Pll|lo)(lo| +P22|up)(2| +pi2|lo)(up| +p2l|up)(lo| , 

where the coefficients pij stand for the ensemble average 



(131) 



(132) 



which represents a two-dimensional, Hermitian, covariance matrix. Bloch J^ presented 
a simple, intuitive, geometrical interpretation of the state in terms of a real three- 
dimensional vector r* with components ri,r2,r3. In the Schrodinger picture they are 
given by 



ri = 2Re(pi2) 
r2 = 2Im(pi2) 

^3 = P22 - Pii , 
where often a forth component is added: 

^0 = P22 + pii = 1 . 



(133) 
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The correspondence between the Bloch vector and the density matrix is a consequence of 
the properties of the respective symmetry groups, namely the correspondence between 
the real orthogonal group 03 and the special unitary group SU2. Note, that also the 
Stokes vector that is often used to represent the polarization state of light (see e.g. P^ ) 
has 03 symmetry. 





Z 


fully ^^ 
excited +1 \ 

^''partly \ 
y^ excited | 






^^ 


y 


Xy^ 


^■■■■..... 


fully excited —1 

_.-■■' 





Figure 6: Bloch-vector representation of the state of a two-level quantum system. Pure 
states are characterized by vectors of unit length. "Fully excited" refers to 
system occupying pure state |up), whereas "fully unexcited" to |lo). 

In the three-dimensional Bloch-vector representation sketched in Fig.lHlpure states are 
characterized by their unit length. For example the pure state llo) corresponds to the 
down pointing vector r = (0, 0, —1), whereas |up) to the up pointing vector r = (0, 0, 1). 
Pure intermediate states point in various directions, in particular, states with an equal 
mixture of upper and lower states (p22 = Pii) lie in the horizontal x,?/-plane. 

As mentioned, a pure state is represented by Bloch vectors of unit length, while for 
mixed or impure states the length is less than unity: 



2,2,2 

r^ + r2 + rg 



4|pi2p + (P22 -/On 
1 -4(p22Pii - |pi2| 



From the Schwarz inequality follows 

P22P11 - |P12|' = (|C2H(|C1H - |(CiC;)|2 > 



(134) 



(135) 



where equality is realized when the ensemble contracts to a single realization, thus a 
pure state. Therefore 

r^ + r^ + r^ < 1 , (136) 

with equality only for a pure state. On the other side the Bloch vector vanishes when 
the ensemble consists of an equally weighted mixture of upper and lower states with 
random phases. 
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With ()1()9|1 we mentioned that the spin operators Rq, form a complete set, so allowing 
to represent any observables in the two-dimensional Hilbert space of the two- level system. 
Thus in this base the density operator p^^"^ has the expansion 

3 
P^"'^=5Z^aR«, (137) 

where the coefficients Qa are time dependent, whereas the operators Rq are time-inde- 
pendent. Evaluating {i\p^^^\j), {i,j = lo, up) with ()133|) we can identify the coefficients 
as Qa = ra, and therefore 

3 

p(^)=^r,R, . (138) 

This result shows that the Bloch-vector components can be interpreted as weights in the 
Schrodinger-picture expansion of the density operator in the spin operator base. 

The symmetry of the Bloch-vector representation (cf. Fig. ^ suggests a representation 
of r in polar coordinates {r,6,(j)) rather than Cartesian coordinates (ri,r2,r3). From 
(fnij) and ET^ we deduce 



(139) 



= i(l-^3) = 


|(1 — rcos 


= 1(1+^3) = 


|(1 + rcos 


= i(ri + ir2): 


= ir sin6'e'' 



resulting in 



^ P22 - Pll , 

cosy = , and 



Im(pi2) 



T 

(140) 



3.1.4 Spin operator expectation values 



In fll38|) we have illustrated the relation between the spin operators R^ and the Bloch- 
vectors components r^- In the following we calculate the expectation values of the spin 
operators for the system in an arbitrary quantum state which is characterized by the 
density operator p^-^': 

(R„)=Tr[p(^)Rj , (a = 0,...,3). (141) 

Inserting p38|) for the density operator in the Schrodinger picture we obtain 

3 

(R«) = ^Tr(r^R^R„) . (142) 

/3=0 
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The product of two different operators R gives either Ri, R2, or R3, which all have 
vanishing trace. The sum therefore only gets a contribution from term a = (3. Since the 
trace of a two-dimensional unit vector is 2, we find with (jlllj) 

(R„) = Tr(ir,l) = ir, , a = 0, . . . , 3 . (143) 

The components of the Bloch vector can thus also be interpreted in terms of the expec- 
tation values of the spin operators. 

We can now calculate the expectation of the energy (|116|) as 

(H^) = Eo + \nuj^r^ . (144) 

In the same way the expectation value of the dipole moment (jl23|) is given by 

(^) = Re(/ii2)ri + Im(/ii2)r2 . (145) 

Prom this equations we can see that the z-component of the Bloch vector is associated 
with the energy of the two-level system, whereas the x- and y-components are related 
to its dipole moment. For a two-level system without a permanent dipole moment the 
expectation value of the operator /x must vanish in the lower, as well as in the upper 
state. In those states the operator n is not well defined and its value fluctuates. From 
(frrn|l or |T^ we can see that 

/x' = |/2i2|'l . (146) 

As a consequence, the variance of the dipole moment in the lower state llo) and in the 
upper state |up) is given by 

(1o|(Am)'|1o> = I/I12P = (up|(Am)>P> , (147) 

and is maximum in those states. Therefore, even if the expectation value of the dipole 
moment of a two- level system vanishes in the ground and upper state, the system can 
interact with an electromagnetic field through the nonvanishing fluctuations of /x. In 
the following we will now discuss these interactions. 

3.2 Interaction of a two-level system with a classical 
electromagnetic field 

Let us assume that the classical electromagnetic field is described by its electric field E, 
and that the two-level system is located at a fixed point in space and its dipole moment 
is characterized by the operator ^. Also, the typical wavelength of the field is much 
larger than the length of the dipoles. The interaction energy H/ is then given by the 
usual expression for the potential energy of a dipole in an electric field. With p23|) we 
can write 

Hi = -At(t) ■ E = -2[Re(/ii2)Ri + lm{fIu)R2] ■ E (148) 
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The interaction of a charge with an electromagnetic field of moderate power^ can also 
be described by the interaction energy 

Hj = -p{t) ■ A{fo,t) . (149) 

m 

For the canonical momentum operator p we can set mfi/e. and with p25|) we obtain 

Hi = -ii{t) ■ A{t) = icuo[/ii2b(t) - ni2^\t)] ■ A{t) . (150) 

We have given two interaction Hamiltonians, ()148|1 and p5()j) . Depending on the problem 
at hand one can choose the more appropriate one, although in most cases they lead to 
similar results. The total energy of the system in the field consists of the sum 

H = H^ + H, , (151) 

where H^ is diagonal in the states |lo) and |up), whereas H/ is off-diagonal. 

In the Schrodinger picture the time dependence of the density operator is given by 
the equation 

' [(H^ + H,), p(^)(t)]. (152) 



dt ih 

3.2.1 Bloch equations 



As next we consider ()152|) in matrix element notation. After inserting the unity operator 
(I102J) between H and p'^^\t) in the commutator relation on the rhs of p52|) . we obtain 
the following Schrodinger picture equations 

Pll = TT [(lo|H/(t)|up)p21 - C.c] 
P22 = -TT [(lo|H/(t)|up)p21 - C.c] 

^^^ (153) 

Pl2 = TT [-^0Pl2 + (l|H/(t)|up)(p22 - Pll)] 

P21 = TT [hjQP21 + (2|H/(t)|lo)(pii - P22)] • 

Using p33|) we can express ()153p in terms of the Bloch vector and the result is often 
referred to as Bloch equations 13j: 

2 
fi = -Im((lo H/ up))r3-u;o'^2 

r2 = -\ Re((l|H,|up)) rg + u^n (154) 

2 I 2 I 

rs = -- Im((l|H7 up)) ri + - Im((l|H7 up)) rg 



Afield power below the regime where multiphoton interactions dominate, so that terms in A and 



higher can be neglected 
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We can see that r^ is constant, when there is no external field. Multiplying the first 
equation with ri the second with r2 and the third with r^, we find 

^^irl + rl + rl) = Q. (155) 

Thus, in the presence of a classical field the length of the Bloch vector remains constant. 
That means a pure state stays pure, and a mixed state remains mixed. The structure 
of ()154|) led Feynman et al. ^H] to propose an interesting geometric interpretation. For 

— * 

that they introduce a vector Q with following components 



(156) 



Qi 


= |Re((lo Hi 
n ^ 


up» 




Q2 


= ^Im((loH, 


np» 




Q3 


= ^0 • 






With this vector the Bloch equations p54() are equivalent 


to 




d. ^ ^ 
-—r = Q X r 
at 







(157) 

In analogy to the mechanics of rigid bodies this equation shows that the Bloch vector 
r precesses around the vector Q with a frequency that is given by the magnitude of Q. 
When Q itself varies with time, the precession motion may become complicated. We 
may evaluate the matrix element ((lo|H/|up)) using ()148p or ()150p . For example, with 
(fTiH|) we obtain 

(lo|H,|up) = -/ii2^(t) . (158) 

Let us assume an electric field with frequency ui that is nearly resonant with the two- 
level system, so that \ui — uqI <t^ uq. Then we can write 

E{t) = e£{t)e'''^' + c.c. 

= 2\£{t)\ {Re(e) cos[uJit - 0(t)] + \m{e) sin[ujit - 0(t)]} , (159) 

where 8{t) = \£{t)\e^^^^'^ is a slowly varying, complex amplitude, and e*is a unit polar- 
ization vector. With this we obtain 

(loJH/lup) = -|^(t)| /2i2 ■ [ee"'["^*-*(*)l + c.c] . (160) 

This matrix element can now be substituted in the Bloch equations (jl54j) . 

3.2.2 The rotating wave approximation 

As an example to illustrate the rotating wave approximation, let us compare now tran- 
sitions with Am = 1 and Am = 0. In the case Am = 1 the vector fli2 has the form 

jli2 = -1=-^ I/I12I , {ex,ey are unit vectors) . (161) 

v2 
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In addition, assume an external field that is circularly polarized and propagating in the 
2;-direction, thus 



V2 
Then we obtain for the scalar products 



(162) 



/2i2-£ = 0, and flu ■ £" = If^ul ■ (163) 

The Bloch equations ()154|) then become 

r-i = -i7(t)r3 sin[a;it - (j){t)] - ujQr2 

7-2 = fi(t)r3 cos[cJit - (f){t)] + uori (164) 

fs = —Q(t)ri sin[u;it — 0(t)] — Q{t)r2 cos[ct;it — 0(t)] , 

where we have introduced the parameter 

nit) = lfi,2-^\£{t)\, (165) 

n 

which is known as the vacuum Rabi frequency. The Rabi frequency in this context is a 
measure for the coupling strength of the two-level system with the external field. For 
the components of the Q-vector ()156|) we obtain 

Qi = -n{t) cos[uJit - <p{t)] 

Q2 = -^(t) sm[uJit - 0(t)] (166) 

Qs = ^0 • 

On the other hand, for a Am = transition /I12 is a real vector, and if the incident 
light is linearly polarized, then £*is also real. In this case the Bloch equations (jl54j) read 
as 

ri = -uJor2 

7-2 = 2n{t)r3 cos[uJit - 0(t)] + uon (167) 

rs = —2Q(t)r2 cos[u;it — 0(t)] , 

and following ()157|) we obtain for the corresponding Q-vector 

Qi = -2n{t) cos[c<;it - 0(t)] 

Q2 = (168) 

Q3 = ^0 • 

Looking at the two cases of the field interaction, the first with a Am = 1-system, the 
second with a Am = 0-system, ()164j) and p66|) seem to be rather different from p67|) 
and p68|) . However, Allen and Eberley 11, sect. 2.4] have shown that they only differ by 
some anti-resonant terms. Comparing the expressions for Q we can see that ()168p can be 
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decomposed into a sum of two vectors, of which the first vector is p(j(j|l . and the second, 
auxihary vector consists of components (— f2(t) cos[uJit — 0(t)], +Qit) sin[ti;it — 4>(t)], 0). 
This auxihary vector rotates around the z-axis with frequency ui — 0(t), whereas the 
Bloch vector r rotates in the opposite direction with frequency uq- Their relative rotation 
frequency thus amounts to ujQ+uJi—(j){t). Thus, when integrating the equations of motion 
over a characteristic time interval, say At > I/cuq, the contributions of the fast precessing 
auxiliary vector Q are small. To a good approximation we can therefore neglect the 
auxiliary vector Q. If we drop the auxiliary Q, then the Am = interaction with linearly 
polarized light is described by the same set of equations as the Am = 1 interaction with 
circularly polarized light. This is known as the rotating wave approximation. 



3.2.3 Bloch equations in a rotating frame 



Let us go back to (|164j) . This set of equations describes a rotation around the 2;-axis 
at an optical frequency. We will now introduce a rotating reference frame, in which the 
motion of the Bloch vector r is slower. At first sight, Eq. ()164|) seems to justify the 
atomic (or system) frequency tuo as a suitable choice. But after we consider that the 
atomic frequency in a spectroscopic sample varies from atom to atom, we realize that the 
frequency of the applied field ui is better suited if we want to refer the time evolution 
of all the differing atoms to the same reference frame. The transformation from the 
stationary frame r = (ri.r2,r3) to the rotating frame f' = {r[.r2,r'^) is given by 

r' = ef, (169) 

where is the 3x3 orthonormal rotation matrix 



©= I —sinujit coscuit (170) 




Let us now transform Bloch equations ()164|) to the rotating (primed) frame. For that 
we rewrite (11641) in matrix form as 



f=Cf, (171) 

where C is the 3x3 coefficient matrix. Inserting p69|) we obtain 

dr ' • ^ 1^ 

—— = r + r 
at 

= 00^0f+0C0^0r (172) 

= (00~i + 0C0^)f' 

After we insert the elements of C from p64|) and p70|) we get the Bloch equations in 
the rotating frame 

r'l = (cji — uJo)r2 + fi sin r'^ 

r'2 = {ujq — uji)r[ + Q cos (j) Tg (173) 

r^ = —Q sin (p r'^ — Q cos rg , 
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which with Q = {—Q cos (j) , Q sin (f) , ujq — ooi) can be expressed as 

^ = Q'xf' (174) 

We can see that the Bloch vector r ' precesses around Q ' with the frequency depending 
on |f2^ + (cui — C(Jo)^|^'^^ and the orientation of Q'. If the initial direction of the Bloch 
vector f' is approximately parallel to Q', then r' precesses around Q ' on a cone with 
small angle, and this cone tends to follow slow variations of the direction of Q '. This is 
called adiahatic following, and can be used to prepare the quantum states of atoms. 

3.2.4 The Rabi solution 

In the following we discuss the interaction with a sinusoidal exciting field (well stabilized 
laser), so that in ()159|) the complex amplitude £ is constant and the phase can be made 
zero by a proper choice of the time origin. Historically the solution of the interaction 
of a two-level system with such a field was given by Rabi jT2] when he studied a spin 
I system in a magnetic field. For our discussion here, let us assume that at t = the 
system starts in the lower state |lo), then r^{Q) = — 1 and ri(0) = r2(0) = 0, and the 
Bloch equations in the rotating frame are given by 

These equations characterize an intricate motion. In ()133|1 we have seen that the rs 
component of the Bloch vector is related to the population inversion of the two-level 
system. Eq. p75|l reveals an oscillation of population characterized by r'.^{t) around 
{ijjQ — ijJiY/[Vl? + {ijjQ — cui)^] with frequency [Vt^ + {uq — cji)^]^/^ and amplitude Vt^ /Vt^ + 
(tuo — cji)^. This phenomenon is known as Rahi oscillation or optical nutation, and can 
be observed in systems, which are well isolated from disturbances, like in a low pressure 
gas fZl^], that means in systems, in which damping time constants are longer than 
the time scale of system motion. 

3.3 Interaction of a two-level system with a quantum field 

According to the postulates of quantum mechanics, a two-level system initially in its 
excited state I up) will stay there for all times in the absence of an interaction. Experience 
shows, however, that normally the system will return to its ground state llo) after a 
certain time. In the last sections we discussed how the interaction with a classical field 
can change the occupation of states. On the other hand, spontaneous emission is a 
quantum mechanical process, and its description requires the quantization of both, the 
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atoms and the field. A result of this theory is that the rate of spontaneous emission is 
proportional to the mode density of the surrounding environment. 

To illustrate this theory we start with the interaction of a two-level system with a single 
mode of the electromagnetic field. The total Hamiltonian thus consists of a component 
describing the two-level system H^, the component of the electromagnetic field Hp, and 
the interaction Hamiltonian H/ 

Htot = H^ + H^ + H, , (176) 

where the two-level system Hamiltonian H^ is given by ()115|) . and where for the purpose 
of this short overview we insert the field Hamiltonian Hp in its semiclassical form ()100|1 
restricted to one mode kg. The interaction is described by a Hamiltonian analogous to 
p48j) . in which the field is expressed as Hilbert space operator 

H, = -/x • E . (177) 

Thus, the quantum properties of the field enter at two places, first through Hp pOOp . 
and second through H/. In Hp 

Hp = huj^^ ■ al^{t)a^^{t) (178) 

the field is represented as the product of annihilation and creation operators. Even 
though in this semiclassical Hamiltonian the zero point energy term ^/icuk^ of the fully 
quantum field Hamiltonian ()54|) is missing, a part of the quantum reality is still repre- 
sented in the noncommutativity of a|^ and ak^; cf. (I^H|) . On the other hand, in (jl77p the 
quantum aspect is introduced by the electric field, which appears as the Hilbert space 
operator (jH^ 

E = i^"'^/^ [ak, e'^^-^^ - aj,^ e^'^^'l . (179) 

In this notation E describes a travelling wave mode^, and e*is the polarization vector. 
In the rotating wave approximation the interaction Hamiltonian ()177|) can be written as 

[cf. dm] 

H/ = ^(aLJlo>(up|+h.c.) 

^ (180) 

(al b + h.c.) , 






where the coupling coefficient f2 represents the vacuum Rabi frequency (|165|) . which in 
general depends on the location r of the two-level system. The total Hamiltonian thus 
is given by 

Htot = ^0 + ^^oRs + ^ke ■ al^ak, + — (aj^^ b + h.c.) . (181) 

In (I18()j) we have represented the interaction Hamiltonian in a form which permits an 
intuitive interpretation: the two-level system can 



5 



Depending on the geometry of the problem, a standing wave mode can be more suitable 19^. The 
corresponding operator is then given by: E = £\/ j^ sinfcr (a + a^) 
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• absorb a photon from the field and make a transition from state |lo) to |up), or 

• emit a photon to the field and make a transition from |up) to |lo). 

A state |\E') of the combined quantum system consisting of the two- level system in state 
a) and the field say in a number state In) can be represented symbolically as |\Ef) = 
A, n). Within the framework of the model described by (|181|) a state |\l/) = |up, n) can 
only couple with state |\1/') = |lo, ?2 + l). Consequently, we are allowd to consider the 
interaction for each manifold of levels lupj-n.), [lo, « + l) independently. Obviously, in 
each manifold the number of excitations amounts to n + 1 and is conserved. Technically, 
this means that for each manifold the problem is reduced to solving a two-level problem. 
The full dynamics is then obtained by summing over the dynamics within the appropriate 
manifold. 

To illustrate the basic concept of spontaneous emission, let us assume that the two- 
level system is initially in its excited state |up) and interacts with only one field mode 
which is in the vacuum state |o), 

|^(t = 0)) = |up,0) . (182) 

As stated above, the two-level-field system remains in the one-quantum excitation man- 
ifold for all time (given the validity of psil) ). and its state at time t is given by a 
superposition of |up,0) and |lo, l). At resonance ujq = cj^, the state is described by 

|^(t)> = cos (^^ |up, 0) - isin (^\ |lo, 1) . (183) 

The probability to find the two- level system in its ground state |lo) after time t is 

P,o=Klo,l|^&(t))p = sin2(^) . (184) 

This result shows, that the initially unoccupied ground state |lo) spontaneously becomes 
occupied, even though the field initially was in a vacuum state. This is due to the effects 
of the quantum nature of the field, which we pointed out at the beginning of this section 
13.31 Note, however, that if the two-level system is initially in its ground state |lo) and 
the field in vacuum state |0), then the two-level system will remain there for all times. 
On the other hand, above result for P\o exhibits an oscillatory behavior at Rabi fre- 
quency Q, which is in contrast to experience where the decay is irreversible. In the simple 
model described by I^181\i the two-level system interacts with only one field mode, which 
also is undamped. Therefore the two-level system can reabsorb the photon from the 
field, so returning to its initial state. Later on we will discuss how recent experimental 
progress has allowed to observe this oscillatory regime. But for the moment, we will 
discuss how spontaneous emission becomes irreversible, when the photon is emitted into 
a multimode vacuum field. 
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3.4 The Fermi golden rule 

A realistic model of a two-level system in free space involves the interaction with a 
multimode field. Instead of ()18H) we therefore consider a Hamiltonian of the following 
form 

Htot = ^0 + ^oRs + XI ^k, ■ aj,^a^ + ^ -^ (aj^^ b + h.c.) , (185) 

kj kj 

where the last term describing the interaction is the multimode interaction Hamiltonian 
H/; cf. the single mode Hamiltonian ()18()|1 . 

Similar as above in the single mode case we assume that all the modes of the field are 
initially in their vacuum state and that the two-level system is in its excited state, so 
that we can write 

|vl/(t = 0)> = |up,{0}>, (186) 

where |{0}) labels the multimode vacuum, and for each mode k^ we have ak^|{0}) = 0. 
The state of the combined system can then formally be described by the superposition 
of states 

|*(t)) =a(t)e-^-«*|up,0) + J]6k^(t)e-^'^''«*|up,lkJ , (187) 

ke 

where the coefficients at t = are set to fulffil the initial conditions, thus a{t = 0) = 
1, by^^(t = 0) = 0. For convenience we have extracted the fast time dependence, and we 
have reset the origin of the two- level system energy so that hw\o = and hw^p = TujJq. 
The field state | Ik^ ) can be represented by 

|lk,>=aLjO>. (188) 

The equations of motion for the coefficients a{t) and h{t) are readily obtained as 

^ = -^E^k,e-K--X(t) (189) 

^M^ = _il]*e^K--«)*a(t) (190) 

If we restrict ourselves to times t close enough to t = so that the initial coefficients do 
not change significantly, we can approximate a{t) in (jl9(jp by its initial value a(t = 0) = 1 
(this approximation is often named ffist order perturbation theory). With this we can 
readily integrate ()190p and we obtain 

4 ^(o^k^-^o)^ 

The probability to find the two-level system in its excited state after time t is then given 
by 

P,p = l-5^|6k,(t)P, (192) 
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where the sum collects the contributions of every mode which the two-level system can 
couple with. Here we can see for the first time, how the structure of the mode space 
affects the spontaneous emission: As the sum only gets positive contributions it clearly 
increases when the number of participating modes increases, so reducing the probability 
Pup- Since the frequency bandwidth over which the dipole moment couples with the 
field is limited, so actually the number of modes per frequency interval is the significant 
parameter governing the size of the sum. And this parameter depends on the geometry 
of the space in which the fluorescent system can radiate into. This is the observation 
on which attempts to modify the spontaneous emission rate will hook on. We will come 
back to this in a moment. 

When the two-level system radiates into a space which is densely populated with 
modes, such as free space for example, then it couples to a whole continuum of modes 
over which we have to extend the summation. Mathematically this means replacing the 
sum in fll92|) over the modes kg by an integral. We discussed this problem in part 12.31 
where we obtained the transformation rule ()28|) . If we equate ( in (j^Hj) with ()191|) where 
we insert (jl65j) for the vacuum Rabi frequency with (j83|) as electric field, we can work 
out the integrals and we obtain 



P -1 ^ f A; ,;3 I,; |2 Sm'[lK-^o)t] 



(193) 



Now, as time t moves on, the term with sin^ [K^^^ — ^o)'t\/[^iuJk^ — cuq)]^ will be 
significantly over zero only for modes with frequency cuk^ ~ uo- This term therefore 
secures energy conservation in the interaction; in fact mathematically 

hm l'^"-' '^ = 2nt 5(^k, - ^o) . (194) 

Inserting this we find for the decay rate of the excited state in a free space vacuum 
environment (which is also known as the Einstein A- coefficient) 

r\P / jS I ,7 12 1,7 |2 

_ Q-Tup _ ^uJq^ |^12| _ _|, |3 IM12I dQSl 

dt c^ Sntoh Sntoh 

This is an example of a special case of Fermi's golden rule, which predicts that for 
times large enough that energy conservation holds, but short enough that first-order 
perturbation theory applies, the excited upper state irreversibly decays according to 
p95p . In summary, when all the approximations described above are considered the 
Fermi golden rule for a general transition can be expressed as 

2n . 

Fermi golden rule: 7 = — N^|^/|H/|i)p 5(u;j — cjj) , (196) 

ft/ 

f 
where 7 characterizes the decay rate from the initial state |z) to the final state |/). 

Let us go back to the initial first-order perturbation theory assumptions, in which we 
agreed to restrict ourselves to short times t so that for the coefficients a and b we can 
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write a(t = 0) = 1, 6k^(t = 0) = 0. Of course, when this holds, then the upper state is 
still nearly fully excited, and Pup ~ 1- Thus we still stay within the limits of our initial 
assumptions when we write 

^^-7fsPup, (197) 

which of course corresponds to an exponential decay with decay rate 7fs. 

In the next paragraph we will show that this is true even for longer times than we 
have considered here. 

3.5 The Weisskopf-Wigner theory 

An other approach to solve the dynamical equations for the coefficients a and b ()189p . 
(I19(J|) was introduced by Weisskopf and Wigner [20] • They start by formally integrating 
(|19U|) . and inserting this result in (|189p . which will produce 

^ = - E ^ /* dt'e-^K--')(*-*')a(t') . (198) 

As above and for the same reason we replace the sum over the modes by an integral, 
which results to 



\— [ do^k, <|/xi2p rdt'e-^K--o)(*-*')a(t') . (199) 

•" "-C Jkf Jo 



da{t) _ 
dt 6eo7r 

The construction of the dependent superposition of states ()187p was such that compared 
with all other time dependent functions the variation in a and b is slow, and therefore 
we assume that the variation of a{t) in (jl99|) is much slower than in the exponential as 
well, and can thus be pulled out of the time integration. At the end we may check, if 
this assumption is consistent with the result. Because of 



lim / ^t>^-K-^,--o){t-t') ^ ^^. -uj^)-V 



^k^ -CJo. 



(200) 



for larger times we have a similar situation as in ()194j) . where here now the exponential 
assures energy conservation. Neglecting the frequency shift due to the principal value 
term (which is analogous to the Lamb shift) we find 

^ = -^a(t), (201) 

dt 2 ^ ^ ' ^ ' 

or 

^ = -7fs Pup , (202) 

which is the same as we guessed in fll97j] . 
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3.6 Reservoir theory and master equation 

In the last two paragraphs we have seen how the irreversibihty of spontaneous emission 
emerged when the source was coupled with a quantized field, more exactly, its vacuum 
modes. Notwithstanding, the Hamiltonians are perfectly energy conserving, so the ap- 
parent nonreversible dynamics of the system comes somewhat surprising and its physical 
reason is rather obscure. With the reservoir theory approach we will discuss as next, 
we will be able to gain a more intuitive understanding of the physical processes which 
lead to the nonreversible decay of the excited level of a two-level system j^H p. 374ff]. 
In reservoir theory we shift the perspective from the particular two-level system coupled 
with the field to a more general situation, which is a small system coupled to a large 
system,. We characterize the small system by its Hamiltonian Hs, the large system by 
H^, and their coupling by the interaction Hamiltonian V. Thus: H = H^ + H^ + V. For 
our particular case the small system can be identified with the two-level system and the 
large system with the continuum of modes of the field. In addition we assume that the 
large system always stays in thermal equilibrium at some temperature T. This means it 
acts as a thermal reservoir. A thermal reservoir is usually described by an equilibrium 
(i.e. time independent) density operator of the form 

p, = -exp{-Hr/kBT), (203) 

where with ks we denote Boltzman's constant, and the partition function Z is given by 
the trace over the reservoir Tr^ 

Z = Tir [exp {-Hr/kBT)] . (204) 

Now, we are interested in the dynamics of the small system only. In that case we can 
find the dynamics in the evolution of the reduced density operator p^ 

P, = TTr{Psr). (205) 

where p^^ is the density operator associated with the full system, i.e. small system and 
reservoir. Thus the reduced density operator p^ is the trace over the reservoir of the 
total density operator. If we know the reduced density operator at any time t we can 
calculate the expectation value of any system operator. The equation of motion for p^^ 
is called a master equation., and this is what we will derive in the following. In order 
to focus directly to the relevant dynamic time scale of the system, we switch to the 
interaction picture, where all the free evolution is eliminated. The interaction between 
the small system and the reservoir is described by the Schrodinger picture interaction 
Hamiltonian V, for which the interaction picture representation \//(t — tg) is obtained by 
the unitary transformation 

V,(t - to) = exp [iHo(t - to)/h] V exp [-iHo(t - t^)/h] , (206) 

and where we have set Hq = H^ + H^. Similarly we can relate the full system density 
operator in the interaction picture P^^ to the Schrodinger picture density operator p^^ 
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by the unitary transformation 

p,,{t) = exp [-iHo(t - to)/h] P,, exp [iHo(t - to)/h] . (207) 

Observing the Schrodinger picture rule p = — ■|[H, p] we obtain the following time deriva- 
tive 



dt h 



exp [-iHo(t - U)/h] |[Ho, P..(t)] + ^1 exp [iHo(t - to)/h] , (208) 



in which the motion of the density operator in the Schrodinger picture is related to its 

motion in the interaction picture. From this we obtain the interaction picture equation 

of motion 

dP i 

~W = ~h [V/(^-^o),P..(t)] . (209) 

We may assume that at t = to the small system and the reservoir do not exhibit 
any correlations. We can therefore approximately solve this equation to second order in 
perturbation theory, and obtain 



Psr{t) = Psr{to) -^ f dt' [V,(t - to), Psr(to)] - 

-^ f dt' I' dt" [V,(t' - to), [V,(t" - to), P..(to)]] + 



(210) 



We trace out the reservoir and obtain the reduced density operator in the interaction 
picture p(t) = Tr,,[Psr(t)] for which we can write 

p(t) = exp [iH,(t - to)/h] p, exp [-iH,(t - to)/H] (211) 

When the time interval r = t — to is long compared to the relaxation time (memory 
time) of the reservoir Tc, but short compared to times in which the small system variables 
show significant changes (for example 7^"^ in spontaneous emission), we can define a 
coarse-grained equation of motion (time derivative) by 

Pit) ^ '^^ - P^' - "^ . (212) 

r 

Applying this to ()210j) we obtain after some algebra 
m=-^f^dr'TT{\l:(T')Psr{t)]- 

- ir r^^' r dr"Tr{V,(r')V,(r")Psr(t) - V,(r')Ps.(t)V,(r")} +adj. . 

(213) 

It can be shown that with the properties of the reservoir given here the first term on 
the rhs of ()213p vanishes, and that the second term is composed of a sum of two-time 
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correlation functions of reservoir operators only ^TT, p. 380]. In fact, with the dipole 
coupling Hamiltonian p85|) and taking into account the transform ()206|) we have 



Vj(t) = nbtF(t) + adj. , (214) 

where the operator 

k 

acts only in the reservoir Hilbert space. As is better visible from their Fourier trans- 
formed form, such operators are usually associated with noise sources and are thus called 
noise operators. The trace over the reservoir involves first order correlation functions of 
the forms (F(t')F(t")), {f\t'p{t")), {F{t')F^{t")) , {F^{t')f{t")) . For example, 

{Hty\t'%. = J2y^^ <"^''^'^^ e^^oit'-n^(^,t'~.yt") ^ (216) 

k,k' 

where (. . . )r denotes the average over the reservoir. As we have 

{akal,)r = {nk + l)Skk', (217) 

where n^ denotes the average number of thermal photons in mode k (n^ = at zero 
temperature), the correlation function ()216|1 reduces to 

{f{t')F\t"))r = Y, ^ (^-t + 1) e'("o-"^)(*'-*") . (218) 

k 

For the other correlation functions appearing in ()213|) analogous expressions can be ob- 
tained. As we assumed a reservoir in thermal equilibrium, thus with stationary statistical 
properties, the reservoir correlations depend only on the time difference At = t' — t", 
and we obtain for an example term appearing in (j213|l 

/ V /V' (F(r')F^(r"))r. = /V f^dAt V i^ (n^ + 1) e^(-o-.)At (219) 

Jo Jo Jo Jo k 

The time dependence of fl219p is governed by the first order correlations existing in the 
reservoir, for which ()218|1 is an example. Given the reservoir properties introduced in 
the discussion of ()212j) (reservoir correlation time much shorter than characteristic small 
system time constants) we can extend the upper integration limit r' in ()219|1 to infinity, 
and we obtain 

PT j*t' J*T poo IQ I 2 

/ dr'/ dr"(F(r')Ft(r"))r= / dr' / dAt J^ ^^ (n^ + 1) e*('^«-'^^)^* 
Jo Jo Jo Jo k ^ 
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where we have replaced the sum over the modes by an integral in the same way as we 
discussed in paragraph 13.41 and 13.51 Obviously this equation has a similar structure as 
()199p . After some algebra, as inserting the delta function representation ()200|) and ignor- 
ing the associated principal part frequency shifts, combining the different contributions 
of (j21Up . and explicitly focusing on the two-level system as the small system, we obtain 
the interaction picture master equation 

Pa = -^{n + 1) [btbp^(t) - bp^(t)bt]- 

(221) 
-^n[p^(t)bbt-bV^(t)b]+adj. . 



In contrast to p99p of the Weisskopf- Wigner theory, ()220j) gives us some insight into the 
processes leading to the decay of the excited state. Together with ()22H) . ()220p relates the 
decay with properties of the reservoir expressed as first order correlation functions. The 
approximation of these correlation functions with delta functions is known as the Markov 
approximation. The decay appears thus as a result of the delta function like memory 
time of the reservoir which instantaneously looses track of the interaction with the two 
level system. So the way of the excitation back to the two level system is forgotten, so 
to speak. Combining the Markov approximation with second-order perturbation theory 
to derive the master equation is usually labeled Born-Markov approximation. 

The master equation (J22H) describes the decay of a two-level system which is coupled 
to an electromagnetic field characterized as a reservoir at temperature T. In this equation 
n represents the number of thermal quanta of the reservoir at the two-level transition 
energy hujo. For a reservoir at zero temperature, the population Pup of the excited state 
of the two-level system is given by 

Pup(t) = Tr[|up>(up|p^(t)] =Tr[(R3 + ^)p^(t)] , (222) 

and exhibits the equation of motion 

dP (f) 

-^^ = Tr [ |up>(up| pAt) ] = -7fsPup(t) , (223) 

which is the Weisskopf- Wigner result (|2(J2|) . 

Let us look again at the master equation of the two-level system (j22ip . In the 
Schrodinger picture and at temperature T = this equation has the form 

PAit) = -{[^A, Pa] - Y[b^bp^(t) + p^(t)btb - 2bp^(t)bt] , (224) 

where H^ denotes the two-level system Hamiltonian pi5|) . If we formally define a non- 
Hermitian effective Hamiltonian Heg as 

Heff = H^-i^^btb, (225) 
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then introducing this in ()224|1 leads to 

PA{t) = --^[Heff/o^W - pAml^] + 7fsbp^(t)bt . (226) 

Going back to (jl87|l we define the nonnormahzed excited two-level system state as 

|0up(t)> = a(t) |up) . (227) 

With this the master equation ()224|1 can be replaced by the following effective Schrodinger 
equation 

i^^|0up(t)> = Heff|0up(t)>. (228) 

Obviously such an equation is technically easier to deal with than the corresponding 
master equation, because here state vectors take over the role of density operators. 
Recently similar effective Schrodinger equations were introduced to develop Monte Carlo 
wave function simulation techniques 



3.6.1 Harmonic oscillator coupled to a reservoir 

In section 12.41 we have shown, that one mode of the electromagnetic field is mathemat- 
ically equivalent to a harmonic oscillator to which in the semiclassical approximation 
we can associate the Hamiltonian H^ = fkjk °k^k (cf- (jlOOp ). where the wavenumber 
k labels the particular mode. The advantage of reservoir theory is that the particular 
structure of the reservoir does not affect the result, as long as the reservoir stays in ther- 
mal equilibrium at all times, and that perturbations of its state decorrelate immediately 
in a Markovian sense. We can therefore imagine the reservoir as a large collection of 
harmonic oscillators which are in thermal equilibrium at temperature T. For this situ- 
ation, namely the harmonic oscillator associated with the field mode, which is coupled 
to a reservoir of harmonic oscillators, we can work out the calculations in analogy to 
section 13. 6^ and we obtain the following master equation 

Pf =t [Hf, Pf] - o (^ + 1) [a^apF(^) - apF^a^ 

n ^ 2 (229) 

- 2 " [Pf Waa"^ - ^^PF^f)^] + adj. , 

where n represents the number of thermal excitations of the harmonic oscillator at 
frequency ujk- This equation shows that the decay of the mean number of quanta (a"''a) 
of the harmonic oscillator occurs with the rate k at zero temperature, whereas the 
expectation value of the annihilation operator a (which is proportional to the positive 
frequency part of the electric field mode) decays with a rate of n/2. 

4 Spontaneous emission in an optical resonator 

Up to now we have shown that spontaneous emission is a consequence of the coupling 
of the atomic or molecular system to its surrounding electromagnetic field which in the 
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last section we assumed to be the "universe" composed of a continuum of plane wave 
modes. Thus, spontaneous emission is not a generic property of the atom or molecule, 
but depends on the mode structure of the surrounding environment. Purcell |2^ found 
out that the spontaneous emission rate can be enhanced by placing an atom into a res- 
onator in which a mode is tuned to the atomic transition frequency, and Kleppner J2S] 
described the opposite case, inhibition of spontaneous emission. Spontaneous emission 
needs not be irreversible. Irreversibility is a consequence of the coupling with a Marko- 
vian Reservoir, i.e. a reservoir with a very short memory. For vacuum modes which can 
not be approximated in this way, spontaneous emission can exhibit considerable different 
properties, as in the case of long memory times, where the excitation can be periodically 
exchanged between the atom or molecule and the field J2S1 I2ZJ • 

In the microstructures discussed in section [3 the atoms or molecules are coupled to 
one, or a few, resonator modes. However, in section IHl we show that in wavelength size 
microresonators it is inevitable that the resonator modes couple with the continuum of 
modes of the external world as well. This, if the resonator cavity is not completely closed 
around the atom (or molecule), it can couple to the continuum of vacuum modes. In 
an optical resonator we therefore have to consider at least tree coupling coefficients: (1) 
The Rabi frequency Q ()165p which characterizes the atom (or molecule) coupling with 
the resonator modes, (2) the decay rate n of the resonator mode, and (3) the rate jr 
of spontaneous emission into the vacuum modes to which the system couples with. In 
wavelength scale resonators 7^ differs from the spontaneous emission rate into free space 
7fs because of the reduced mode density in the resonator. If the system couples with 
more than one field mode or when more than two levels are involved, then we have to 
consider the associated coefficients as well. 

To illustrate this concept we consider a linear resonator build with two ideal mirrors 
spaced by the distance L (Fabry-Perot interferometer). The space between the mirrors 
defines the resonator. If we neglect diffraction losses at the mirrors, the modes of the 
electromagnetic field in this resonator of volume V consist of a discrete set of standing 
waves with frequency Un = ckn = cvrn/L, where n is an integer. In section IHl we will 
show that also transversal modes exist, which we neglect in this example. If the cavity is 
small enough the mode frequencies are well separated, so that only one of them couples 
with the transition of the two level system. We label this mode frequency with Ur and 
with uq we denote the frequency of the two-level transition (nearly resonant mode). 
In reality the mirrors are not perfect and the field is diffracted, so that the mode is 
damped. We characterize the damping by the damping rate k. Damping of the field 
mode is tantamount to widening the mode frequency to a width proportional to k. In 
analogy to the free space mode density (j^ we can associate an effective mode density 
Pr with the damped mode situation in which the two-level system now radiates into. If 
we assume that the mode width is less than the separation of the mode frequency, the 
effective mode density can be approximated by a Lorentzian as 

p(uj) = -^ 5 . (230) 

Usually a resonator is characterized by its quality factor Q [211 P- 430] which is propor- 
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Figure 7: Fabry-Perot resonator with a standing wave mode of frequency Ur which is 
excited by the two- level system with dipole moment fiu- The coupling of the 
two-level system to the resonator mode Ur is characterized by the coefficient 
Q{x) which depends on the position of /2i2 with respect to the standing wave 
field distribution. The coupling coefficient 7' describes the spontaneous emis- 
sion rate into exterior vacuum modes (because of the mirrors, the solid angle 
which determines the coupling to exterior modes differs from the free space 
case), and k describes the decay rate of the field mode due to diffraction losses, 
mirror absorption and transmission. 

tional to the number of resonator round trips of a photon, and which is related to the 
damping rate k, by 

g = - . (231) 

K 

We can calculate the spontaneous emission rate in the resonator in a similar way as 
before in sections 13.41 13.51 or 13. 6t except that we have to replace the free space mode 
density (PK|) by ()23()|1 which represents the resonator geometry. For a resonator tuned 
near the atomic transition uq we obtain 



m resonance: 



7r ~ 7fs X Q 



V 



(232) 



It is interesting to observe what happens, when the resonator is tuned out of the reso- 
nance. Note that the radiating dipole is still in the resonator, and not in free space, only 
the oscillating dipole has no mode to radiate into. It can be shown that for coupling to 
the off-resonance mode the decay rate can be approximated by 



off-resonance: 7r ~ 7fs x Q 



-1 \ 
V 



(233) 



Although the situation is not utterly realistic, these equations illustrate the effect of en- 
hanced or inhibited spontaneous emission. The equations are very crude approximations 
in which we have neglected a number of things, for example that the two-level system 
interacts with free modes associated with the open sides of the resonator volume. In the 
following we will fix some of the approximations. 
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4.1 Master equation of a two-level system in a resonator 

As above we assume a two-level system in a resonator which exhibits sufficiently small 
dimensions so that the mode frequency separation c/2L is large compared to the charac- 
teristic coupling parameters of the two-level system fi, -, — . In this way the two level 
system interacts only with one resonator mode. Let us also assume that this resonator 
mode is tuned to near resonance with the two- level transition: ujr ~ cjq- In analogy to 
section Eini we now identify the small system with the two-level system combined with 
the resonator mode it couples to. Corresponding to this definition of the small system 
we define the small system Hamiltonian H^ as 

H, = Ha + H^ + H, , (234) 

which corresponds to Htot of (jl76p . The suitable mode set for describing the field in a 
resonator are the standing waves modes. For standing waves the coupling constant Vt 
becomes a function of the position x inside the resonator, VL{x) = ficos{kx) (cf. Fig. [Tj), 
where k = oj^/c and 

^ = 2/^ (-^12 ■ ^ , (235) 

with e denoting a unit vector in the direction of the electric field polarization; cf. ()165|) 
and footnote on p. EHl Transforming into a reference frame which rotates with the 
resonator mode frequency uo^ the Hamiltonian becomes 

H, = -MR3 + Hj , (236) 

with 

H, = ^ (atb + bta) , (237) 

and where 5 = u,. — ujq denotes the detuning. 

As we have mentioned before, and illustrated in Fig. [7J the small system two-level- 
system-cavity-mode is subject to two dissipative processes. The ffist one consists in 
the coupling with a cone of free space modes through the open sides of the resonator 
structure. The second dissipative process is the loss of the standing wave resonator 
mode through the mirrors to the exterior world, which is characterized by the mode 
decay rate k. Those two mechanisms are statistically independent and can therefore 
be combined by a simple addition, which in a reservoir theoretical approach looks as 
follows: The small system consisting of the two-level system and the single resonator 
mode is coupled to two thermal reservoirs that represent the cone of electromagnetic 
modes at the open resonator sides, and the mirror losses. The thermal reservoirs are 
both modeled as a continuum set of harmonic oscillators. In accordance with the Born- 
Markov approximation the coupling of the small system to the two reservoirs is described 
by a master equation in which the dissipative (non-Hermitian) part consists of the sum 
of (J22H) and (j229|) . where we now have to take into account that the two-level system 
only interacts with free space modes contained in a cone defined by the openings in the 
resonator walls. As a consequence we replace the free space spontaneous emission rate 
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7fs in ()221|1 by the smaller rate 7'. We then can derive the small system master equation, 
which is similar to ()226|) . and which at zero temperature looks like 



Psit) = ~ [HeffP, - p.Hls] + /tap.at + ybp^b^ . (238) 

where 



HeflF = Hs + Hioss , (239) 



and 



Hioss = -i/lyb^b-i^^ata. (240) 

As in section 13.31 [cf. p82|) ] we assume that initially at t = the state of the system 
is characterized by an excited two-level system and an empty resonator mode 

|^(0)) = |up,0). (241) 

As next we have to identify the possible states the system can occupy. The energy 
conserving part of the small system evolution consists in the exchange of quanta between 
the two-level system and the resonator mode. As we found out in section 13.31 the total 
number of quanta in the energy conserving process is constant - in this case here the total 
number is one quantum, or excitation respectively. On the other hand, the dissipative 
processes - that are the coupling-to-the-reservoir processes - are characterized by the 
irreversible loss of the small system excitation. In summary we therefore distinguish 
three relevant states: two states corresponding to full excitation (one quantum states), 
|up,0) and |lo, l), and the state corresponding to the total loss of excitation |lo,0) 
(zero quantum state). The evolution of the total quantum system takes place in the 
fully excited state subspace and in the zero quantum state subspace. In analogy to 
the way we derived ()227|) and (j228|) we can define the following unnormalized single 
excitation state 

|^(t)) = c,p(t)eT|up,0)+cio(t)et|lo,l) (242) 

where 5 = ujj. — uiq denotes the detuning. From ()238|) we can derive the effective 
Schrodinger equation describing the system evolution 

i;i^|vl>(t)) = Heff|v&(t)>, (243) 

where W^s is given by (j239|l . Inserting ()242j) we obtain the equations of motion for the 
coefficients Cup(t) and cio(t) 

^ = -|^c„pW-i§cio(t) (244) 

^^ = _(i5 + k/2)cM - i|c,p(t) , (245) 

which obviously is a system of coupled first order differential equations. 

In the following we will discuss two important hmit cases: the had cavity limit or weak 
coupling regime where k, 7' > fi, and the good cavity lim,it or strong coupling regime, 
where k, 7' < Vt. 



52 



4.2 Bad cavity limit (weak coupling) 



In the last section we have derived the evolution of the system, which was described in 
with ^M^i and (f^i^j) . Eq. (f^I3|) can be formally integrated: 



Clo(t) 



-i — / dt' Cup(t')e~ 
^ Jo 



{i6+K/2){t~t') 



(246) 



For partially open resonators the weak coupling regime is characterized by k, 7' > Q. 
Inspecting fl244j) we can see that Cup(t) is slowly changing if Q/2 and 7' are smaller than 
\6\ + \k\/2, what corresponds to weak coupling. In this case we can thus pull Cup(t) in 
front of the integral, which for t ^ 1/k, gives 



C\o(t) 



Substituting this in ()244|1 we obtain 



-iQ 
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Obviously this differential equation is solved by an exponential ansatz, and the occupa- 
tion probability Pup for the excited level decays as 



P, 



up 



-up 
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where again S = Ur — ujq denotes the detuning. The decay rate 7 is thus composed of 
a sum of terms of which the first can be attributed to the contribution of the vacuum 
field, and the second to a resonator contribution. In the limit where the resonator is 
completely open (free space) k — > cxd, and as a result 7' -^ 7. Thus the decay rate 7 
reduces to the Weisskopf-Wigner coefficient 7fs. 



^'=2-n-T 



4.2.1 Enhanced spontaneous emission 

Let us imagine a completely closed resonator that surrounds the two-level system. Then 
the loss mechanism which in Fig. [7|is characterized by the coefficient 7' is absent, thus 
7' = 0. Let us also assume resonance between the two-level system and the field mode, 
so that the detuning 6 = 0, and still assume weak coupling n > Q. Then the decay 
coefficient 70 becomes maximal, 70 = Vt^ /2k. With p95|) . ()231|) . ()235|) . and Aq = 27rc/u;o 
we obtain 

^ 3g 
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(252) 
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We can see that 7max has the same wavelength dependence as the approximate result 
fl232j) . By taking into account the coupling of the dipole moment with the field polariza- 
tion ()235|) we now also obtained the correct geometrical dependence. We can see that for 
resonators of good quality Q with sizes approaching V '^ \^ a. substantial enhancement 
of the free space spontaneous emission rate 7fs is expected. 

4.2.2 Inhibited spontaneous emission 

Let us consider the same two-level system in the same resonator which, however, now 
is tuned far off the two-level system resonance uoq so that the detuning can be set to 
\5\ = uq. For a good resonator with Q ^ 1 ()251|) then becomes 

1 3 A^ 

7inh ~ 7max " ^ = J^^ " ^ " ^f- (253) 

The spontaneous emission of a two-level system in a one-mode-resonator with a large 
quahty factor Q can therefore be nearly completely suppressed by detuning. 



We now have derived three different rates, ([TQoj), ()252|) . and ()253|) . by which an excited 
two-level system can decay. One might ask, how the two-level system decides which 
one to choose [221 HOI • Obviously the choice depends on the electromagnetic boundary 
conditions defining the mode structure which surrounds the two-level system, and the 
question therefore is, how does the two-level system explore its mode environment. To 
answer this question we have to consider the very initial processes of the decay, which 
thus involve time scales short compared to the inverse decay rate. Unfortunately this 
is a range, which is not consistent with some approximations we have made in our 
previous discussion, especially it is not consistent with the assumptions for using the 
coarse grained derivative ()212p . Thus, to answer this question we need a more precise 
theory. Nevertheless, let us try a qualitative argument: Let us assume the two-level 
system spontaneously starts to emit, radiating a wave packet into the mode. If the 
emission occurs into a free space mode, nothing comes back. However, if the emission 
occurs into a high-Q resonator mode, then the emitted packet is refiected at the resonator 
boundary and returns to the emitting dipole, i.e. the two-level system, where it is picked 
up by the dipole. Coded in its phase, the refiected packet carries the information about 
the geometry of the mode, as well as the state of the emitting two-level system at earlier 
times. Thus, the question now is, so to speak, what the phase is of the refiection arriving 
at the generator driving the dipole. If the refiection is received with a phase opposite to 
the actual emitting phase, then further emission is inhibited, whereas in-phase reception 
enhances the emission. We can see that a valid theory for the early emission stages must 
involve the phases of the quantum systems involved. On the other side the phase is not 
a Hermitian observable, therefore we must accept that the early stages of the emission 
are hard to divulge. 
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4.3 Good cavity limit (strong coupling) 

Strong coupling of the two-level system with the resonator mode is observed when the 
photon which is emitted by the two-level system into the resonator mode lives long 
enough so that it can be reabsorbed by the two-level system. In this case we have to 
consider the exact solution of the system of coupled differential equations (J244|) - (j245|) . 
The general solution of a system of two linear, first order differential equations is given 
by the exponential ansatz 

Cup(t) = c,pie"^* + c,p2e°^*, (254) 

where the constants Cupi, Cup2 are chosen to fit the initial conditions, for example c^pit = 
0) = 1, and the exponential are given by 



2V22/ 2 



J^ + ^ + isX-a^ 
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(255) 



For strong coupling we have Q/2 ^ 7', k, 6, and then the exponents reduce to 



ax. = -5(^ + 5 + WJ±.2. (256) 

In this expression the imaginary part is larger than the real part, so that the time 
evolution of the upper level probability Pup is characterized by oscillations at the vacuum 
Rabi frequency Q, which are slowly damped. The emission spectrum of course is no 
longer a simple Lorentzian, but a doublet of Lorentzian lines, each with a width of 
(7' -|- /€)/4, split by the vacuum Rabi frequency fi |3Tj . 

4.3.1 Dressed states 

We can generalize the situation if we merge the small system consisting of the two-level 
system and the resonator mode into a single quantum system. As we have done before, 
we associate the Hamiltonian (jl76|) to this system.^ As we have discussed before, this 
Hamiltonian H conserves the total number of small system excitations, more explicit, 
H only couples state |up,'n,) with llo, ra-|-l). Thus in a system characterized by H 
transitions only occur inside the {n + l)-quanta manifold {|up,n), |lo,ra + l)}. It is 
therefore possible to decompose H into the sum 

H = ^ H„ (257) 

n 

where H„ only acts in the (n -|- l)-quanta subspace. In this subspace the states |up,n) 
and llo, n-|- l) form a basis in which H„ can be represented by a 2 x 2 matrix. The 
eigenvalues of H„ are obtained by diagonalization as 

Ein = h [n+ ^)UJr + hVt^ 



^In the literature the Hamiltonian H176() is often referred to as Jaynes-Cummings Hamiltonian |32| . 
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Figure 8: Eigenenergy of two manifolds of dressed states (dots) and the corresponding 
bare states (lines) as a function of the two-level system transition frequency 
ujq. The anticrossings occur at resonance uq = ojr- 



where with Vt'^ we denote the n-photon Rabi frequency 



= V('^r-^0)' + fi'(^+l), 

and where the eigenstates are given as 

|2n) = — sin6'„|up,'ri) + cos^„|lo, n + l) 
|ln) = cos6'„|up,n) + sin6'„|lo, n + l) 



(259) 



in which 



tan 29r, 



VL^/n + 1 
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(261) 



The states fl260|) are called the dressed states of the two-level system, which refers to 
the picture of a two-level state dressed by the strongly coupled resonator mode. In 
this terminology the states |up,n) and |lo, n+ l\ are called hare states. Note that the 
zero-quantum manifold has a single eigenstate up, O) with eigenvalue Eq = 0. The 
dressed state energy spectrum given by ()258|) and ()259|) . and the corresponding bare 
state eigenenergies are shown in Fig. |H1 as a function of the resonator detuning. The 
energy of the bare states |up,?7,) and |lo,n-|- l) cross at resonance. However, for the 
dressed state eigenenergy this degeneracy is removed by the the interaction of the two- 
level system with the field, causing dressed states of the same manifold to repel each 
other, an effect that is often called anticrossing. 

Spontaneous emission of the two-level system occurs in the one-quantum manifold, 
where two transitions are allowed, |l,0) -^ |lo, O) and |2,0) -^ |lo, O), corresponding to 
the frequencies —6/2 + Q'^ and —6/2 — Q'q. At resonance 6 = the separation becomes 
equal to the vacuum Rabi frequency Q; cf. Fig. IHl 
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Figure 9: Spectrum of the spontaneous emission doublet in the one-quantum manifold. 

5 Stimulated emission in an optical resonator and 
thresholdless lasing 



In section 14.21 we discussed the modification of the two-level system decay when a sur- 
rounding resonator reduces the density of field states. Equation ()25H) was derived for 
a two-level system that interacts either with an empty, or with a one photon, damped 
resonator mode. As a consequence, three system states had to be considered: |up,0), 
|lo, l), and |lo,0). In (12511) we distinguished two contributions, the first is attributed 
to the mode density, and the second, which was characterized by the vacuum Rabi fre- 
quency Q, to the field. Stimulated emission occurs when the excited two-level system 
interacts with a mode which contains several photons. In analogy to section 14. ![ and 
observing the approximations for the Fermi golden rule, cf. ()196|) . we can derive the fol- 
lowing result for the decay rate of a two-level system in a damped resonator, interacting 
with an n-photon field 

where with fi^ we denote the n-photon Rabi frequency f2^ = y5M-TF(n + l); cf. ()259|) . 
We can see that for a resonator tuned to resonance {6 = 0) the stimulated emission decay 
rate scales with the number of photons n interacting with the two-level system. Equation 
(|262j) can be written as 

7st = 7o(r^ + l), (263) 

where 70 stands for the spontaneous emission rate which is modified by the resonator, 
cf. fl251|) . For fluorescent dyes we can assume that the lower laser state is rapidly depop- 
ulated, and that nonradiative processes, inversion saturation and collective spontaneous 
emission are absent. With this assumptions the number of photons in the resonator 
mode can be expressed as 

n = 7o(ra + 1)A^ - ku , (264) 



57 



where N denotes the inversion A^ = A^'up — A^io, that is the difference between the popu- 
lation of the upper laser level iVup and the lower laser level Nio. The rate by which the 
inversion N is changed depends on the rate with which the upper level is populated, 
that is the pumprate W, and by the rate of depopulation given by the photon emission 
rate, thus 

N = W --fo{n + l)N (265) 



The steady state solution of the system of differential equations ()264p and ()265|) has the 
simple form 



and 



W 
n= — (266) 

K 

^ = -iZ^ (267) 

7o(W^ + /t) 

From (j266|) we can see that the light output n increases for any pump rate ly in a linear 
way, thus the laser operates without threshold. This is a phenomenon called thresholdless 
lasing or lasing without inversion. This is the consequence of the exclusive coupling of 
the laser active two-level systems to one single resonator mode. On the other hand, 
the inversion N approaches the limit of Nt = n/'Jo when the pump rate W increases. 
This limit can be interpreted as the inversion at which the photon emission changes 
its character from spontaneous emission to stimulated emission. This is the transition 
which traditionally is associated with the onset of lasing. Note that this phenomenon 
occurs when a collection of two-level systems interacts with a single field mode, which is 
thus different from the interaction of only one two-level system with a single mode, i.e. 
a one atom laser ^^. In summary, the conditions for inversionless lasing as expressed 
in (j266j) and (j267|) are that the collection of two-level systems interacts only with a 
(damped k) single mode field. In particular there are no open resonator side walls 
through which the two-level system can couple to the mode continuum (characterized 
by 7'); cf. Fig. [7| Of course there is a continuous transition in the threshold behavior 
from the open resonator case shown in Fig. [7|to the closed resonator discussed here [21] 
and thresholdless lasing was observed in a dye solution placed a plane-plane resonator 
in which the mirrors were spaced by half an emission wavelength [33] . In section [7| we 
will discuss threshold reduction in dye loaded molecular sieve microlasers. 

6 The mode structure of microresonators 

6.1 The concept of dielectric constant, dielectric interfaces 

Maxwell's equations are the foundations on which both classical and quantum optics are 
built. The possibility of creating a double-slit diffraction pattern with a single photon 
- a famous manifestation of the wave-particle dualism - is a result of the fact that 
photons occupy the modes of the classical electromagnetic field, cf. the definition of the 
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bosonic operators in terms of normal modes in sectional These modes can be discrete in 
frequency, as they are in lossless closed resonators, or have a continuous spectrum as is 
the case in scattering experiments such as the diffraction example [36j; there, the whole 
information about the diffraction process is encoded in the modes themselves because 
they must satisfy the boundary conditions on the diffraction slits. 

Boundary conditions are a way of introducing rapid spatial variations of a medium 
into the macroscopic Maxwell equations without having to give up the concept of a 
dielectric constant or permittivity, by which the microscopic properties of the medium 
are taken into account in a very efficient mean- field manner |SZ]. Even when details 
of the quantum- mechanical light-matter interaction are of interest, such a mean- field 
approach is a good starting point in order to define an appropriate modal basis set in 
which to expand the relevant matrix elements. This applies to dielectrics of infinite 
extent without losses [2H] or with losses ^Hl, but especially when fiuorescent or laser 
emission in the presence of a microcavity are concerned. 

The dielectric "constant" is moreover a function of field strength, if the polarizabil- 
ity - which describes the microscopic matter-field interaction - depends nonlinearly on 
electric field. When such nonlinearity occurs in combination with boundary conditions, 
the relative importance of the two depends on the size scales and field strengths in- 
volved. At low pump levels, phenomena such as vortex formation which has long been 
known in nonlinear media 03], can cross over to linear vortices ^^, an example of 
which are the whispering-gallery modes which will be discussed in section 16.31 In small 
cavities, boundary effects become dominant, and nonlinear effects can be understood as 
interactions between modes of the linear cavity 03] . 

As we discussed above (cf. section 0}, a well known microcavity effect is due to Pur- 
cell j231 who argued based on Fermi's Golden Rule that the Einstein A coefficient for 
spontaneous emission can be enhanced in a microcavity due to its highly peaked density 
of modes. For this perturbative approach to the quantum electrodynamic problem of 
spontaneous emission, the modes into which photons are emitted are determined com- 
pletely by Maxwell's equations with the boundary conditions defining the cavity. In 
the strong-coupling regime, this ceases to be correct when the radiant matter and light 
field are mixed in comparable proportion in the eigenstates of the total system, as can 
occur, e.g. in cavity polaritons in a quantum well microstructure jJU]; or for the Rabi 
oscillations of an atom in close proximity to a dielectric microsphere [H] . Nevertheless, 
the first step in all these cases is to determine the modes of the electromagnetic field for 
the resonator goemetry at hand. 

The conclusion from these preliminary remarks is that especially in microcavities, an 
understanding of the effects of the resonator geometry on the field distribution, neglect- 
ing nonlinearities of the medium, is of central importance. Even in the conceptually 
simple problem that remains, we shall see how the well-understood fundamental equa- 
tions of electrodynamics lead to solutions that are at present only partially understood, 
in the sense of predicting their dependence on system parameters, or even giving con- 
ditions for the existence of certain solutions - an important example again being the 
whispering-gallery modes. 
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6.2 Fields at dielectric interfaces 

6.2.1 Matching conditions 

We assume that all fields have the stationary time dependence e~*'^*. Then Maxwell's 
equations become 

1 rlT-f 

V X E = 7- = ikH (268) 

c ot 

V xH = - ^— = -ikn^E , (269) 

c dt ' ^ ^ 



where the wave number is given by 



k = u/c. (270) 



Since we want to illustrate the effects of boundaries, let us make the further simplification 
of considering the refractive index n to be piecewise constant, but not necessarily real. 
Combining these equations, we obtain the wave equations 

V X V X ^ = {nk)^E , (271) 

V X V X if = {nkyH . (272) 

Since charge density can only appear at the surface of the dielectric, we have V ■ E = 
in each domain of constant n. In these regions, (|271|) therefore becomes 

-V^E = {nkfE . (273) 

The conditions to be satisfied at a dielectric interface are deduced from the requirement 
that no current flow is possible along the interface |12l! implying that the tangential 
components of E and H must be continuous. At vertices or edges, however, the tangent 
is undefined. In this case we must invoke the additional requirement that the energy 
contained in any volume element of the fields should not diverge near such singularities 
of the surface |121 • 

These conditions do not constitute true boundary conditions in the traditional sense, 
but are instead matching conditions. The field on one side of the interface is determined 
by the field on the other side, and only if we already know the latter, can the former 
be obtained by solving a boundary-value problem. For some simple problems, it is in 
fact easy to eliminate the field on one side of the interface, if one knows its form a 
priory. The simplest example is a plane wave in air, impinging with wavevector fc on a 
planar interface with a lossless dielectric of refractive index n. The knowledge that the 
transmitted wave in the medium is again a plane wave allows us to straightforwardly 
obtain Snell's law of refraction and Fresnel's formulas for the reflectivity. The latter 
depend on polarization, which in this case can be chosen either perpendicular to (TE) 
or in the plane of incidence (TM). With this we then can obtain decoupled scalar wave 
equations. The essential difference between TE and TM polarizations is that the latter 
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exhibits the Brewster angle xb at which perfect transmission occurs. Independently of 
polarization, Snell's law relates the incident angle xo (which we measure with respect to 
the surface normal), to the transmitted angle x of the plane wave in the medium by 

sin Yn 

sinx= —. 274 

n 

Note that for large n, transmitted waves in the dielectric are allowed to propagate only in 
a progressively narrower interval around the perpendicular direction x ~ 0- Conversely, 
the critical angle Xc = arcsin(l/n) therefore defines the "escape cone" for plane waves 
inside the material. Total internal reflection prevents all waves with x > Xc from escaping 
to the optically thinner medium, i.e. the condition for confinement is 

sin X > - ■ (275) 

n 

6.2.2 Impedance boundary conditions 

The fact that high index contrast leads to a narrow escape cone is at the heart of a large 
body of literature summarized in Ref. ^3. The aim is to replace the continuity require- 
ments for the fields at a dielectric interface by approximate boundary conditions called 
(normal) impedance conditions. There, an approximate knowledge of the field on the 
high-index side allows one to eliminate it from the problem. Yet, the method becomes 
unreliable at low refractive-index contrast and at non-planar interfaces; in particular, the 
Brewster effect is not correctly reproduced unless additional corrections are introduced 
which, however, do not permit a simple physical interpretation, because they involve 
derivatives of the fields that are of higher order than the original continuity conditions 
for the fields themselves. We shall therefore make no use of impedance boundary condi- 
tions, except to point out that they reduce to the familiar Dirichlet boundary conditions 
(vanishing fields on the surface) if the index contrast becomes infinite, as would be the 
case in an ideal metal. 

6.3 Scattering resonances versus cavity modes in the dielectric 
cylinder 

One way to probe the interaction of dielectric bodies with light is by elastic scattering. 
Atmospheric phenomena such as the rainbow, the halo, or the glory arise from light 
scattering |3H] , and in fact atmospheric science relies on scattering experiments as diag- 
nostic tools. When scattering experiments are carried out with high spectral resolution, 
a ripple structure is observed in scattering cross sections, which cannot be understood 
in a purely ray-optics framework. These ripples are resonances which occur when the 
incident light couples to long-lived cavity modes. In order to illustrate the relationship 
between resonances and modes, we consider here the example of a dielectric rod with a 
circular cross section. 

Dielectric cylinders are of great practical interest because they are the archetypical 
model for an optical fiber. In the context of this review, the approximate cylindrical 
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symmetry of the molecular sieve microcrystals makes it desirable to establish some 
fundamentals of cylindrical systems. The modes of an optical fiber can be divided 
into two classes: guiding modes and "leaky" modes. As is well known, guided modes are 
responsible for the ability of silica fibers to carry long-distance communication signals, 
corresponding in the ray picture to zigzagging trajectories traversing the fiber core in an 
almost planar motion - the plane of propagation coincides with the cylinder axis. The 
confinement mechanism in the ray picture is total internal reflection at the interface to a 
lower- index medium, either the cladding of the fiber or - as we will assume for simplicity - 
the abrupt interface with air. However, over short distances a significant power transport 
can take place through the leaky modes as well, which must be taken into account in 
any comprehensive treatment of optical waveguides [IH]. It is also known jlHl that 
one has to distinguish between tunneling leaky modes which are confined by frustrated 
total internal reflection, and refracting leaky modes whose attenuation is even larger 
because they correspond to rays that violate the condition for total internal reflection. 
Leaky modes correspond to rays spiraling down the fiber in a helical trajectory, and are 
therefore also called "spiral modes" [HUj . 

Of these three types of modes - guided, tunneling, and refracting - for propagation 
along a fiber, only the leaky ones remain if the radiation is incident with propagation 
vector perpendicular to the fiber axis. This is the situation we shall focus on, because 
the applications we have in mind are not intended for power transport but for power 
storage, i.e. resonators. The corresponding ray trajectories are then confined to a plane 
in which they perform a circulating motion. An example where this mode structure 
has in fact been observed directly is a laterally structured cylindrical VCSEL |51j where 
such a mode supports lasing action at an unexpectedly low pump power. Because the 
circulating mode structure is analogous to the acoustic "whispering-gallery" phenomenon 
in which sound waves cling to the curved walls of certain buildings, modes with a ring- 
shaped intensity distribution are commonly termed "whispering-gallery (WG) modes". 

The range of material parameters relevant to the design of dielectric resonator struc- 
tures is increasing continuously, as novel materials enter device applications. For glasses, 
refractive indices between n ^ lA and 2 are available (glass for application in fiber am- 
plifiers has n ^ 1.8). Semiconductor materials extend this interval to even larger indices, 
while organic compounds border on the lower end of the index range. The fundamental 
wave equations for elastic light scattering in the case of no propagation along the axis of 
the dielectric cylinder are much simpler than for arbitrary oblique incidence, and conse- 
quently treatments of plane wave scattering from dielectric cylinders at normal incidence 
can be found in many textbooks jlHl 113 Ell- Because of the ever increasing range of 
applications, we review them here, emphasizing, however, what in the literature is usu- 
ally missing, namely the aspect of the relation between resonances and cavity modes. 
One important piece of notation that we introduce here is conspicuously absent from 
the classic texts: the formulation of the light scattering problem in terms of S-matrix 
theory, as it is used in quantum mechanical scattering. Investigations of leaky cavities 
may well profit from a more unified scattering-theoretical terminology. 
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6.3.1 Metastable well in the effective potential 

To understand the origin of long-lived whispering gallery mode resonances we examine 
Maxwell's equations for an infinite dielectric cylinder. Following an argument used by 
Johnson j33] for the dielectric sphere, we can rewrite (j273j) in a form similar to the 
Schrodinger equation of quantum mechanics, 

-W^E + k\l~n')E = eE . (276) 

This shows that dielectric regions {n > 1) correspond to an attractive potential well in 
the quantum analogy, except that here the potential is itself multiplied by the eigenvalue 
P. 

The reason why extremely long-lived resonances are created when the scatterer has 
rotational symmetry is that after separation of variables in cylindrical or spherical coor- 
dinates, ()276|) gives rise to a radial equation with a repulsive potential term due to the 
angular momentum barrier, as well as the attractive term just noted. For the case of a 
dielectric cylinder studied here the resulting equation reads 



d' 1 d 
dr"^ r dr 



where the effective potential is 



E{r) + Kff(r) E{r) = k^ E{r) , (277) 



2 

m 



V,^{r) = e{l-n') + '^ + kl. (278) 






The additional centrifugal potential term has appeared as a consequence of conservation 
of the z-component of angular momentum, and the offset k1 results from the conserva- 
tion of the z-component of linear momentum. As inidcated above, we focus on planar 
propagation perpendicular to the 2;-axis of the cylinder, so that fc^ = 0, and the inci- 
dent wave has its /c-vector in the x — y plane. The problem then becomes effectively 
two-dimensional and we will use the two-dimensional polar coordinates r = \Jx'^ + y^, (j). 
The radial "potential" which results from the sum of the attractive well due to the 
dielectric and the repulsive angular momentum barrier is shown in Fig. EH cf. also Ref. 
jS3] . At nonzero angular momentum one sees that a metastable well is formed within the 
dielectric. The electric field inside is separated from the propagating field outside by a 
tunnel barrier, and one may expect solutions of the wave equation in which the intensity 
inside the well is exponentially larger than outside the well. In order to complete the 
formal analogy to the one-dimensional Schrodinger equation we make a substitution of 
variables to eliminate the first derivative term in the radial equation ()277|) . This can be 
achieved by introducing a new coordinate ^ := ln(A;r); the resulting equation is 



;^ +,=(?) 



E{r) = , (279) 



where g(^) = (?T,^exp(2^) — m^)^/^, and g(^) is the effective wave vector of the rescaled 
problem. In order that there are any resonant solutions, kR must be sufficiently large so 
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Figure 10: Effective potential picture for whispering gallery resonances of the cylinder; 

kmin = m/{nR) and k^nax = m/R. 



that q is real for some value of r within the dielectric. The largest value of the argument 
of the square root within the dielectric occurs when r = R, so we deduce the condition 



m 
nR 



;280) 



On the other hand, if infinitesimally outside the dielectric where tt, = 1, g is real, then 
there is no region of exponential decay, and no sharp resonances exist. The maximum 
value of k such that q is imaginary just outside the dielectric is given by 



m 

"'max "^ 



(281) 



Consequently, we expect narrow resonances of angular momentum m, broadened only 
by tunnelling decay (evanescent leakage), for wave vectors satisfying 



m , m 
—^ < k < — 
nR R 



or equivalently 



1 m 

n nkR 



< 1 . 



(282) 



(283) 



We will label these metastable well states by a discrete index v which corresponds to 
the number of radial nodes. If the refractive index is changed, the barrier top in Fig. ^J 
remains the same, but the well grows deeper so that the lowest allowed k will give rise 
to a narrower resonance. Later we will show that essentially all of the metastable states 
satisfying ()282j) are discrete in the sense that their spacing in k is larger than their decay 
width through the barrier. In this way they will give rise to isolated resonances in the 
scattering functions which are well set apart from the background. 

It is worth quoting here the main features of the numerically observed resonances in 
plane-wave scattering off a dielectric cylinder at normal incidence, as given in Ref. 
(but using our notation): 
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. . . their widths decrease as m [angular momentum] increases for a given u 
[radial node number], and their widths increase as u increases for a given 
m. Also, as the index of refraction [n] is increased, the positions of the 
resonances shift to lower kR, and their widths become narrower. Resonances 
having relatively narrow widths occur in the range x to nx. 

All of the above phenomena can already be understood qualitatively by inspecting the 
effective potential picture. 

The inequality (j283j) can be brought into a very suggestive form which can be justified 
by considering the short- wavelength limit in the WKB approximation |H^ : 

TTl 

sinx:=^^, (284) 

then (j^ leads to 

- < sinx < 1 • (285) 

n 

The right inequality is trivial, but the left relation is precisely the condition for total 

internal reflection which must be satisfied by the angle of incidence x oii the high-index 

side at an interface to air, cf. (j275p . At this point we have provided no reason why x, 

as defined here, should be the classical angle of incidence. Let us try a non-rigorous 

argument: If we interpret the ray trajectories as classical particles with momentum 

H{nk) inside the resonator, and angular momentum L = hm, then the classical relation 

between angular and linear momentum would imply L = f x p---^ L = R (hnk) sin x '^ 

m = Rnk sinx, where sinx arises from the vector product of linear momentum and 

radius vector r, which in turn has magnitude R at every reflection from the boundary. 

At such a reflection, x is then just the angle of incidence. 

The introduction of h into the present optical context is somewhat arbitrary. The 
question is, does the optical field in the cavity truly acquire angular momentum. We 
know that photons, as massless fundamental particles, carry a unit spin angular mo- 
mentum, but no orbital angular momentum in the physical sense. Therefore, when we 
refer to m as an angular momentum quantum number, we invoke a formal analogy that 
arises between the mathematical structure of bound photonic states and bound quan- 
tum systems: There are, in our case, two degrees of freedom which are associated with a 
finite region of space, giving rise to a discrete frequency spectrum. The modes are then 
labeled by as many discrete indices, or "quantum" numbers as there are confined degrees 
of freedom. For a rotationally symmetric problem, the indices, or quantum numbers, are 
simply the number of radial and azimuthal nodes of the wave function. The azimuthal 
quantum number is what we call m. Its relation to the angle of incidence is rigorously 
derivable from the eikonal (or WKB) limit [54j. 

There is an additional complication to this argument if the cavity is not closed. Due 
to the openness of the system, it is then not clear how to define the photons in a resonant 
mode. However, the cavity can approach the closed limit as n ^ oo. Before we discuss 
this problem, we will elucidate the concept of modes in the leaky cavity setting, where 
discrete resonances appear but have a finite width in k. We will return to this issue 
under the heading of quasihound states in section] 
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For this moment let us nevertheless dwell on ()285j] . which we can use to interpret 
some results. For example, it indicates that states at the bottom of the well correspond 
to a ray motion tangential to the boundary, whereas states at the top correspond to 
rays colliding with the boundary at exactly the critical angle. As next we turn to a 
more detailed discussion of the actual wave solutions for the dielectric cylinder where 
the incident wave propagation is normal to the axis. 

6.3.2 Matching conditions for TM polarization 

Our general considerations so far have not addressed the boundary conditions at the 
interface between dielectric and vacuum (or air) - conditions which also depend on the 
polarization of E. The arguments given above for the existence of narrow resonances 
are, however, independent of the boundary conditions, because we rely only on the 
fact that a tunnel barrier is formed where the field decays. Thus there should exist 
resonances corresponding to two different polarization states. Henceforth we will focus 
on the transverse magnetic polarization for which the boundary conditions lead to the 
simplest matching conditions for the electric field inside and outside the dielectric. For 
this polarization state we assume an incident plane wave of the form E = E{x,y) e^, 
where the unit vector e^ points parallel to the cylinder axis (and hence the magnetic 
field is transverse to the cylinder axis). From now on we will focus on the scalar function 
E{x, y) = E{r, 0) which is the amplitude of the electric field in 2;-direction. Since a pure 
dielectric is an insulator, there is no current flowing on its surface in response to the 
incident field and H is continuous everywhere (neglecting any variation in the magnetic 
permeability). From ()268|) . we derive the azimuthal component of the magnetic field, 
which in polar coordinates has the form 

i7, = - [V X Ee- ], = -^ ^ . (286) 

Therefore the radial derivative of E{r, 0) is continuous at the dielectric surface. Since 
the tangential component of E is always continuous at such an interface, we arrive at 
the simple matching conditions that both E and dE/dr are continuous at the cylinder 
surface. TE polarization [H along e^) does not result in these familiar requirements 
|46| ITF] , however one can still obtain a scalar wave equation, which can be treated with 
a generalization of the methods used below. Most conveniently, one writes the wave 
equation for the magnetic field H, which is again continuous, but exhibits a jump of 
magnitude n^ in the normal derivative at the interface. The effective potential picture 
with purely tunneling escape for the whispering gallery modes remains unchanged by 
this, since our arguments above did not involve polarization. Note that for TE polar- 
ization the electric field is in the plane of incidence. Transmission out of the dielectric 
becomes unity at the Brewster angle which neglecting the finite curvature of the interface 
is given by 

sinxB = ^==. (287) 
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The whispering gallery mode criterion sinx > 1/n is clearly not modified by this effect 
since sinx^ <'\-/n. 

An important point is that the same matching conditions remain valid even when 
the cross section of the cylinder is non-circular, because for TM polarization E remains 
tangential to the surface. As long as we consider only convex surfaces, the continuity of 
dE/dr in polar coordinates implies the continuity of the normal derivative of E. These 
boundary conditions are identical to those for quantum scattering from a potential well. 
Hence, we are allowed to use the terms familiar from quantum theory, such as tunneling, 
as we did in the previous sections. However, as was already noted with ()276|) . it should be 
kept in mind that the analogy is incomplete, because the dispersion relation in quantum 
mechanics isuj = h k"^ /{2m) and not uj = nek asm optics. The results presented here can 
therefore not be obtained by simply copying known quantum mechanical calculations. 
An important difference between the quantum and optical wave equations can be read 
off ()278p . where the effective potential well is seen to become progressively deeper with 
increasing k. As a consequence, above barrier reflection remains signiflcant even for 
large /c, whereas it would become negligible in an attractive circular quantum well at 
high energies. This effect is just the well known classical Fresnel reflection at a dielectric 
interface, which independently of the wavelength has the value R= {1 — nY /{I + n)^ at 
normal incidence, i.e. sinx = m/{nkR) — ^ 0. 

6.3.3 Resonances in elastic scattering 

The wave equation for TM (or TE) polarization as defined above is of the scalar form 

V^V^ + n2fc2^ = 0, (288) 

where if) denotes the electric field for the TM case, and the refractive index ra is a constant 
different from 1 inside the dielectric. The wavenumber k is not position dependent 
because it enters the problem simply through the definition uo = ck for the harmonic 
time dependence. All the scattering functions for the elastic scattering of light (i.e. 
scattering at the fixed frequency uj) by a dielectric cylinder can be calculated if we know 
the scattering states 



AmJm{nkr) {r < R) 

H^^\kr) + SmH^\kr) {r > R), 



E^{r) = <( ,^(2),, , , ^ „(i)., . . . „. (289) 



for all m {Jm and Hm denote the Bessel and Hankel functions, respectively). These 
are solutions to the radial equation (|277|) . assuming that the refractive index is unity 
outside the cyhnder and n{r) = 1 + (n — 1)6(-R — r) inside. The amplitudes Am, S^ are 
determined by the matching conditions for E and its radial derivative E', 

AmJminkR) = H^\kR) + SmH^\kR) 
AmnJ'minkR) = HJ^^\kR) + SmH^^\kR) , (290) 

where, the primes denote differentiation, e.g. J'{x) = dJ/dx. This can be solved for 
the scattering amplitude Sm as a function of the size parameter x = kR. In terms 
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of quantum mechanical scattering theory, Sm is a diagonal element of the scattering 
or S-matrix which for the rotationally invariant scatterer is diagonal in the basis of 
angular momentum states m. Using these solutions, one can for example calculate the 
scattering properties for an incident plane wave. A plane wave can be decomposed into 
a combination of Bessel functions 1521 as 



-ikx 



J2 i-ire'^^^Jmikr). 



(291) 



m=— oo 



This then determines the coefficients with which the partial waves ()289p are superim- 
posed to get the full solution. The conventional representation of this solution is 



E(r^ 



g gim0 (^-i)m^^ Jm{nkr) (r < R) 

Jm{kr)-brnH^\kr)\ {r>R) 



m=— oo 

oo 



_ m=— oo 

r(l) , tt(2) 



(292) 



Using Jm = ^{Hm + Hm ) in the incoming wave, this becomes for r > _R 

^ oo 

m=—oo 

SO that we can extract the relationship between S- and the 6-coefficients, 



(293) 



(294) 



Hence the hm are analogous to the quantum mechanical transition matrix T, defined by 
S = 1 — 27ri T. The S-matrix satisfies an additional condition as a consequence of fiux 
conservation, namely ^1] 

S^S = 1. (295) 

This unitarity relation also holds for deformed resonators, where S is not diagonal. It 
provides an independent equation with which to check the accuracy of the wavefunction 
matching. For the circular cylinder, unitarity just implies that all Sm have modulus one. 
For numerical calcultions borrowing this formalism from quantum-mechanics is therefore 
an added advantage. 

If a plane wave is incident on the cylinder, the scattered intensity has the angular 
dependence 1^1 ITF] 



/w« 



E "" 



^im(f> 



(296) 



where is the angle with the beam. From the matching conditions ()290p . one finds that 
the expansion coefficients can be written in the form 

1 



P. 



l + i(3m ' 

nJ^_i(nx) Yn,(x) - :im{nx)Ym-i{x) 

nJ„_i(nx) J^(a;) - J^(nx) J„_i(a;) 



(297) 
(298) 
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where we have used the abbreviation 



X = kR 



(299) 



and the function Ym(x) is the Bessel function of the second kind. This provides an 
exphcit formula for the scattering intensity I{4>), an example of which is plotted in Fig. 
im As is seen in the figure, the scattered intensity shows rapid resonant variations. 
The resonances are not simply Lorentzian peaks, but are of the Fano shape jSHj with 
varying values of the Fano asymmetry parameter. This asymmetric shape is due to 
interference between the nonresonant scattered waves and the resonant scattering. As 
long as the resonance is sufficiently isolated the resonance widths can be extracted by 
fitting a Lorentzian to the Fano function. It can immediately be observed that the 
resonances have widely varying widths corresponding to how close they are to the top 
of the barrier in the effective potential well. 




Figure 11: Irradiance I which is scattered off a circular cylinder with refractive index 
n = 2 at = 50° with respect to the direction of the incoming wave which is 
assumed to be plane and TM polarized. 



The decomposition ()297|) is analogous to the case of a scattering sphere [5Bj, and in 
the same way one observes that intensity maxima occur near the value of k at which 
one of the bm has a maximum (if there is no interference with other nonresonant waves 
then it will be at exactly that value of k). This in turn occurs when Pm = 0. Thus, we 
set (I298|l to zero and obtain an implicit equation for the resonance positions x = kR in 
each angular momentum channel m: 



n3rn-iinx)Ym{x) = 3m{nx)Ym-i{x). 



(300) 

It should be noted that the denominator of ()297|) is nonzero only for the regime of below 
barrier resonances, x <. Otherwise, additional resonances arise whenever j3m becomes 
infinite, i.e. when 

nJrn-l{nx) Jm(x) = Jm(nx) Jm-l(x). (301) 

Above-barrier resonances thus occur if either ()300j) or ()3()1|1 is satisfied. Yet, these 
resonances will typically be so broad that they do not give rise to isolated peaks in the 
scattering functions. 
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6.3.4 Quasibound states at complex wavenumber 

Having obtained the condition for when a resonance is excited in elastic scattering, we 
now turn to a complementary approach in which no incident wave is present. Resonances 
can also be excited without any incoming wave at that frequency, e.g. if the light is 
generated inside, such as in a laser. Also the specific lineshape of a resonance in elastic 
scattering depends on the nature of the incident wave and is not a property of the 
scatterer alone. It is thus useful and important to consider an alternative definition of 
the resonant state which is independent of the manner of its excitation. We will first 
introduce the concept of the quasibound state. We will discuss this concept using the 
example of the dielectric cylinder we have have considered above, and then turn to the 
more general significance of such states. 

Solutions of the wave equation with no incoming wave are termed quasibound states 
(or in the recent literature on resonances of dielectric spheres, quasinormal modes [5^ , 
hinting at the fact that they are not an orthonormal set in the usual sense). Due to 
the unitarity of the ^-matrix for real k, there can be no such solutions for real values 
of k. Instead quasibound states are directly connected with the complex poles of the 
S"- matrix (or in our simple case the matrix element Sm), the real part of k giving the 
resonance frequency, and the imaginary part giving the resonance width. This is because 
the condition for having a solution to the matching equations with no incoming wave is 
just the condition for a pole of the S'-matrix 15^: If in (j289|) we give the incoming wave 
{Hm an arbitrary amplitude Aj, then the outgoing wave will have amplitude 

Ao = SraA. (302) 

If we now set Ai ^ 0, there can be a finite outgoing radiation if at the same time S^ 
has a divergence, requiring us to tune k to the complex value where Sm has a pole. 
Thus the complex wavenumber of the quasi-bound state is determined by 



Am Jm{nkr) (r < R) 
SmH^\kr) {r>R), 



Em{r) = <( ^ „{1),,.. ^ ,__ . „^ (303) 



and the accompanying condition for continuity of the derivatives. These matching con- 
ditions for E can be satisfied only for a discrete set of complex k at which there exists a 
nontrivial solution to the homogenous part of the linear system ()290p . The latter occurs 
only at the zeroes of the determinant 



D 



JminkR) -HJn\kR) 
nJ'^inkR) -H^^'ikR) 



(304) 



which leads to 



nJ'^{nkR) HJ^\kR) = Jm{nkR) H^m^'{kR) . (305) 

Using the recursion relations for the Bessel functions to eliminate the derivatives, we 
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obtain the resonance condition 

n Jm~i{nkR) — 



m 



Jm{nkR) H^^\kR) 



(306) 



J 'n 



m 



nJr, 



nkR) Hi^l^{kR) 



nkR 
(nkR) ' 

i{nkR) H^^\kR) = JminkR) H)^i-^{kR) . (307) 

This can be solved numerically to find the real and imaginary parts of kR at which a 
metastable state occurs. Equation ()307j) is a complex equation for a complex variable. If 
we restrict kR to be real, then the real and imaginary parts of the equation are precisely 
()3()()j) and ()3()H) which we obtained for the scattering resonances. This shows the direct 
relationship between the quasibound states and the scattering resonances: For real kR, 
()307p cannot be exactly fulfilled, but the closest one can get to satisfying it with real 
kR is given by the scattering resonance conditions. 

The exact results obtained numerically for the real parts of the quasibound state 
wavevectors are summarized in Fig. El For each m, there exists an infinite sequence of 
discrete points in this grid, corresponding to an increasing number of radial nodes in the 
effective potential well. The quasibound states are hence labeled by angular momentum 
and radial quantum numbers (which we call m and v) , as in the closed circle. We have 
drawn in Fig. El a line of slope m/kR = 1 corresponding to the condition sinx = ^; 
resonances to the right of that line are true whispering gallery modes and we expect 
their widths to be much less than the resonance spacing, with the narrowest resonances 
corresonding to the points with the largest distance from this line. 




10 20 30 40 
Figure 12: Resonance positions in the circle of %fractive index n = 2 with TM polariza- 
tion for the lowest 44 angular momentum numbers m. Each dot corresponds 
to one resonance, and the solid line kR = m represents the dividing line be- 
tween broad and narrow resonances, with narrow widths expected below it. 
The horizontal dashed lines enclose a kR interval in which we count a total of 
9 resonances below the critical line. Compare this observation with the scat- 
tered intensity in Fig. ^2 The three narrowest resonances there correspond 
to the three lowest kR points at m = 17, 18, 19. The closer the points are to 
the critical line, the broader the resonances get. 

Note that the spacing in k of resonances for a given value of m is roughly constant 
above the line m/kR = 1, whereas the spacing of resonances below the line increases 
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and hence is not independent of kR for a given m. The latter are the whispering gallery 
mode resonances of primary interest to us. 

The uniform spacing of the above barrier resonances can be demonstrated by using 
the large argument expansions of the Bessel functions in ()307p . This is justified at large 
enough size parameter nkR ^ m, noting that nkR > m always holds. One obtains 



n cos I nkR — —m — j H)^\kR) 

= cos i nkR — —m — —] H^_i{kR). 
Noting that cos(a — 7r/2) = sin a, this simplifies to 

1 H';l\kR) 



tan nkR m 

2 4 



n 



(30^ 



(309) 



An important limiting case is n ^ cxd, for which the large argument expansion always 
applies if m and kR are fixed. Then the righthand side becomes negligible, so that the 
resonance position is given by 
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(310) 



where the nonnegative integer v is the radial quantum number. This is just the limit of a 
closed system with Neumann boundary conditions, J^{kR) = 0, in the limit nkR ^ m, 
and it clearly leads to a purely real kR, i.e. to truly bound states. 

On the other hand large n is not always necessary to get an approximate result: one 
obtains the same expression for the resonance position independently of n, provided the 
quotient of the Hankel functions is purely imaginary. This is precisely the case for above- 
barrier resonances, where kR > m, and hence we can use the large argument expansion 
for the Hankel functions, too. The righthand side in ()309|1 then becomes equal to —i/n, 
and one can solve exactly for 
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Thus the real and imaginary parts of kR = x + iy can be given analytically in the limit 
kR > m: 
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(312) 
(313) 



Since the righthand side in ()309j) has no real part when kR -C m, it does not affect the 
condition for the real part of kR. The approximations made here are thus valid either 
for large n or for above barrier resonances. 
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Of greatest interest to us, however, are the long-hved resonances with kR > m. For 
this case, the ratio of Hankel functions will acquire an m- and /c-dependent real and 
imaginary part, which will emerge in both resonance positions and widths. Approximate 
explicit expressions for resonance widths and positions can be obtained, but are not 
generally valid for the entire range of interesting parameters n, kR, m |5S- - therefore, 
we can for now be content with the simple limiting cases considered here. They are 
already of practical use because a numerical solution of the implicit equation ()307|) for 
the real and imaginary part of k at given m and n is made more efficient by using the 
approximate resonance positions from ()312|) as an initial guess in finding the roots of 
pTTjl . 

6.4 Quasibound states in lasing and fluorescence 

The previous considerations have brought us to the quasibound state concept via the 
resonances in elastic light scattering. Even if we allow the cross section of our dielectric 
cylinder to be deformed, TE and TM polarization can still be decoupled (see the remarks 
in section lB.3.2j) . provided we continue to restrict ourselves to propagation normal to the 
cylinder axis. After we clarified the role a complex wavenumber k can play in the 
original wave equation, we deduce a rather general statement about quasibound states. 
To arrive at the time- independent wave equation ()273|) . we assumed a monochromatic 
time variation e"'^*^*, with k = u/c according to (|27(J|) . With the quasibound state 
boundary condition (no incoming wave), this same wave equation admits solutions only 
at discrete complex k = u/c — ij/c, where we have split the frequency into its real and 
imaginay part. Returning to the ansatz for the time dependence, this implies that the 
quasibound state decays in time as exp[—iujt — •yt]. 

This decay is of course a consequence that the system is open and radiates energy 
away to infinity. Consequently, the corresponding "radiation boundary condition" is 
composed only of outgoing waves which are present far away from the cavity. This is 
the type of boundary condition that occurs in many emission problems, in particular in 
lasing. To see how these "leaky modes" , which we found in the treatment of the passive 
cavity, arise in lasing, we have to take a new point of view: Consider a scalar wave 
equation of the form 

V^ilj + n^k'^ilj = 0, (314) 

as it arises for each of the coupled polarization directions in the cylinder. This is similar 
to ()273|) where we made the ansatz of a steady-state time dependence with real k. If we 
do indeed require k to be real, then in order to find a solution satisfying the radiation 
condition, we permit the refractive index to be complex inside the resonator, which we 
denote by n = n — in', where n is the real part. The imaginary part allows to introduce 
an amplifying medium. Consider the simple example of a plane wave exp[inA;x] which 
will clearly grow with x in the direction of propagation. This is the "energy source" , but 
note that the detailed mechanism of this energy production, for example the pumping 
mechanism, is not defined in this way. Nevertheless, outside the cavity, we again assume 
air with n = 1 with matching conditions that are appropriate for the given polarization. 
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Finally, let us assume TM polarization as before in the scattering problem. The solutions 
for the exterior and interior field are ipext and ipint- 
Let us now recast ()314|) as 

V'^iIj + n'^k'^ ip = , (315) 

where n is the real part of n as defined above, and k = k — ikn'/n is the complex 
wavenumber inside the cavity, which reduces to k = k outside. This recast equation is 
almost of the form of ()288|) . except that k has different values inside and outside the 
cavity. If instead, we also had k = k — ikn'/n outside, the solutions of (j315p were 
exactly the quasibound states of the passive resonator (i.e. solutions of the Helmholtz 
equation at complex k) if the deacy rate was defined as as 

7 = ckn' /n. (316) 

We now may argue that adding or dropping the imaginary part of k outside the resonator 
makes only a small difference if 7 is small: In fact, as can be checked (cf. ^J) by 
inspecting the large-argument asymptotics of the Hankel function if^^\ the field of a 
quasibound state at distances larger than ^ 0/(27) from the cavity grows exponentially 
due to retardation, but within this physical range ipext vanishes as 7 — >■ 0. Therefore 
one can write tl^extix) ~ 7C(r)- If we expand 'ipint{,l) and ^(7) in a Taylor series in 7, 
then to linear order the 7-dependence of C, but not that of ipint, can be dropped in the 
full solution. Therefore, the stationary state of the active medium described by fl314|) ) 
or ()315|) . as well as the metastable decaying state obtained by replacing the real valued 
outside wavenumber k hj k — ik n' /n are identical to first order in 7 within an area of 
order 7"^. 

In conclusion, we have seen that long-lived quasibound states appear not only as 
resonances in elastic light scattering but also as stationary states in emission problems 
at low pump rate (characterized by 7 = c/c n' /n). When asking for the intrinsic properties 
of such long-lived states, it is often a good starting point to consider first the limiting case 
of no losses, i.e. the closed cavity. In the next section, we discuss the familiar concept 
of Gaussian beams as freely propagating waves, and how they can be used to describe 
resonator modes. Although since the work of Fox and Lee jHl] this is the standard 
approach to describe laser resonators 28], this approach has some severe limitations for 
microresonators. 

6.5 Paraxial approximation and the parabolic equation 

In free space, the spectrum of electromagnetic waves is not only continuous but highly de- 
generate, as can be seen, for example, by choosing plane- wave solutions of the Helmholtz 
equation and noting that at each frequency an infinite number of propagation direc- 
tions can be chosen. There exists a class of approximate free space solutions of the 
Helmholtz equation, the Gaussian beams. Gaussian beams are obtained by introduc- 
ing the slowly-varying envelope approximation in the wave equation. This is a special 
case of a frequently encountered approximation which is labelled with different names. 
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for example adiabatic approximation, paraxial approximation, or Born-Oppenheimer 
approximation. In the theory of partial differential equations one also finds the name 
"parabolic-equation" method, because the main step is to neglect the second derivative 
of a slowly-varying "envelope" function, which results in an equation of the type of the 
time-dependent Schrodinger equation for the latter. In optics, this approximation leads 
to the Fresnel propagator describing diffraction in the limit of small angles with the 
optical axis. 

6.5.1 Gaussian beams and the short-wavelength limit 

Gaussian beams play a prominent role in traditional laser optics, because they arise 
naturally as solutions in parallel-mirror resonator structures, and they also represent 
diffraction-limited propagation. However, they are not in general acceptable functions 
to describe the modes of a microcavity. In contrast to free space, it is not even clear a 
priori how to define the equivalent of the paraxial condition in a general microresonator 
structure, such as for example a circular disk. To clarify this point, we will discuss the 
paraxial approximation as next. 

Behind the paraxial approximation is the picture of a "light beam" propagating in 
2;-direction. Let us assume that the beam is described by the envelope function 

ilj{x,y,z)=u{x,y;z)e''' (317) 

which we insert into the Helmholtz equation 

V^V^ + k'^i; = . (318) 

As aresult we obtain 

Uxx + Uyy + Uzz + 2iku^ = . (319) 

With this step we have formally reduced the wavenumber exponent by one, without loss 
of information. 

In the paraxial approximation we take ()319|) . but neglect the second order derivative 
Uzz = d'^u/dz^. In what situation is this appropriate? The name "paraxial" comes from 
a comparison between length scales over which wave solution extends along z and in the 
transverse directions, respectively. Let L be the relevant length in the z direction. Now 
let us introduce rescaled coordinates which measure the transverse space in dimensionless 
form as jSH] 

x' = x ^fk/L, y' = y y/k/L, , (320) 

while for the longitudinal dimension we set 

z = z/L . (321) 

With this ()319|) appears as 
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Ux'x' + "Uy'y') H n Uz'z' + 2i/i; — Uz' = . (322) 
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With multiplication by L/k we obtain 

1 
kL 



Ux'x' + Uy'y' + -j-y Uz'z' + 2iUz' = . (323) 



Thus, the second order z' derivative can be neglected, if 1/kL <^ 1. In cases in which 
this is fulfilled the paraxial approximation is appropriate. 

As 1/kL = \/2ttL -C 1 implies A <^ L, the paraxial approximation corresponds in 
fact to a short-wavelength approximation. The approximation consists in the assumption 
that the envelope u is slowly-varying compared to the wave fronts along z. The resulting 
equation in the original coordinates reads as 

^' + 1^ V = 2iA:|^ . (324) 



Qj.2 Qy2 J Q^ 

This equation is mathematically of a fundamentally different type than the Helmholtz 
equation: whereas the former is an elliptic partial differential equation, the latter is 
parabolic. The difference is that solutions of an elliptic equation are determined by 
solving a boundary-value problem, while solutions of ()324j) are specified by initial con- 
ditions along a line z = zq, for example. If we try to impose boundary conditions on 
()324p by specifying the values of m on a closed surface in x, y, z (e.g. a microresonator), 
there will in general be no solution for any k, whereas the original eigenvalue problem 
(j318|) always exhibits solutions satisfying the boundary conditions for a discrete set of 
/c-values. This well-known consequence of the theory of characteristcs jHI] by which 
second-order partial differential equations are classified explains why paraxial approxi- 
mations in general do not allow to find all the modes of a resonator. 

Let us put it in other words: Equation ()324|) is formally analogous to the time- 
dependent Schrodinger equation in which we substitute z ^^ t and k ^^ m/h. We know 
that this equation uniquely determines the time evolution of any initial wave packet 
prepared at time t = to- This corresponds to specifying the solution of (|324p at some 
z = zq. Clearly, if in addition to the initial condition we also want to impose an arbitrary 
"final distribution" of u at some other time t = ti, then, in general, a contradiction will 
arise to the unique time evolution. Only for special choices of initial distributions can 
the final condition be met. Boundary conditions on a closed surface in the "spacetime" 
spanned by x, y, z are therefore not generally consistent with (j324j) . Now, based on this 
parabolic equation how can we find the modes of a closed resonator? 

6.6 Gaussian beams in free space 

There are various schemes by which this can be accomplished for some subset of modes, 
and the vast majority of conventional resonator calculations in optics are based on such 
methods [2H1- In order to discuss these approaches and their limitations for micro cavities, 
we first recall how the propagation along z (or the "time evolution") can be expressed 
in terms of the retarded Green function G: 

u{x',y';z') = G{x',y',z'] x,y,zo)u{x,y;zo)dxdy . (325) 
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The Green function of the free particle time dependent Schrodinger equation, also called 
its propagator [^0], is well known from quantum mechanics [^ as well as from the 
theory of thermal conduction [61j (where it appears with an imaginary time scale): 

r h 1 '^/^ 
G{f',t';rM)= ^ ..I . e^'^^^^'-^^'/^i^'-^) {k^m/h). (326) 

Zn 1 [t — ij J 

Here, d = 1, 2, 3 is the number of spatial dimensions. This result reminds us of an- 
other important feature of the parabolic equation: The above function describes how 
a wavepacket localized at a single point r at time t spreads with time. Such pulse 
spreading is completely absent in the time dependent Green function of Maxwell's wave 
equation, (V^ — d'^/d{ctY)u = 0, where the time dependence is simply given by the 
effect of retardation, G oc 5{t' — t — ^Jx^ + y^ + z^/c) (in three dimensions). 

The discussion in quantum mechanics textbooks usually applies to three space dimen- 
sions (c? = 3). However, we require the analogous objects for the case of two spatial 
dimensions (d = 2) - in the paraxial approximation, z takes the role of a time param- 
eter and hence leaves only x and y as spatial coordinates. The reduced dimensionality 
introduces different prefactors in Green's function but leaves the basic implications of 
the paraxial approximation (including pulse spreading) unchanged. With d = 2, ()325|) 
describes the transverse profile of a paraxial beam, as it propagates along the z axis. 
The equation is identical to the Fresnel wave propagation formula [HH] for the field ip, if 
we use the definition of u in ()317|1 . 

The Green function in 13 2 5^ alias the Fresnel propagation kernel, defines an integral 
equation for the transverse mode functions u{x,y), which is a frequent starting point 
for determining resonator modes in the paraxial approximation. Before we turn to the 
problem of resonator modes, it is useful to consider first the eigenf unctions of ()325|) in 
free space which are called Gaussian beams. The derivation of the fundamental Gaussian 
beam profile follows the analogy between (j324j) and the time dependent Shrodinger 
equation. An important solution is the minumum-uncertainty wavepacket, 

(327) 






X 



Aa^ + 2iz/k 



where a is the width of the wavepacket in x direction at "time" z = 0. An analogous form 
can be found for the other transverse direction, y, and the full wave function requires 
taking the product of both profiles. This separability is a consequence of the fact that we 
have imposed no boundary conditions on the free-space propagation. We can therefore 
restrict ourselves without loss of generality to wave functions that depend only on one 
transverse coordinate (x) for now. 

Note that ()327p is a special minimum uncertainty wavepacket with zero average x- 
momentum. By the uncertainty relation, a is related to to the spread in wavenumber 
k;j. through 

aAA;^ = l/2. (328) 

The wavepacket therefore has a nonzero width Akx in Fourier space, which will lead to 
divergence as a function of "time" z. By decomposing prefactor and exponent into their 
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real and imaginary parts, it is straightforward to show that ()H27|1 describes a standard 
Gaussian beam [HI], if we reinstate the definition ip = u exp{ikz) and interpret 2 a as 
the minumum spot size at z = 0. The divergence, or angular beam spread, is given by 

e = arctan ( 7— ) ■ (329) 

It must be emphasized that the Gaussian beam is not an exact solution of the free-space 
Helmholtz equation ()318p . but becomes exact only in the paraxial limit. The beam 
spreading is, however, not an artefact of this parabolic equation method. It correctly 
describes the conjugate relation between spot size and transverse wavenumber which is 
a consequence of diffraction. 

6.7 Gauss-Hermite beams 

The spot size 2 cr is a parameter at our disposal, by which we can determine the relative 
uncertainties in position and momentum according to ()328|) . This freedom is important 
when one attempts to construct approximate resonator modes by properly adjusting 
Gaussian beams, as we will see below. Further generalizations of the fundamental Gaus- 
sian beam are possible in free space, by allowing the transverse amplitude to vary. An 
important special case are the Gauss-Hermite beams which are constructed as follows. 
If we set 2; = in ()327|) . then the resulting instantaneous Gaussian 

can be formally interpreted as the ground state wavefunction of a one dimensional har- 
monic oscillator, 

in which formally identify the oscillator length with the beam spot size, 



h 

V2cx. (332) 



mu 



Higher order eigenfunctions can be generated by acting on uq with the formal creation 



operator^ 



at = i=(Q-iP), (333) 



where we defined the operators 



1 dx 1 ax 



^Note that the operators we define now are motivated solely the mathematical structure of the problem. 
They do not represent quantum mechanically relevant entities like position or momentum, nor are 
they physically related to the operators of section |5] 
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The solutions are products of Hermite polynomials and the Gaussian uq. From this 
construction at z = 0, however, it is not directly obvious how the beam shape changes 
as it propagates along z. To see this we will define a generalized creation operator 
c^ which properly describes the spreading (divergence) of the beam, and which, when 
applied to the z-dependent u{x,z) of (|327p . generates functions which are solutions of 
the parabolic equation (j324|) . 

For propagation in the x — z plane to which we restricted our attention, the transverse 
beam coordinate is x, whereas z is the propagation direction playing the role of a time 
variable. Let us define the following dimensionless quantities 

q(kz) = pkz — 2ip{ka)'^ , 
where k and a are the parameters from ()327|) for which we thus obtain 

{kxf 



^"<--)^^\/^- 



ip- 



2q{kz) 
We can now define a generalized form of the creation operator ()333j) as 



(337) 



which in the limit z = reduces to ()333j) because then q = —ika. One can verify 
that ()338p and ()339|) lead to the correct commutation relation [c, c^] = 1, as a direct 
consequence of the fact that the quantities defined in ()336p satisfy the relation 



q*p-qp* = i (340) 

for all z. We included complex conjugation of p even though it is real, in order to 
highlight the similarity to the quantum commutation relation between conjugate oper- 
ators - and indeed, p can be interpreted formally as the conjugate momentum of the 
complex trajectory q{kz) because it satisfies the equation of motion p = dq/d{kz) with 
dimensionless time kz. 

For c^ generating Gauss-Hermite beams, this operator must exhibit two additional 
qualities which we check in the following. 



First, we verify by direct application of ()339|) that c annihilates the fundamental 
beam u{x,z) given in ()327p . 

Second, we verify that {c'')"^u{x, z) with m = 1,2. . . is again a solution of the 
free-space parabolic equation 
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This latter property requires that c^ commutes with the differential operator 

which in fact corresponds to ()H24j) written in operator form. Thus, in operator notation 
the parabolic equation simply reads as 

lu = 0, (342) 

where we still assume that the y dependence can be separated off. Using the definitions 
in Eqs. fj338p and ()336p . we obtain 

''''-''' = ' iw^) ah -"'ak))' ''"'' 

where the righthand side vanishes because of the equation of motion dg*/d(A;2;) = p*. 
This result means that if m is a solution with L m = 0, then c^u is also a solution, because 

l{c^u) = c\lu)=0. (344) 

Thus, we have shown that the set of Gauss- Hermite beams are solutions of the free space 
wave equation in paraxial approximation. 

6.8 Resonator modes in the parabolic equation approximation 

The construction of resonator modes using the results obtained so far requires us first to 
identify what is meant by the z axis along which the beams propagate, if the system is 
bounded by some surface. For simplicity, we first assume that this boundary is perfectly 
reflecting, what corresponds to Dirichlet boundary conditions. We can approach this 
situation from two logically distinct directions. 

Traditionally it is assumed that the mirrors are shaped such that they precisely match 
the wavefronts of a Gaussian beam for some value of spot size a and wavenumber k. 
On the other logical side one can start with a desired field configuration for which the 
shape of suitable mirrors is sought. The goal is to find modes with a predefined form 
which is optimally adapted to a given practical problem. In contrast to the traditional 
approach this procedure has the form of an inverse boundary-value problem. It has been 
addressed in Ref. [HH], and plays a role whenever beams have to be matched to a given 
situation. 

The molecular sieve crystals on which this review focuses, as well as all other self- 
assembled realizations of dielectric microresonators such as aerosol or polymer droplets, 
pose instead a well defined boundary value problem with little or no freedom of adjusting 
the cavity shape. This forces us to take a different route if we wish to extend the results 
of the previous subsections. 

The application of the parabolic equation method in the diffraction theory of res- 
onators is presented in great detail in Ref. jSH], where, however, the salient features are 



80 



at times obscured by the mathematical apparatus that is necessary to treat wave prop- 
agation in three dimensions with inhomogenous refractive index. In the folloeing we do 
not take such comphcations into account, but approach the solution of the Helmholtz 
equation ()318|) with only x and z as variables, and Dirichlet boundary conditions on the 
cavity wall which in view of the situation with microresonators, we assume to be a two 
dimensional closed surface. The idea of this approach is that since the Gaussian beam 
is an approximate short wavelength solution near some z axis, we can piece together 
a resonator mode by using Gaussian beams which follow a dosed ray trajectory in the 
cavity. An example for such a closed orbit is shown in Fig. ^1 It is a nontrivial example 
where a ray based mode analysis of the present type leads to accurate results jHZI. Each 






Figure 13: (a) A stable periodic orbit in an oval cavity. The path can be parametrized 
by the arc length s, i.e. the longitudinal distance along the ray from the 
(arbitrary) starting point A. reflections occur at lengths Sk- The stability of 
this bow tie shaped pattern is apparent in (b) where a ray is launched with 
initial position and direction slightly off the periodic path: the trajectory is 
nevertheless confined to the vicinity of the bow tie, performing an oscillation 
transverse to the trajectory in (a). The transverse deviation in position can be 
measured by a coordinate q as shown in (c), and likewise p denotes a deviation 
in momentum (or direction) of the initial ray. The outcome of the reflection 
for such a perturbed ray is shown dashed. Deviations are exaggerated for 
clarity. 



straight segment of the trajectory between two successive reflections defines a z axis for 
which the analysis of the previous subsections can be carried out. Assuming a solution 
of the form 13371 in each segment, we determine p and g in a generalization of (j336|) 
such that two requirements are satisfied: The boundary conditions must be met to the 
accuracy of the paraxial approximation, and the full solution must be single valued (i.e. 
reproduce itself) after a complete round trip around the ray orbit (constructive inter- 
ference condition). We retain the form given in ()336|) . with one modification. Without 
modification, the minimum width of the Gaussian beam is located always at z = 0. We 
can shift this location by allowing p to be complex, and still maintaining the relation 
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q*p—pq* = i underlying the construction of the Gauss- Hermite modes. Thus we define 

q{kz) = Aikz + Aq , 

where Aq^i are complex but restricted to Rs^AqAi) = 0. 

To follow the wave around the closed trajectory, we introduce a new longitudinal 
coordinate s which measures the distance along the whole ray loop from some arbitrary 
starting point, say point A in Fig. E| Then the individual reflections occur after the 
lengths Sk, k = 1, . . . N, where A^ is the number of bounces. Furthermore, let L be 
the length of the close path. In order for the wave function to be single- valued, the 
parameter q governing the spreading in the transverse wave function must be periodic 
as a function of s with period L, for example 

q{s + L)= q{s) (346) 

for all s in the interval of straight propagation we consider. 

6.9 Monodromy matrix 

From the boundary condition of vanishing wavefunction at the point of reflection, we 
get additional relationships for p and q. Firstly, the phase of the Gaussian beams before 
and after reflection must exhibit a jump by tt. Moreover, considering a reflection at the 
point s = Sk, one flnds that the vector {q,p) before and after reflection (at s = Sk — As' 
and s = Sk + As, respectively) must obey the equation 

where, M^'^^As, As') is a 2 x 2 matrix which depends on the distances As, As' from the 
reflection point. It is called the monodromy matrix of the given reflection k. Its explicit 
form can be stated in terms of the local radius of curvature R = R^k) of the reflecting 
surface and the angle of incidence X = X^k) with respect to the local normal [02] as, 

/ 2As _ 1 (As)(AsO \ 
M('=)(As,As')=M™'^ 2AS™''' ' (348) 



where 



R cos X R cos X 



, 2 cosx cosx\ .„.„^ 



The linear matrix equation ()H47|1 holds only within the paraxial approximation, which 
for the reflections means that the transverse spreading at the reflection points shall be 
small, as measured by the real part of the exponent of ()337p . which is 

^ikx)Hm^ = hkx)^\q\-\ (350) 
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where we used q*p — qp* = i. A small extent in x relies therefore on \q\/k being small. 
Away from the beam waist |g| grows with z, and the growth is smaller when k is large. 
Hence, the required limit is again that of large k. For details of the calculations, the 
reader is referred to Ref. 



6.10 Round trip stability of periodic orbits 

As is expected from a short-wavelength aproximation, the monodromy matrix itself has 
a ray interpretation. Equation ()347|) is precisely the transformation that determines how 
a small deviation from the periodic ray path changes the outcome of the reflection. As 
noted earlier, q and p behave like conjugate variables, and here this fact emerges again. 
If we interpret formally g as a cartesian coordinate transverse to the ray direction along 
which s varies, andp as the conjugate transverse momentum, then the true periodic orbit 
would be described hj q = p = and s = . . . L. Nonvanishing p and q at some point 
s = Sk — As' denote transverse deviations from the periodic orbit. This leads to straight 
line trajectories hitting the surface slightly off the locations where the exactly periodic 
orbit has its reflections. The outcome of the reflection of this displaced ray is then a ray 
deviating from the periodic path by g(sfc + As) in position and p(sfc + As) in momentum. 
To show that these classical quantities are related by ()347p in the linear way of a Taylor 
expansion for small deviations requires straightforward but tedious trigonometry [JHl ISHI • 

The connection between the two requirements of (|346|) and (j347p becomes apparent, 
if we extend the monodromy matrix from a single reflection to include all N reflections 
of a full round trip, and ask for the cumulated deviation of the ray, as a function of the 
initial displacement in position q and direction (momentum) p. As a result we can see 
that in a linear approximation we simply can replace M*^'^^(As, As') in ()347|) with the 
product M[s) of the N individual monodromy matrices. Assume that we start a ray at 
longitudinal coordinate s in the first leg of the ray path, Si < s < S2, with deviation q 
and p from the periodic orbit [cf. Fig. ITST a)]. We are then interested in the deviation 
incurred along one round trip s + L. 

The distance from s to the first reflection is As' = si — s, and the distance between the 
first and second reflection accordingly is As = S2 — Si. If we denote the length between 
reflections at s^ and s^+i by l^, then the round-trip deviation is obtained from 

where for the product M{s) of the A^ individual monodromy matrices is given by 

M{s) = M^''\s - sn,In-i) X M^''-'\Im-i,In-2) 
X ...xM(2)(/2,/i) xM(i)(/i,si-s) . 

The periodicity condition ()346|) requires that the vector d = {q{s), p{s)) is an eigenvector 
of matrix M{s). Up to now we retained the dependence on the starting coordinate s. 
However, the investigation of the eigenvalue problem becomes simpler when we eliminate 
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this dependence. This is possible because we know how any given trajectory depends on 
s: Rays in the cavity follow straight lines, so that g is a linear function of s, whereas p is 
in fact a constant. (This can be compared to the similar assumption we made in ()345|) 
with kz instead of s as the longitudinal coordinate). For the vector d we can therefore 
write 

:)- GO (a) -"«(':)• '^=^) 

where the matrix n(s) = ( J f ) appearing here obviously has unit determinant, indepen- 
dently of s. Thus we can rewrite the periodicity condition ()35H) as 

(;^;) = n-^(. + L) M{sMs) (;^°) , (354) 

which is an eigenvalue equation. We know that the coefficient vector ( yi° ) ^^^^ is 
by construction s-independent, and as a consequence, the s dependence of the matrix 
product E = n~^(s + L) M{s)Il{s) must cancel out. The matrix E is the monodromy 
matrix of one whole periodic orbit, because it characterizes the propagation of small 
deviations along one round trip irrespective of where exactly we started. Recall, however, 
that we set the starting point s on a specific branch of the periodic orbit. Clearly, the 
same arguments would apply with s starting in any other branch, but then the resulting 
matrix equation will in general be different. 

6.11 Eigenvalues of the monodromy matrix 

Since ii^ is a 2 x 2 matrix, its eigenvalues and eigenvectors are easily classified, especially 
since E can also be shown to have unit determinant 1^2] . This follows from det(n) = 1 
together with classical area preservation (Liouville's theorem) for the ray dynamics; the 
latter implies that also det(M(s)) = 1. As det E = 1 the two eigenvalues Ai,2 satisfy 

A1A2 = 1 and Ai + A2 = Tr(^) = t (355) 

so that determinant and trace of E fix its eigenvalues via the quadratic equation 



Ai,2 = -t±-Vt2-4. (356) 

Three cases can be distinguished depending on the value of t: 
t > 2: The eigenvalues are real and have different magnitude, 
t = 2: Both eigenvalues are degenerate, Ai = A2 = 1. 
t < 2: The eigenvalues are complex and both have unit magnitude, 

X^,2 = lt±^VA~^ = e^"^. (357) 

In this case the eigenvectors hi 2 are also complex conjugates of each other. 
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The periodic ray orbit under consideration is called linearly stable if and only if case 
t < 2 holds. If this is satisfied in one branch of the path, one can show that it also holds 
in all others. 

The difference between these cases can be appreciated if one writes a general vector 
d = {q, p) as a linear combination of the eigenvectors, 

d = ahi + hh2 (358) 

and asks for its evolution under repeated application of the round trip mapping E. For 
a v-iold round trip we obtain 

E^rf = aE^hi + bE''h2 = a\\hi + 6A^/^2 = ae''"''hi + he-''"^h2- (359) 

If c? is a real vector - as it should be if it describes real ray trajectories - then E^d is 
also real because ii^ is a real matrix, describing only the trigonometry of the successive 
reflections and the straight line ray motion in between (cf. the definitions of 11 and 
M[s)). Therefore, we know that in the last equation there can be no imaginary part, 
and thus nothing is changed by taking the real part. If hi = [Ao, Ai), then we have 

™ r^ fRe[{a + b*)Ao] cos z/0 - Im[(a + b*)Ao] sin u(f)\ , . 

\Re[{a + b*)Ai] cos z/0 - Im[(a + b*)Ai] sin uc^J ' ^ ' 

Remember that is a fixed eigenphase, characteristic of the given branch of the pe- 
riodic orbit. If we assume for the moment that u is not an integer but instead varies 
continuously, then the above equation describes an ellipse in the plane spanned by q 
and p. For integer z/, the trajectory simply visits discrete points on this elliptical curve 
every time it completes a loop. This means that the small initial deviation d from the 
periodic ray orbit leads to a trajectory which stays close to this original periodic path 
for all times. The perturbed ray only performs a small oscillation in transverse position 
and momentum around this original path. 

In case t > 2, the curve described by E'^d is also a conic section, but in this case a 
hyperbola instead of an ellipse. To see this, we go back to (j36(jp but with real eigenvalues 
of the form Ai = I/A2, cf. (I355p . Since for real eigenvalues the eigenvectors of the real 
matrix E can also be chosen real, we obtain with ^ = A'' 



i.f 



E''d = ^ahi + -bh2, (361) 

which describes a hyperbola if ^ is formally allowed to vary continuously. This con- 
struction tells us that the perturbed ray deviates more and more from the periodic path 
with increasing v. Hence a periodic orbit whose monodromy matrix has trace t > 2 is 
unstable, and for t < 2 it is stable. 

The question of stability is crucial for the applicability of the paraxial approxima- 
tion in constructing resonator modes. Only if a ray orbit is stable can the paraxial 
approximation be justified over the whole round trip. In this case one can say that the 
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mirror configuration as seen by the closed ray trajectory and determined by the radii 
of curvature R at the successive reflection points acts in a focussing way. Recall that 
the spreading of the Gaussian beam in each ray segment is given by ()350|) . This will 
be unchanged upon a round-trip if q acquires only a phase, not a change in absolute 
value. If we determine q according to (j353j) with h = {Aq, Ai) being an eigenvector of the 
monodromy matrix E, then q{s) indeed changes only by the eigenvalue e^'^ when after a 
round trip we return to the starting point s. This is a mere phase change provided that 
the orbit is stable. For unstable orbits, the eigenvalues are real and different from unity, 
making it impossible to recover a beam with the original spread after one period of the 
orbit. 

6.12 Resonator eigenfrequencies in the parabolic approximation 

The use of the parabolic (Gaussian beam) approximation requires that we identify the 
stable periodic ray orbits of the resonator at hand. This can in general be a difficult 
task in itself, because their number and location can depend rather sensitively on the 
shape of the boundary. Furthermore, it must be kept in mind that there may also be 
modes which are not in any way associated with stable ray orbits at all. 

Assuming that we have found a set of stable orbits, we can then use the Gaussian 
beam solutions to construct the corresponding eigenfrequencies (within the parabolic 
approximation). This is done by recalling that the longitudinal coordinate s along the 
orbit plays the role of a time-like variable due to the approximation made in (|324j) . In 
order for the modal wave function to be single- valued upon a round-trip, the dependence 
on s has to be periodic. The imposition of periodic boundary conditions introduces a 
longitudinal "quantum number" m counting the nodes along the ray path. The equations 
determining the wave solutions are themselves periodic in s, and hence the problem is 
analogous to that of a periodically driven time dependent system to which Floquet's 
theorem applies [SHj, or equivalently to Bloch's theorem governing the band structure of 
spatially periodic solids. The central observation there is that the wavefunction acquires 
only a phase factor when a translation amounting to the periodicity interval is carried 
out. This is precisely the content of (|357|) which however holds only when a stable orbit 
is present. 

Allowing for transverse excitations according to the harmonic oscillator ladder as 
described by the generalized operators ()338|) . an additional transverse mode number n 
labels the number of nodes in the two-dimensional Gauss- Hermite beam perpendicular 
to the propagation axis. The resulting quantized wavenumber is found to be 



kL = 2'K (m + —\ + (n + ]A 4). (362) 

where L is the round-trip path length of the stable orbit. A^ is the number of reflections 
the ray undergoes, and this number enters here to take into account the phase shifts 
associated with reflection. The phase multiplying the transverse excitation state is the 
eigenphase of the monodromy matrix as defined in (|357|) . A simple application of this 
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formula would consist in determining the longitudinal mode spacing. The intuitively 
expected result is that an integer number of wavelengths must fit around the ray path. 

6.13 Polygonal resonators 

6.13.1 Marginally stable orbits 

The stability considerations of section IB. 1 1! provide a conceptual background helping us 
to understand why pathological effects can occur when the resonator is bounded by a 
polygon with straight sides. The example most relevant here are the hexagonal facets 
of a molecular sieve crystal we will present in section [3 In case of the hexagon, the 
boundary allows neither stable nor unstable orbits. We have precisely the marginal case 
t = 2, where the round-trip monodromy matrix exhibits Tt{E) = 2, so that both of its 
eigenvalues are unity. The sides of the cavity are then neither focussing nor defocussing, 
and the rays deviating slightly from a closed ray path do not describe an elliptical 
transverse envelope as they would for a stable ray. One obtains this result by evaluating 
the monodromy matrix for a plane interface. The resonator modes for a straight sided 
cavity are thus not well approximated by Gaussian beams, as can also be seen from 
the simple example of a rectangular cavity, for which the solutions of the scalar wave 
equation are (for propagation in the xy plane) 

ipix, y) oc sin /c^(x - a;o) sin ky{y - yo) (363) 

with suitably chosen discrete wavenumbers kx^y and offsets xq, yo. 

This simple example bears no resemblance to Gaussian modes but can be solved ex- 
actly because the problem is separable in cartesian coordinates. However, not all polyg- 
onal resonator geometries permit a separation of variables. This applies in particular 
to hexagonal resonators - the geometry exhibited by the facets of the molecular sieves 
microcrystals with which microlasers were realized; cf. [7| For a polygon with precisely 
120° angles between adjacent sides, any ray lauched at some angle to the surface will go 
through only a finite number of different orientations [02]! just as in the more familiar 
rectangular resonator where there are at most two non-parallel orientations for any ray 
path. In the hexagon, a ray encounters the interface with at most three different angles 
of incidence. Despite this apparent simplicity, there exists no orthogonal coordinate sys- 
tem in which the wave equation for the hexagonal cavity can be solved by separation of 
variables. This property is not surprising, since there is only a finite number of suitable 
coordinate systems available, and on the other hand there are infinitely many shapes 
one can think of. 

6.13.2 Triangular and hexagonal resonators 

Even when separability is lacking, the resonator problem may sometimes be solvable by 
an appropriate ansatz. An example for this is a resonator composed of three straight 
sides forming an equilateral triangle on which the fields are assumed to vanish (Dirichlet 
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boundary conditions). The wave solutions in cartesian coordinates are given by 
V'm,n(a;, y) = sin ( —(2m -n)x] sinl—nyj- 

sin I — [2n — m) x \ sin I —= my \ + (364) 

• /^-27r, , \ . /27r, , 

sm I [m + n) X \ sm I —=[m — n)y 

This solution is slightly more complicated than the one for the rectangular resonator 
of (j363p . In (j364|) there is a degenerate set of solutions V'm.n of different symmetry 
type which are obtained by replacing sin with cos in the x-dependent factors. We get 
identical wavefunctions if we exchange the quantum numbers m and n, replace n by 
m — n,OT replace m, n by —m, —n (simultaneous sign change); the solutions differ only 
in sign so that they do not represent different eigenfunctions. There are, however, true 
degeneracies as well. They are characterized by combinations of m and n giving the 
same wavenumber 

km,n = — Vm? + n^ — mn . (365) 

o 

Since any linear combination of eigenfunctions with the same k is also an eigenfunc- 
tion, there is a considerable ambiguity in the spatial distribution of the wave function 
whenever there are degeneracies. For example, for any given m, n we can form the 
superposition 

^m,n = ^m.n + i ^m,n • (366) 

Shown in Fig. HH are the three functions V'm,n, V^m.n and '^m,n for the particular case 
m = 14, n = 21. 



Figure 14: In the equilateral triangle with Dirichlet boundary conditions, degenerate 
wave functions (a) and (b) with m = 14, n = 21 are superimposed according 
to ()366j) to obtain a completely different spatial intensity pattern (c). Plotted 
is the absolute square of the wave function in a grayscale representation with 
black as the highest intensity. 
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The wavefunction shown in Fig. El (c) has nodal hnes along the edges of a set of smaller 
triangles inscribed into the original one. The self similarity exhibited by this solution 
bears an intriguing resemblance to the intensity patterns recently reported in unstable 
laser resontors with a mirror in the shape of an equilateral triangle [70] . Although further 
investigations of this emerging field of open, unstable resonators are needed, it seems 
possible that the transverse mode structure observed in the above laser experiment is 
given by a two dimensional wave equation with solutions similar to the ones shown in Fig. 
IT^ The particular linear combinations appearing in any given experimental realization 
can be a result of the measurement setup or perturbations that lift the degeneracies. 
In particular, complex-valued superpositions as given by ()H(J(ij) become relevant in open 
resonators such as those in Ref. jZO]- Whereas for a closed resonator with Dirichlet 
boundary conditions one can always choose the solutions to be real without losing any 
modes, this is not generally possible in open systems, when travelling waves are assumed 
to describe the field outside the resonator. 

The subdivision into smaller triangles in Fig. IHT c) suggests that modes of other 
polygonal resonators can be constructed in an analogous way if the boundary can be 
decomposed into eqal equilateral triangles. This is the case in particular for the hexagon, 
which is composed of six equilateral triangles. A hexagon mode can hence be pieced 
together from six copies of any given solution of the triangle problem. However, solutions 
obtained in this way do not represent all the possible modes of the larger hexagon 
resonator. This is because the triangle solutions necessarily satisfy Dirichlet boundary 
conditions on all the sides - also on those which form the diagonals in the hexagon 
composition. Clearly, we can expect that there should also be solutions with even parity 
under reflection at such a diagonal. These would then have vanishing derivative on the 
respective diagonals, i.e., satisfy Neumann boundary conditions. Unfortunately, there is 
no exact ansatz such as (j364p for these different symmetry classes of the hexagon. 

The hexagon with Dirichlet boundary conditions on its faces therefore is an example 
of how not only separation of variables breaks down, but also the more general concept 
of integrability fails. The latter essentially means that one can find as many (globally 
valid) "good quantum numbers" as there are confined degrees of freedom; we managed 
to do this in the equilateral triangle but cannot label the solutions of the hexagon 
using only these same numbers. Integrability carries over to the classical ray dynamics 
in the form of constants of the motion. Nonintegrability is known to be one of the 
central characteristics exhibited by wave equations whose classical (short wavelength) 
limit exhibits chaos [00] • When the ray dynamics is chaotic, this does not imply that 
the rays move in a stochasitc way. They are still governed by the deterministic laws of 
total internal reflection and refraction, and in particular infinitely many periodic orbits 
can be found. However, if for the extreme case of a fully chaotic system we calculate the 
monodromy matrices for these periodic orbits, they all turn out to yield TrE > 2 and 
hence describe unstable (hyperbolic) orbits. 

Hexagonal resonators are therefore an example of a "pathological" intermediate class 
of microresonators for which the ray dynamics is not chaotic, but not integrable as well. 
These systems have been termed pseudo-integrable [68j. It is interesting to observe that 
the modes with which the pathological nature of the problem manifests itself are precisely 
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those which are not required by symmetry to vanish on the diagonals of the hexagon. Or 
with other words, parity with respect to the diagonal does not in itself predetermine the 
value of the wave function at the corners of the hexagon. Thus, the corners are the root 
of pseudointegrabilty in these polygons with Dirichlet boundary conditions, because they 
give rise to corner diffraction. Corners are singularities at which the tangent direction 
changes discontinuously. Diffraction arises here because the length scale over which the 
tangent direction changes (the radius of curvature) is ideally zero, and hence certainly 
much smaller than the wavelength |7lJ. The angle subtended by the corners determines 
whether or not a pseudointegrable cavity is created. In the rectangular cavity, as well as 
in the equilateral triangle, the effects of diffraction at the corners cancel each other out. 

An interesting wavelength effect was reported in Ref . ^^ . Sharp corners of a polyg- 
onal billiard are smoothed over the distance of one wavelength. This means that the 
system becomes indistinguishable from a slightly "streamlined" billiard that would be 
obtained by rounding the corners in such a way as to obtain an everywhere smooth 
boundary. This leads to an apparent contradiction: classically chaotic billiards with a 
continuously varying tangent ( "smooth" walls) can be infinitesimally well approximated 
by polygonal domains that are themselves never chaotic, but at most pseudointegrable. 
From a stability analysis of the periodic orbits in either the original chaotic system or the 
inscribed polygonal approximant, one finds that the former exhibits unstable periodic 
orbits, whereas the latter displays only marginally stable periodic paths. Although one 
should thus expect the resulting mode structure to be qualitatively different, the spectral 
structure in a fixed frequency interval can become indistinguishable for the two cases. 
In section 17.21 we will discuss numerical solutions of the Helmholtz equation, and how 
this smoothing on a wavelength scale manifests itself in specific realisations of molecular 
sieve microresonators. 

Because we are interested in dielectric materials forming open resonators, we shall not 
go into further detail concerning the intriguing problem of semiclassically quantizing 
pseudointegrable cavities, but instead we discuss briefly the relevance of corners in the 
presence of leaky boundary conditions. The breakdown of the short wavelength approx- 
imation near corners will also manifest itself in the emission from hexagonal dielectric 
cavities. The concept of pseudointegrability loses its significance in the case of an open 
resonator, because dielectric interfaces can destroy integrability as easily as sharp corners 
do. An example is the rectangular cavity made up of a homogenous dielectric surrounded 
by air. The dielectric constant can be written as e(x, y) = 1 + em0(|a — x\) 9(|6 — y\) 
where e^ is the permittivity of the medium. Due to the term 1, the Helmholtz equation 
cannot be reduced to two one- dimensional problems for the x and y coordinate sepa- 
rately, as was the case in (j363j) . In the limiting case of large permittivity the offending 
1 can be dropped, which restores separability, because this corresponds precisely to the 
case of a closed resonator. As was observed in Ref. ^2) ^^e boundary conditions of a 
polygonal resonator (determined in this case by the spacial distribution of the refractive 
index) can be of greater significance than diffraction for essential (in that study statis- 
tical) properties of the spectrum. It is thus not clear at present if dielectric resonators 
are suitable to search for signatures of pseudointegrability as it was defined for closed 
resonators. 
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7 Actual realizations of microlasers based on molecular 
sieve-dye compounds 

After the discussion in previous sections we can now venture to give a definition of 
what a microlaser is. A microLaser is a structure in which only a few (in the hmit 
one) mode of the hght field interacts with a collection of atoms or molecules, so that 
spontaneous emission processes are enhanced or inhibited, or in which lasing occurs 
without a visible threshold. It is in principle possible to observe these cavity effects 
in large (10'^ . . . 10^ A) resonators with small mode volume (confocal resonators). But 
because of the small frequency spacing (free spectral range) of the longitudinal resonator 
modes, cavity effects in large resonators can only occur when the linewidth of the atomic 
transition is smaller than the free spectral range of the cavity - a condition which can 
only be met by very dilute atomic systems such as an atom beam, but not by atoms or 
molecules in a condensed state. As the condensed state linewidth of atoms or molecules 
is in the order of 1 nm to 100 nm, cavity effects will only be observable with resonators 
exhibiting a correspondigly large mode frequency spacing, that means resonator sizes 
of a one half to few wavelengths. In the past, wavelength size resonators were realized 
with semiconductors as microdisks [711172] or VCSELs (Vertical Cavity Surface Emitting 
Lasers) [7H]-|Hni, and with organic dyes in planar resonators as Langmuir-Blodgett |8l] 
or liquid films [33, or embedded in polymer spheres [H21^IHlj- This work was discussed 
recently in several reviews |HH]-|Hni so that here we will concentrate on microlasers based 
on organic dye molecules embedded in molecular sieve microcrystals [3] . Before we turn 
to the proper optical and laser properties, we will shortly discuss the synthesis of the 
compounds. 

Table 1: Lattice constants and free pore diameter (p of molecular sieves with linear chan- 
nels lUj. 
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7.1 Molecular sieve crystals as host material for microlasers 

As is reviewed in the previous chapter Nanoporous compound materials for optical ap- 
plications - Material design and properties, molecular sieve materials are characterized 
by a crystallographically defined framework of regularly arranged pores. In Table Q we 
list some sieves with wide channel pores, in which optically effective organic molecules 
can be inserted. Among the listed materials especially the aluminophosphate AIPO4-5 
(molecular mass 1463.4 g/mol) can be synthesized with good optical transparency and 
low internal scattering losses. In addition, the AIPO4-5 channel pores exhibit a diameter 
of 0.73 nm, fit to accommodate a variety of laser active organic dye molecules. The struc- 
ture is shown in Fig. ^J AlP04-5-crystals are crystallized from aqueous or alcoholic 
solutions under hydrothermal conditions, with the addition of an organic structurizing 
agent, called template. The template is a necessary device in the synthesis of molecular 
sieves to direct the chemical reactions towards the desired crystal structure. For the 
synthesis of the AIPO4-5 laser crystals tri-n-propylamine was used as template. The 
preferred pH range for the synthesis is mildly acidic to mildly basic, while the source of 
phosphor is mostly orthophosphoric acid, and the sources of aluminum are pseudoboeh- 
mite or alkoxides [HI]. It was shown that single crystals with nearly perfect morphology 
can be grown using hydrofluoric acid |.92j, and with a specially prepared and aged alu- 
minum hydroxide gel crystal sizes around 1 mm in c-axis direction were obtained P^ . 
In addition, microwave heating proved to increase the crystallization rates by more than 
one order of magnitude |9"21 IM] . 

Pure AIPO4-5 crystals are optically transparent from below 400 nm to above 800 nm 
and exhibit a refractive index of n(500nm) = 1.466. After removing the template 
(usually by heating) they show practically no birefringent properties. X-ray diffraction 
revealed patterns which are consistent with space group P— cc as well as P6cc. The 
latter group, however, corresponds to a polar structure: In fact, the 4 in the formula 
AIPO4-5 is the result of the strict alternation of Al and P in the tetrahedral nodes of 
the framework, which prevents the corner-sharing oxygen tetrahedra to occur with odd 
numbers, and which leads to an alternating stacking of Al and P in the direction of the 
channels (c-axis); cf. Fig. El This alternance is assumed to cause the crystallographic 
polar nature ^Hl of the framework ^Q^. The macroscopic polar nature of AIPO4-5 single 
crystals was proven recently in scanning pyroelectric microscopy investigations [HZI, and 
it was observed that AIPO4-5 crystals are usually twinned. The murky stripes inside 
the pyridine 2- loaded crystals shown in Fig. [THland their slightly bowed side faces could 
well be a result of this kind of twinning. At this moment it is not clear, however, to 
what extent twinning should affect the properties relevant for luminescence and lasing 
discussed here. 

Up to now laser action was obtained in two different compounds, namely in AIPO4-5 
loaded with l-ethyl-4-(4-(p-dimethylaminophenyl)-l,3-butadienyl)-pyridinium perchlo- 
rate (pyridine 2 [HHl), and with a modified rhodamine B dye. In the following we will 
discuss the two materials in more detail. 
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Figure 15: Structure of the hexagonal AIPO4-5 molecular sieve. The corners of the 
polygons are occupied alternately with aluminum and phosphorus, while the 
polygon sides represent oxygen. The hexagonal c-axis is oriented along the 
channel pores. 

ClOi" 



Figure 16: Structure formula of the dye l-ethyl-4-(4-(p-dimethylaminophenyl)-l,3- 
butadienyl)-pyridinium perchlorate (pyridine 2 [HHl); molecular mass 
378.9 g/mol. 



7.1.1 AIP04-5/pyridine2 compound 

These AlP04-5-dye compounds were synthesized ^] following a procedure in which the 
dye l-ethyl-4-(4-(p-dimethylaminophenyl)-l,3-butadienyl)-pyridinium perchlorate (trade 
name pyridine 2 [HHI), cf. Fig.^1 is added to the template or to the aluminum hydroxide 
suspension [1()()[ llOJt ll()2j . After 1 h of hydrothermal synthesis dark red crystals with 
a length of up to 100 /im are obtained. As the linear dye molecules fit snugly into the 
0.73 nm wide channel pores of the nanoporous AIPO4-5 host, they become aligned along 
the crystal c-axis. As a result the compound exhibits strong dichroism, and the emit- 
ted fluorescence light is polarized parallel to the c-axis. This is documented in Fig. ^1 
It is also observed that with the inclusion of pyridine 2 the entire compound acquires 
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Figure 17: Morphology of typical pyridine 2/AIPO4-5 crystals; left: without lasing 
properties, right: with lasing properties. 




Figure 18: Transmission micrographs of the dichroism in dye-loaded AIPO4-5 crystals; 
the rod-shaped crystals contain ca. 0.1 wt-% or 1 pyridine 2 molecule per 
260 unit cells, while the barrel-shaped ones (shown in the inset) enclose rho- 
damine BE50 (ca. 0.5 wt-% or 1 molecule per 75 unit cells). Only the polar- 
ization component parallel to the optical transition moment of the molecules 
is absorbed. In the rod-shaped crystal the pyridine 2 dyes are completely 
aligned, whereas with the rhodamine BE50 dye in the barrel-shaped crystals 
we observe only a weak dependence of the color upon the incident polar- 
ization. Left: incident light horizontally polarized; Right: incident light 
vertically polarized. 

pyroelectric properties and an optical second order susceptibility [3]. Even though the 
size of the dye molecules is well compatible with the channel diameter, the dye content 
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visibly affects the crystal morphology. Regular hexagonal crystals with a rodlike form, 
as e.g. the ones shown in Fig. ITTT left). or Fig. ^[ are obtained when the dye content is 
low, around 0.1 wt% or 1 molecule per 260 unit cells. At higher concentrations the dye 
accumulates in the middle of the crystal, and the disturbances grow: at a dye content 
of ^ 0.2 wt%, crystals with a characteristic fascicular shape resulted; cf. Fig lTTT right). 
Given the small size of the crystals the dye content was determined by chemically dis- 
solving them. With this method, however, it is not possible to accurately determine the 
spatial distribution of the dye. Therefore the content was also evaluated qualitatively 
by comparing the depth of the color. It is only with these fascicled pyridine 2-loaded 
crystals where laser emission was observed. Apparently, in the undisturbed rod-shaped 
crystals the low concentration of dye does not spoil the growth, but it is not sufficient 
to provide the necessary optical gain for compensating all losses either. 

7.1.2 AIP04-5/rhodamine BE50 compound 




M(Cli3)£ 



Figure 19: Structure formula of the new dye rhodamine BE50 (ethanaminium, A^-[6- 
(diethylamino)-9-[2-(iV, iV-dimethyl-3-amino-l-propoxycarbonyl)phenyl]-3H- 
xanthen-3-ylidene]-A^-ethyl-chloride) molecular mass 564 g/mol [103]. 



The second molecular sieve dye compound which exhibited lasing properties was ob- 
tained by including a modified rhodamine B dye in an AIPO4-5 crystal. The size of 
rhodamine molecules exceeds the dimensions of the channel pores in AIPO4-5 crystals. 
As a result the rhodamine B is normally not accepted inside crystal pores but adheres at 
the exterior. The group of Schulz-Ekloff and Wohrle started an attempt to accomplish 
rhodamine inclusion in AIPO4-5 crystals, even so the pores are too small. Their idea 
is to modify the dye molecule in a way that makes the dye molecule electrically resem- 
ble a template molecule. Their new derivative rhodamine BE50 (ethanaminium, A^-[6- 
(diethylamino)-9-[2-(iV, iV-dimethyl-3-amino-l-propoxycarbonyl)phenyl]-3H-xanthen-3-y- 
lidene]-A^-ethyl-chloride; cf. Fig. [T^ was synthesized by esterification of rhodamine B 
(Rh B) with 3-dimethylamino-l-propanol |l()3j . It was shown that the concentration of 
rhodamine BE50 (Rh BE50) achievable by crystallization inclusion in AIPO4-5 exceeds 
the possible Rh B concentration by a factor of 3-4. This was attributed to the different 
molecule structures, i.e., to the zwitterionic nature of Rh B on one hand, and to the ad- 
ditional positive charge of a protonated aliphatic amino group of Rh BE50 on the other 
hand |103j . The latter molecules with the localized positive charge are more compatible 
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with the AIPO4-5 framework than the delocahzed charge of Rh B. As a consequence they 
observe that at a given dye concentration Rh BE50 inclusion leads to a better crystal 
morphology than Rh B inclusion |I04j. 





Figure 20: Morphology of typical rhodamine BE5O/AIPO4-5 crystals with lasing prop- 
erties; dye content ca. 0.5 wt-% or 1 dye molecule per 75 unit cells. 

The gel for the synthesis of the AIPO4 crystals is prepared according to recipes jHSl 
I105J which are modified for the purpose of crystallization inclusion of dyes |106j . To 
a suspension of 61.6 mmol AI2O3 as aluminum source |107j and 75 g deionized water, 
61.6 mmol P2O5 |108j in 11.3 g deionized water is added under mechanical stirring. 
After 5 min a uniform gel formed and then 92.4 mmol tripropylamine |109j is added 
slowly. Subsequently, the appropriate amount (0.1-10 mmol) of Rh BE50 dye powder 
is mixed with the gel. The synthesis of the Rh BE5O/AIPO4-5 crystals is performed by 
microwave heating |110j , which proved to be superior in respect of avoiding damage of 
sensitive dyes like coumarines |lllj as well as of reducing the time of synthesis |103j . 
Unlike the pyridine 2 molecules, which with a diameter of 0.6 nm fit into the 0.73 nm wide 
pores of the AIPO4-5 host, the Rh BE50 molecules with dimension of 0.91 x 1.36 nm^ 
are accommodated in defect sites, or mesopores, of the host crystal. Remarkably, up to 
concentrations of 1 molecule per 75 unit cells this remains without any visible negative 
consequences for the crystal morphology, as is documented in Fig. EHl 

In Fig. ^1 the dichroic properties of the compound are illustrated. In comparison with 
the pyridine-2/AlP04-5 compound, the dichroism of the Rh BE50 compound is reduced 
and the fluorescence emission is only partially polarized. This is an indication that the 
electrical anisotropy of the host structure does not fully carry through the mesopores, 
resulting in a weak dye alignment in the statistical average. 

7.2 Microresonator structure 

As is visible in Fig. ^J- and as the polarization of the emitted fluorescence of the com- 
pounds indicates - the absorption, as well as the emission dipole moment of the included 
dyes are oriented preferentially along the crystal c-axis (in the pyridine 2/AIPO4-5 com- 
pound the orientation is complete). As dipole emission along the dipole axis is not 
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possible, the emission parallel to a plane perpendicular to the prevailing dipole orienta- 
tion (i.e. the hexagonal axis) is enhanced. Here a bundle of emission directions meets 
the condition for total internal reflection (TIR) at the hexagonal side faces inside the 
crystal. In a whispering-gallery-mode-like way the corresponding emission can circulate 
sufficiently often to accumulate the gain required to overcome the lasing threshold; cf. 
Fig. 121 in which this intuitive model is illustrated for a particular ray bundle of high 
symmetry. However in section [HI we have shown that for resonators with a sizes in the 
order of a few wavelengths the naive ray picture does not correctly represent the field 
modes. E.g. the ray picture insinuates a mode concentration in the center of the faces 
but field-free corners. This, however, is not consistent with the experimental evidence, 
which clearly shows that the emission occurs at the corners; cf. Fig. ISEl As we already 
pointed out in section \t).l'A\ the main feature that distinguishes the hexagonal resonator 
from other common whispering-gallery type cavities such as microdroplets 112j or semi- 
conductor disk lasers [HZl I113J , is that the latter do not exhibit sharp corners and fiat 
sides. There we described that portions of the boundary in convex resonators act as fo- 
cussing or defocussing elements, whereas the straight sides of a hexagon are neither one 
nor the other. As a result we have shown that ray paths in the hexagon display a degree 
of complexity that cannot be classified as chaotic and was called pseudointegrable. 

The closed ray path underlying Fig. |21 is only one member of an infinite family of 
periodic orbits of the hexagon billiard that all have the same length, L = 3 x WoF. 
Long-lived cavity modes exist only if the corresponding rays satisfy the condition of TIR 
at the interface, sin x > l/""-! where x is the angle of incidence with respect to the surface 
normal. In a naive ray approach one would furthermore obtain the spectrum of modes 
by requiring an integer number of half wavelengths to fit into L, leading to constructive 
interference on a round-trip. As we shall see shortly, this estimate is justified, even 
though a proper treatment of the ray-wave connection has to take into account that any 
given mode is in fact made up of a whole family of different ray paths. The ultimate 
breakdown of the ray model, however, occurs when a ray orbit hits the corners where the 
classical laws of refraction and reflection are no longer deflned. To characterize the size 
of the samples, we can specify either the radius R of the hexagon at the corner points 
or - more conveniently - the width over flats (WoF) characterized hy R = WoF/y/S. 

7.2.1 Wave picture: spectral properties 

Because the hexagonal faces are neither focussing nor defocussing, there is no obvious 
way of determining the weight that should be given to individual members of a ray 
family in order to predict the spatial structure of the resulting mode. Full solutions of 
Maxwell's equations have therefore been carried out for the TM polarized modes of a 
dielectric hexagonal prism, using methods previously applied in jHIl I113J and discussed 
in sectional In anticipation of the experimental spectra discussed in section 17^ we focus 
on three different sample sizes characterized by a WoF of 4.5 /xm, 7.5 /im, and 22 /im. 
The aim is to understand the observed laser line spacings and the emission directionality. 
To document that orbits of different families determine the characteristic mode spac- 
ing of the cavity, Fig.l2T]shows scattering spectra for different sample sizes in the vicinity 
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Figure 21: Calculated scattering intensity spectra of a hexagonal cylinder for plane- 
wave incidence at 15° to a side face and detection at 60° from incidence, (a) 
corresponds to a spectral interval A ~ 653... 8 16 nm for width over fiats 
(WoF) 4.5 /im; (b) covers the interval A ~ 605. . .680 nm for WoF 7.5 /zm. 
Vertical lines are guides to the eye, indicating narrow resonances. The spacing 
between resonances is A{kR) ^ 0.84 in (a) and A{kR) ^ 0.83 in (b), in good 
agreement with the characteristic mode spacing A{kR)c ~ 0.83 of a closed 
hexagonal orbit. Expected resonances not clearly seen in the above spectra 
are marked by dashed lines; they appear at other detection angles. 



of the experimental wavelengths. Intensity is plotted versus dimensionless wavenumber 
kR, where k = 2tt/X. This is the natural scale for comparison with semiclassical predic- 
tions because modes differing by one node along a closed path should then be equally 
spaced, with a characteristic separation A{kR)c = 271 R/{nL) = 27i/{3\/3n) = 0.825 
independent of the sample size. The expected wavelength spacing of the modes (free 
spectral range FSR) is AA = A^ x A{kR)c/{2TTR) ^ 23 nm in (a) and AA ^ 11 nm in 
(b). For WoF = 22yum, we obtain AA ~ 4.9 nm. Figure I?I1 indeed shows a series of 
resonant features with approximately the predicted wavevector spacing. 

Each of the peaks marked in Fig. l^lT b) is in fact a multiplet which is not resolved, 
because the splittings of the individual modes comprising the multiplet are smaller than 
their passive linewidths. There is evidence for this because several of the peaks are 
very asymmetric, and in particular exhibit a steep slope on one side. For an isolated 
resonance, the most general lineshape that could arise is the Fano function (of which 
the Lorentzian is a special case), which however does not yield satisfactory fits here. 

To further expose the multiplet structure, we modeled deviations from the ideal hexag- 
onal shape which could lead to narrower individual linewidths and increase the multiplet 
splitting. Shape perturbations were chosen that preserve the Dq^ point group symmetry 
and hence remove only "accidental" quasi-degeneracies. The actual perturbation that 
is present in the samples of Figs. |21] and |2S1 eluded experimental characterization, so 
that a model calculation can reproduce only generic features which are insensitive to 
the precise type of perturbation. One such feature is the average mode spacing after 
degeneracies have been lifted sufficiently. 

Figure |22la) shows the spectrum of a rounded hexagon where the radius of curvature 
at the corners is p ~ 0.9 A (assuming A ~ 610 nm for definiteness) . No qualitative dif- 
ference to Fig. l^lT b) is seen, except that the resonant features have become somewhat 
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Figure 22: Calculated scattering intensity spectra for slightly rounded hexagonal cavities 
(shapes depicted as insets). The incoming plane wave is at an angle of 15° to 
a facet in (a) and 30° in (b); detection occurs at 60° from incidence. Spacings 
between modes of the same color agree well with A{kR)c ~ 0.83, cf. Fig. |^ 
All resonances in (a) appear as doublets. At kR ^ 42.5 the doublet structure 
is seen most clearly. In (b), stronger deviation from hexagonal shape leads 
to further lifting of degeneracies. Dashed lines mark expected resonances not 
seen at this observation angle. 



narrower, thus enabling us to identify two distinct series of modes with caracteristic 
spacing A{kR)c- This indicates that departures from sharp corners are not resolved in 
the wave equation when their scale is smaller than A. A qualitatively different spec- 
trum is observed in Fig. Elfb) where p ^ 3.7 A. Here, the perturbation reveals three 
well-separated, interpenetrating combs of modes, again with period A{kR)c- There are 
21 distinct resonances in the wavelength interval of Fig. Elfb), which translates to an 
average mode spacing of AA ~ 3.6 nm for a WoF = 7.5 /im resonator, well in agreement 
with the experiment; cf. Fig. |2S1 

In order to verify that no further modes will be revealed by other choices of deforma- 
tion, an independent estimate of the average density of modes can be made based on 
semiclassical considerations 
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Here, dN is the number of modes in the interval d{kR). The result is ( '^^ 
hence we expect ~ 22 modes in the interval of Fig 
the actual count. 
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d(kR)i -4.6, and 
b), again in good agreement with 



7.2.2 Wave picture: intensity profile 

There is one class of quasi-degeneracies that is not removed by any of the perturbations 
in Fig. 1221 Their physical origin is time reversal symmetry for the ray motion inside the 
cavity. Any of the periodic orbits can be traversed either clockwise or counterclockwise. 
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and the same holds for more general ray paths. The different propagation directions can 
be linearly combined in various ways to obtain nearly-degenerate standing-wave patterns 
that differ only in their parity with respect to some of the crystal's reflection axes. A 
minute splitting does exist because the nonintegrability of the ray motion implies that 
the propagation direction itself is not a "good quantum number" , i.e. reversals of the 
sense of rotation are unlikely but not impossible in the wave equation. This is analogous 
to quantum tunneling and hence leads only to exponentially small splittings that can 
be neglected on the scale of the individual resonance linewidths |115j . These multiplets 
have been counted as one resonance in ()367|) . 




Figure 23: False-color representation of the cross- sectional intensity in the ideal hexagon 
for a mode with (a) kR = 22.89 [cf. Fig.|ni(a)] and (b) kR = 42.78 [cf. Fig.lHl 
(b)]. The resonance width is 6{kR) = 0.10 in (a) and 6{kR) = 0.04 in (b). 



Following this reasoning, in Fig. |221 the traveling-wave patterns belonging to one of 
the resonances in Fig. l^lT a) and (b), respectively, is plotted. High intensity ridges inside 
the resonator form a whispering-gallery-like pattern that decays from the interface into 
the cavity center. The number of ridges in the radial direction (perpendicular to a side 
face) provides an approximate analogue of a transverse mode order, however upon closer 
examination one sees that the number of ridges and nodal lines is not uniquely defined, 
in particular along a diameter joining opposite corners. The modes can therefore not be 
properly labeled by "good quantum numbers" characterizing the number of radial and 
azimuthal nodes - this is a direct consequence of the nonintegrability of the problem. 
The most significant difference to the whispering-gallery modes of a circular cavity is 
clearly the anisotropic emission. High intensity is seen to emanate predominantly from 
the corners and is directed almost parallel to an adjacent crystal facet. The overall 
emission pattern is very similar in both modes despite the large difference in size (or 
kR) between the two hexagons. 
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7.3 Laser properties 



The dye loaded molecular sieve microcrystals presented in section [7.1.1I and 17.1.21 were 
excited with 10 ns pulses from the 532 nm second harmonic of a Nd:YAG-laser, and 
the emitted luminescence was collected with a 42° aperture lens relaying the microlaser 
emission to a spectrometer and an imaging system consisting of a cooled low noise 
CCD-camera. 



7.3.1 Pyridine2 AIP04-5 compound 



In Fig. 123] the emission and lasing spectra of three pyridine 2-loaded compounds with 
different dye content are compared. Fig. H^ represents the class of regularly shaped 
samples with a low dye content around 0.1 wt% [cf. also Fig. ITTT left)]. In this class the 
fluorescence emission maximum was observed between 645 nm to 665 nm, where the 
shift to longer wavelength correlates with the increase of dye content that was assessed 
by the saturation of the red color. In none of these rod-shaped crystals laser emission 
was observed. 




'- ^.ft,y£LENl3TH J nm 

Figure 24: Emission and lasing spectra of different pyridine 2-loaded compounds with 

from a to c increasing amounts of included dye. The dye concentration was 
estimated empirically based on the sample color depth. The width over flats 
of sample b is 4 /im. The free spectral range (FSR) of this resonator is so 
large (ca. 25 nm) that only one emission mode acquires the available gain, 
resulting in single line emission. On the other hand, sample c is ca. 5 times 
larger. In this case the FSR is around 5 nm so that laser emission can occur 
on a multitude of modes simultaneously. 

However, narrow laser emission peaks were observed in fascicled samples with a de- 
tected linewidth of ca. 0.3 nm which was limited by the spectrometer resolution. Here 
the emission maxima were observed at wavelenghts up to 695 nm, and again, increasing 
dye concentration correlated with increasing redshift (cf. Figs. |^Jd and c). As already 
mentioned, together with the increasing dye content the crystal morphology becomes 
more and more disturbed. The observation of a disturbed morphology and the red-shift 
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of the emission spectrum are consistent with the hypothesis of a host-guest interaction 
which increases with dye content. In fact both, the pyridine 2 molecules, as well as 
the AIPO4-5 framework, carry a static dipole moment [HHlinni- As a consequence, the 
increased buildup of electrostatic energy in the crystal lattice has to be compensated 
by an increasing amount of stacking faults. On the other side, the mechanism of the 
redshift is not unequivocally identified, yet, but is probably related to the one discussed 
in [TTHllTTTj . 
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Figure 25: Lasing spectrum of a rhodamine BE50-loaded AIPO4-5 microcrystal with a 
concentration around 75 unit cells per dye molecule and size of 7.5 jum width 
over fiats. 



7.3.2 Rhodamine BE50 AIP04-5 compound 

Also with these samples the same correlation between the emission wavelength and dye 
concentration was observed; cf. also |lU3j . In contrast to the pyridine 2-loaded samples, 
the fluorescence emission was not completely polarized. The observed polarization con- 
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^p — J ,j was around 10%, indicating that in the average the Rh BE50 molecules 
are only weakly aligned with respect to the host crystal. Laser emission was observed in 
samples with a dye concentration around 75 unit cells per dye molecule, corresponding 
to 0.5 wt%; cf. Fig. ESI 



7.3.3 Laser properties of molecular sieve dye compounds 

Independent of the type of loading, in most microcrystals with WoF ^ 8 yum lasing 
was observed to occur on several sharp lines with instrument resolution limited width (a 
typical example is shown in Fig. 03j) . In general the laser emission lines are not equally 
spaced. In fact, the free spectral range [FSR) of 11 nm corresponding to the resonator 
size of WoF = 8 /im is far above the line spacing of 3.2, 4.3 and 3.4 nm shown in 
Fig. 123 This is in agreement with the theoretical model of paragraph 17.2.11 in which 
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the average lasing mode spacing (after lifting the quasi- degeneracies in the ideal hexagon) 
was estimated to be AA = 3.6 nm. Also in agreement with the theoretical discussion 
are the regions where the laser light leaves the hexagonal resonators. Figure I^BT left) 
shows the laser emission as bright spots. Clearly the emission is concentrated along the 
crystal edges. The complex emission distribution is compatible with the simultaneously 
recorded spectrum (cf. Fig. E3|) which reveals emission on four modes. 

While the average line spacing of 3.6 nm observed in the sample of Fig. EHl is not 
compatible with the free spectral range of AA = 12 nm of the corresponding resonator, 
the 4.2 nm spacing of the 3 dominant peaks in the sample shown in Fig. 1^^ is in accord 
with the FSR resulting from the 22-/im-WoF hexagonal resonator. On the other hand, 
the theoretical model for the 22/im-WoF hexagonal resonator (cf. Eq. ()367|l ) yields an 
average mode spacing of AA ~ 0.5 nm which is close to the spectrometer resolution. This 
high spectral density explains the large background in the lasing spectrum of Fig. |211 It 
is likely that not all the individual lasing modes in this sample are resolved, and hence 
part of the shoulder on which the three peaks of curve c sit is probably shaped by a 
series of closely spaced lasing modes. 

On the other hand, samples with smaller resonator [WoF ^ 4 /im), as e.g. the one 
shown in Fig. 1^^ . emitted one single laser line. Thus the emission is unadulterated by 
hole burning induced multimode beating and interference, and this results in the simple 
emission pattern shown in Fig I^UT right), where two ca. 1 /xm-spots (~ microscope res- 
olution limit) mark the region of laser emission, which, again, is located at the crystal 
edges. Remarkably, compared to larger samples, the ratio of the line peak to the un- 
derlying fluorescence shoulder of these small lasers is an order of magnitude higher (cf. 
Figs. El and 123). 




Figure 26: Patterns of the laser emission show that the emission originates from regions 
along side edges. Left: rhodamine BE50/A1PO4-5 compound; width over 
flats 7.5 fim. An electron micrograph of the sample is shown in Fig. QUI with 
horizontal c-axis. Here the c-axis orientation is nearly vertical. The corre- 
sponding emission spectrum is shown in Fig. 123 Right: pyridine 2/AIPO4-5 
compound; width over flats 4.5 yum. The corresponding sample and emission 
spectrum is represented in Fig. |^ b. 
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7.3.4 Laser threshold 

Figure 123 illustrates the differential efficiency behavior of a typical microlaser with 
WoF < lOyum and one with WoF > 10 /im. Lasing threshold for the latter size samples 
was around 0.5 MW/cm^, regardless of the type of dye loading. On the other side, 
crystals of smaller size {WoF = 4 /im) from the same synthesis batch revealed a con- 
siderably smaller threshold (0.12 MW/cm^) and a factor of > 7 larger differential gain. 
We assume that this is a consequence of the quantum size effects which we described in 
section El 

It is informative to compare the threshold of molecular sieve microlasers with vertical 
cavity surface emitting lasers (VCSELs). For that we convert the pump irradiance / 
(incident optical power per unit area) into a flux of photons (p (number of pump photons 
per second incident on the mode cross section A) 

, lA lAXr, , , 

= — = — -^ . (368) 

nu he 

If 1] denotes the efficiency with which an incident pump photon actually contributes to 
an excitation, then (p x rj corresponds to the number of quantum processes per second 
required to reach a given level of inversion. Multiplication of this number with the 
charge of an electron q^ = 1.6 x 10~^^As gives the electric current which corresponds to 
the pump current flowing in a comparable electric device, such as a VCSEL. According 
to ()368|) the threshold power density of 0.12 MW/cm^ incident on the surface of the 
molecular sieve laser shown in Fig. |2D3 of 1 x 4 /im (cf. Fig. OHb) corresponds to a flux 
of 1.4 X 10^^ s~^ 532-nm-photons. If we take into account that the pump field was not 
polarized, but the dye molecules are aligned, then only 50% of the pump photons could 
contribute to the inversion. In addition the absorption length for 532 nni radiation is 
longer than the crystal size so that only a fraction of the incident pump is absorbed 
by the dye molecules. If we assume that 50% of the pump are absorbed, then the 
effects of polarization and absorbtion length together result in a quantum efficiency of 
7] = 0.25. With this quantum efficiency we obtain a comparable electric threshold current 
of 560 fiA, which compares well with the threshold current of VCSELs of comparable 
size. Thus, in terms of elementary (quantum) pump processes needed to reach lasing 
threshold the molecular sieve lasers are as good as actual VCSELs. 

7.4 Photostability 

Photostability is usually a critical issue with organic dyes. Corresponding investigations 
were carried out with pyridine 2, as well as rhodamine BE50 loaded samples, revealing 
some unexpected results. 

7.4.1 Photostability of pyridine 2 compounds 

The photostability of pyridine 2 compounds was studied with samples exhibiting an 
undisturbed morphology, similar to the one shown in Fig. ^| The samples were irradi- 
ated with 10 Hz trains of 10 ns pulses of the 532 nm second harmonic of a Nd:YAG-laser 
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Figure 27: Lasing threshold and differential efficiency of typical AlP04-5/dye com- 
pounds. Shown is the peak of the laser emission spectrum normalized by 
the fluorescence shoulder as a function of the pump power density for the 
sample shown in Figs. Elb and | 



and a power density of 5 MW/cm^ . Figure |2H1 illustrates the wane of fluorescence ac- 
tivity of a pyridine 2-loaded AIPO4-5 sample under such bleaching irradiation. After 
a bleaching period of 140 seconds the exposure was interrupted for 18 minutes. Then 
the bleaching procedure was resumed. The figure shows that the fluorescence recovers 
during the intermission. 

The physical origin of this unexpected fluorescence recovery is not clear yet. If we 
assume that bleaching consists in breaking bonds of the dye molecules, then we have 
to consider bond energies in the eV-range. Even if we assume that the dye debris stay 
encaged in their pores, spontaneous or thermally activated self healing of broken eV- 
bonds seems not very probable. Considering the stereometrically restricted possibilities 
inside the molecular sieve framework together with diffusion distances of several fun, 
and the observed recovery time in the range of minutes, we believe that diffusion of 
new, intact dye molecules into the irradiated crystal volume is more plausible than self 
healing. 

As bleaching reduces the concentration of dye, a blueshift of the fluorescence is ex- 
pected with increasing photobleaching [HHl ESI • However, the 656- nm- fluorescence emis- 
sion maximum of this sample is already at the shortest observed wavelength (cf. Fig.l24j). 
corresponding to a low dye concentration, and consequently, to weak dipole interactions. 
Thus, the blueshift under these circumstances must be rather small. This explains that 
a blueshift was not observed with these pyridine 2-loaded samples. 

7.4.2 Photostability of rhodamine BE50 compounds 

The rhodamine BE50-loaded samples under investigation contained dye at a concentra- 
tion of around one Rh BE50 molecule per 75 unit cells, and therefore bleaching caused a 
detectable 4 nm shift of the fluorescence towards the blue. The 532 nm bleach irradiance 
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Figure 28: Fluorescence activity of a pyridine 2-loaded AIPO4-5 sample under bleaching 
laser irradiation. Left: Micrograph of the sample crystal. Center: Fluores- 
cence activity as a function of time: After a first bleach period of 140 seconds 
the bleach beam is interrupted for 18 minutes. During this intermission the 
fluorescence recovers to start the second bleach period with 3x stronger emis- 
sion. Right: The bleaching laser is incident from the left and concentrated 
in the center of the crystal. Shown is the fluorescence distribution at the end 
of the first bleach period of 140 s. Clearly visible is the bleached hole in the 
center, where the bleach beam was concentrated. 



was 0.5 MW/cm^ with this samples. In contrast to the pyridine 2 samples these crystals 
exhibited laser emission. Observing the laser emission spectrum while bleaching the 
samples, a further consequence of the blueshift was revealed: Blueshift of the fluores- 
cence reduces the overlap of the fluorescence band with the absorption spectrum, and as 
a result, laser modes at lower wavelengths will suffer less losses with increasing bleaching. 
This is documented in Fig. 123 where the intensity of mode b with the shortest oscil- 
lating wavelength increases, while longer-wavelength-modes a and c decrease during the 
bleach procedure. At the same time a 0.2 nm blueshift of the oscillation wavelength was 
detected. We attribute this to a weak decrease of the refractive index of the resonator 
material due to the smaller polarizability of the dye debris. In contrast to the pyridine 
2 samples, however, recovery of the fluorescence was not detected. As we mentioned in 
section 17. 1.21 the Rh 50BE molecules are considerably larger than pyridine 2 molecules, 
and are encaged in mesopores. As a result their mobility in the molecular sieve frame- 
work is severly hampered. So, diffusion of intact molecules into the bleached volume 
occurs - if ever - on timescales considerably larger than minutes. This is probably the 
reason why fluorescence recovery was not observed in rhodamine BE50 compounds. 



8 Conclusion 

This chapter was centered on microscopic lasers realized with nanoporous materials, 
especially molecular sieve dye compounds. We discussed the effects of a microresonator 
on spontaneous and stimulated emission properties, and we showed that wavelength 
scale laser resonators may exhibit lasing without threshold. Such effects are not just 
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Figure 29: Effect of photobleaching on the laser emission intensity: short wavelength 
line b grows with progressive bleaching, while longer wavelength lines a and 
c decrease. 

thoughtful observations, but were investigated in real devices as well, and we reviewed 
observations of threshold reduction in molecular sieve lasers. 

The microresonators of the molecular sieve lasers which were presented up to now 
exhibit a hexagonal resonator geometry. Hexagonal microresonators in which the light 
field is confined by total internal reflection at the dielectric boundary define a new class 
of pseudointegrable optical structure of which we reviewed the properties in some detail. 

As mentioned, molecular sieve microlasers can be fabricated in large amounts with 
the same large scale processes which are used to produce the molecular sieve crystallites 
in the petrochemical industry. Although applications of such laser powders have not 
been reported yet, many practical uses can be envisaged, such as for example in efficient 
fluorescent, or even lasing paint pigments or lasing pixels. 
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